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Abstract 

We show the existence of global-in-time weak solutions to a general class of coupled 
Hookean-type bead-spring chain models that arise from the kinetic theory of dilute solu- 
tions of polymeric liquids with noninteracting polymer chains. The class of models involves 
the unsteady incompressible Navier-Stokes equations in a bounded domain in R"*, d = 2 
or 3, for the velocity and the pressure of the fluid, with an elastic extra-stress tensor ap- 
pearing on the right-hand side in the momentum equation. The extra-stress tensor stems 
from the random movement of the polymer chains and is defined by the Kramers expres- 
sion through the associated probability density function that satisfies a Fokker-Planck-type 
parabolic equation, a crucial feature of which is the presence of a center-of-mass diffusion 
term. We require no structural assumptions on the drag term in the Fokker-Planck equa- 
tion; in particular, the drag term need not be corotational. With a square-integrable and 
divergence-free initial velocity datum uq for the Navier-Stokes equation and a nonnegative 
initial probability density function ipo for the Fokker-Planck equation, which has finite rel- 
ative entropy with respect to the Maxwellian M, we prove the existence of a global-in-time 
weak solution t {u(t),%l>{t)) to the coupled Navier-Stokes-Fokker-Planck system, satisfying 
the initial condition («(0),^(0)) = (uo,^o), such that t h->- u{t) belongs to the classical Leray 
space and t i— >■ xp{t) has bounded relative entropy with respect to M and t h->- ^(t)/M has 
integrable Fisher information (w.r.t. the measure du := M[q)Aq Ax) over any time interval 
[0,r], r > 0. If the density of body forces / on the right-hand side of the Navier-Stokes 
momentum equation vanishes, then t ^ {u{t),%j)(t)) decays exponentially in time to (0, M) in 
the L^ X L^ norm, at a rate that is independent of (uo, V'o) and of the centre-of-mass diffusion 
coefficient. 

An abbreviated version of this paper has been submitted for publication m Mathematical Models 
and Methods in Applied Sciences (M3AS). 

Keywords: Kinetic polymer models, Hookean, Rouse chain, Navier-Stokes-Fokker-Planck 
system. 



1 Introduction 

This paper establishes the existence of global-in-time weak solutions to a large class of bead-spring 
chain models with Hookean-type spring potentials, — a system of nonlinear partial differential 
equations that arises from the kinetic theory of dilute polymer solutions. The solvent is an in- 
compressible, viscous, isothermal Newtonian fluid confined to a bounded open Lipschitz domain 
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il C M'^, d = 2 or 3, with boundary dil. For the sake of simpHcity of presentation, we shall 
suppose that D, has 'solid' boundary dfl; the velocity field u will then satisfy the no-slip boundary 
condition w = on dfl. The polymer chains, which are suspended in the solvent, are assumed not 
to interact with each other. The conservation of momentum and mass equations for the solvent 
then have the form of the incompressible Navier-Stokes equations in which the elastic extra-stress 
tensor r (i.e., the polymeric part of the Cauchy stress tensor) appears as a source term: 

Given T e R>o, find u : {x,t) e H x [0,T] u{x,t) e R'^ and p : {x,t) £ n x {0,T] 
p{x, t) gM. such that 

du 

-7 + (u-V^)u-:/A^M + V^p = / + V^-T innx{0,T], (1.1a) 

at ~ ~ ~ ~~ ^ ^ „ 

V^-u = mnx{0,T], (1.1b) 

u^O ondnx{0,T], (1.1c) 

u{x,0) = uo{x) Wxen. (l.ld) 

It is assumed that each of the equations above has been written in its nondimensional form and 
the momentum equation (jl.lap has been normalized so that the Strouhal number is equal to 1; 
u denotes a nondimensional velocity, defined as the velocity field scaled by the characteristic flow 
speed Uo] V € K>o is the reciprocal of the Reynolds number, i.e. the ratio of the kinematic 
viscosity coefficient of the solvent and LqUq, where Lq is a characteristic length-scale of the flow; 
p is the product of the nondimensional pressure and the Euler number; and / is the product of 
the nondimensional density of body forces and the Richardson number. 

In a bead-spring chain model, consisting oi K -\- 1 beads coupled with K elastic springs to 
represent a polymer chain, the extra-stress tensor r is defined by the Kramers expression as 
a weighted average of V'l the probability density function of the (random) conformation vector 
q = {qi, . . . ,qK)'^ G R^'^ of the chain (cf. (jl.Sp below), with qi representing the d-component 
conformation/orientation vector of the ith spring. The Kolmogorov equation satisfied by ■0 is a 
second-order parabolic equation, the Fokker-Planck equation, whose transport coefficients depend 
on the velocity field u. The domain D of admissible conformation vectors D C M^'* is a K-fo\d 
Cartesian product Di x ■ ■ ■ x Dk of balanced convex open sets Di C W^, i = 1, . . . , if; the term 
balanced means that qi G Di if, and only if, — e Di. Hence, in particular, G Di, i = 1, . . . , K . 

Typically, Di is the whole of or a bounded open c?-dimensional ball centred at the origin G M'' 
for each i = 1, . . . , K . When K — 1, the model is referred to as the dumbbell model. 

Let Oi C [0,oo) denote the image of Di under the mapping qi G Di i— ^ ^|<Zjp, and consider 

the spring-potential Ui G W^^^ {Oi;'R>o), i = 1,. . .,K. Clearly, G Oi. We shall suppose that 
Ui{0) = and that Ui is monotonic increasing and unbounded on Oi for each i = 1, . . . , K . The 
elastic spring-force Fi : Di C M.'^ ^ of the ith spring in the chain is defined by 

F^{q^)^U^{^\qJ^)q^, l^l,...,K. (1.2) 

Example 1.1 In the Hookean dumbbell model K — 1, and the spring force is defined by F{q) — q, 
with q G D = R'^, corresponding to U{s) — s, s E O — [0, oo). o 

Unfortunately, we are not able to deal with the pure Hookean model. The compactness ar- 
gument that forms the core of our current existence proof requires that the potentials Ui have 
superlinear growth at infinity. For example for any Soo > and i? > 1, we can deal with a 
potential of the form 



which approximates the Hookean potential U{s) — s 



for s G [0,Soo], 

for S> SrycA 
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We shall assume in what follows that D is a Cartesian product oi Di = R'^, i — 1, . . . ,K, with 
K > 1. We shall further suppose that ioT i = 1, . . . , K there exist constants Cij > 0, j — 1, 2, 3, 4, 
such that the (normalized) Maxwellian Mi, defined by 



e ^ ~ dqi 



Di 



and the associated monotonically increasing potentials Ui G W^^^{[0, oo); R>o) satisfy, for a t9 > 1, 

U^il\q^\^) = Cai^\q^\Y as|g,|^oo, (1.4a) 

C/;(i|'Z.P) < + c,3 C^\q^\^r-' Vg, e A, (1.4b) 

and hence 

-=.1(119.1')" 

Mi{qi) = Cj4e ~ as \qi\ oo. (1.4c) 

Hence our use of the words Hookean-type model throughout the paper (instead of Hookean model, 
which would have corresponded to taking = 1 in the above). 
The Maxwellian in the model is then defined by 

K 

M{q):^\lA'U{qr) "^q '.^ {qi, . . . ,qK) e D -.^ Di x ■ ■ ■ x Dk- (1-5) 

i=l ~ ~ ~ 

Observe that, for i = 1, . . . , ii', 

M{qjY,AM{qX' = -[M(g)]-ly,,M(q) = V M^g^l^) = U^i^Qj') (1-6) 
Since [U^{^\q^\^)]^ = {-logM,{q,) + Const. )2, it follows from (fTibl c) that 

'l + [C/,(i|g,|2)]2 + [[/;( i|5,|2)]2]M,(g,)dg,<oo, i^l,...,K. (1.7) 

The governing equations of the general class of Hookean-type chain models with centre-of-mass 
diffusion are ((TTlal^d), where the extra-stress tensor r is defined by the Kramers expression: 

Tix, t)^k(j2j^ ^fe: 2' ^) 1^ 1^ (^Ig^l") - Pix, t) l)j , (1.8) 
with the density of polymer chains located at x at time t given by 

p{x,t) = / V(s,g,Odg. (1.9) 

J D ~ ~ 

The probability density function ?/> is a solution of the Fokker-Planck equation 

K 

-I- (ii ■ xy^'ii/i -I- 

dt 



dtp 



= eA,^ + ^ ^^A,,V,. . (^AfV,^.(^^^^ mnxDx{0,T], (1.10) 

with q{v) = YxV, where (ya;w)(j;,i) G M'*^'* and {V^ — The dimensionless constant 

k > featuring in (jl.Sp is a constant multiple of the product of the Boltzmann constant ks and 
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the absolute temperature T. In (|1.10p . e > is the centre-of-mass diffusion coefficient defined 
as e := {£q / Lq)'^ / {A{K + 1)A) with £o := y/kBT/H signifying the characteristic microscopic 
length-scale and A := {(/4H){Uo/ Lq), where C > is a friction coefficient and 77 > is a spring- 
constant. The dimensionless parameter A e K.>o, called the Weissenberg number (and usually 
denoted by Wi), characterizes the elastic relaxation property of the fluid, and A = {Aij)fj^^ 
is a symmetric positive definite matrix, the Rouse matrix, or connectivity matrix; for example, 
A ~ tridiag[— 1. 2, —1] in the case of a linear chain; see, for example, Nitta [33] . 

Definition 1.1 The collection of equations and structural hypotheses (|l.la| -d) - (|1.10p will be re- 
ferred to throughout the paper as model (Pe), or as Hookean-type (bead-spring chain) models with 
centre-of-mass diffusion. 

A noteworthy feature of equation (|1.10p in the model (Pe) compared to classical Fokker-Planck 
equations for bead-spring models in the literature is the presence of the a;-dissipative centre-of-mass 
diffusion term £ A^ip on the right-hand side of the Fokker-Planck equation (|1.10p . We refer to 
Barrett & Siili ^ for the derivation of (|1.10p in the case oi K = 1; see also the article by Schieber 
[37] concerning generalized dumbbell models with centre-of-mass diffusion, and the recent paper 
of Degond & Liu [M] for a careful justification of the presence of the centre-of-mass diffusion term 
through asymptotic analysis. In standard derivations of bead-spring models the centre-of-mass 
diffusion term is routinely omitted on the grounds that it is several orders of magnitude smaller 
than the other terms in the equation. Indeed, when the characteristic macroscopic length-scale 
Lq k 1, (for example, Lq = diam(f7)), Bhave, Armstrong & Brown TTI estimate the ratio (q/Lq 
to be in the range of about 10^^ to 10~^. However, the omission of the term e A^ip from (|1.10p in 
the case of a heterogeneous solvent velocity u{x,t) is a mathematically counterproductive model 
reduction. When e Axip is absent, ()1.10p becomes a degenerate parabolic equation exhibiting 
hyperbolic behaviour with respect to {x,t). Since the study of weak solutions to the coupled 
problem requires one to work with velocity fields u that have very limited Sobolev regularity 
(typically u e L°°(0, T; L'^{n))nL'^{0, T; ffj(r2))), one is then forced into the technically unpleasant 
framework of hyperbolically degenerate parabolic equations with rough transport coefficients (cf. 
Ambrosio [2] and DiPerna & Lions [16]). For these reasons, here we shall retain the centre-of-mass 
diffusion term in ()1.10p . In order to emphasize that the positive centre-of-mass diffusion coefficient 
e is not a mathematical artifact but the outcome of the physical derivation of the model, in Section 
[5] and thereafter the variables u and ip have been labelled with the subscript e. 

Following the introductory section in the companion paper [10 , which is concerned with anal- 
ogous questions to the ones considered here in the case of bead-spring chains with FENE-type 
potentials, we continue with a brief literature survey. Unless otherwise stated, the center-of-mass 
diffusion term is absent from the model considered in the cited reference (i.e. e is set to 0); also, in 
all references cited, except [TU], K = 1, i.e. a simple dumbbell model is considered, with a single 
spring and a pair of beads, rather than a general bead-spring chain model. 

An early contribution to the existence and uniqueness of local-in-time solutions to a family of 
dumbbell type polymeric flow models is due to Renardy [36 . While the class of potentials F{q) 
considered by Renardy 06] (cf. hypotheses (F) and (F') on pp. 314-315) does include the case 
of a Hookean dumbbell, it excludes the practically relevant case of the FENE dumbbell model. 
More recently, E, Li & Zhang [19] and Li, Zhang & Zhang [27j have revisited the question of local 
existence of solutions for dumbbell models. A further development in this direction is the work of 
Zhang & Zhang 42 , where the local existence of regular solutions to FENE-type dumbbell models 
has been shown. All of these papers require high regularity of the initial data. Constantin }13j 
considered the Navier-Stokes equations coupled to nonlinear Fokker-Planck equations describing 
the evolution of the probability distribution of the particles interacting with the fluid. Otto & 
Tzavaras [35] investigated the Doi model (which is similar to a Hookean model (cf. Example 11.11 
above), except that D — S^) for suspensions of rod-like molecules in the dilute regime. Jourdain, 
Lelievre & Le Bris [24 studied the existence of solutions to the FENE dumbbell model in the case 
of a simple Couette flow. By using tools from the theory of stochastic differential equations, they 
established the existence of a unique local-in-time solution to the FENE dumbbell model in two 
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space dimensions {d = 2) when the velocity field u is unidirectional and of the particular form 

In the case of Hookean dumbbells {K — 1), and assuming e = 0, the coupled microscopic- 
macroscopic model described above yields, formally, taking the second moment of g i— > ip{q^x,t)^ 
the fully macroscopic, Oldroyd-B model of viscoelastic flow. Lions & Masmoudi [29] have shown 
the existence of global-in-time weak solutions to the Oldroyd-B model in a simplified corotational 
setting (i.e. with cr(u) — YxH replaced by ^{YxU ~ (Yxu)"^)) by exploiting the propagation in 
time of the compactness of the solution (i.e. the property that if one takes a sequence of weak 
solutions that converges weakly and such that the corresponding sequence of initial data converges 
strongly, then the weak limit is also a solution) and the DiPerna-Lions [16] theory of renormalized 
solutions to linear hyperbolic equations with nonsmooth transport coefficients. It is not known 
if an identical global existence result for the Oldroyd-B model also holds in the absence of the 
crucial assumption that the drag term is corotational. We note in passing that with e > the 
coupled microscopic-macroscopic model above yields, taking the appropriate moments in the case 
of Hookean dumbbells, a dissipative version of the Oldroyd-B model. In this sense, the Hookean 
dumbbell model has a macroscopic closure: it is the Oldroyd-B model when e = 0, and a dissipative 
version of Oldroyd-B when e > (cf. Barrett & Siili [5^). Barrett & Boyaval j6j have proved a 
global existence result for this dissipative Oldroyd-B model in two space dimensions. In contrast, 
the FENE model is not known to have an exact closure at the macroscopic level, though Du, Yu 
& Liu [T7: and Yu, Du & Liu [4T] have recently considered the analysis of approximate closures of 
the FENE dumbbell model. Lions & Masmoudi ^30^ proved the global existence of weak solutions 
for the corotational FENE dumbbell model, once again corresponding to the case of e = and 
K = 1, and the Doi model, also called the rod model. As in Lions & Masmoudi [29], their proof 
is based on propagation of compactness; see also the related paper of Masmoudi [21] . Recently, 
Masmoudi [32] has extended this analysis to the noncorotational case. 

Previously, El-Kareh & Leal [20] had proposed a steady macroscopic model, with added dissi- 
pation in the equation satisfied by the conformation tensor, defined as 



in order to account for Brownian motion across streamlines; the model can be thought of as an 
approximate macroscopic closure of a FENE-type micro-macro model with centre-of-mass diffu- 
sion. 

Barrett, Schwab & Siili [7j established the existence of, global- in-time, weak solutions to the 
coupled microscopic-macroscopic model ([TTT^d) and (|1.10p with e = 0, if = 1, an a;- mollified ve- 
locity gradient in the Fokker-Planck equation and an j-moUified probability density function ip in 
the Kramers expression — admitting a large class of potentials U (including the Hookean dumbbell 
model as well as general FENE-type dumbbell models); in addition to these mollifications, u in 
the x-convective term {u ■ Yx)tp in the Fokker-Planck equation was also mollified. Unlike Lions & 
Masmoudi [29], the arguments in Barrett, Schwab & Siili [7] did not require the assumption that 
the drag term • {a{u) q ip) in the Fokker-Planck was corotational in the FENE case. 

In Barrett & Siili [8] , we derived the coupled Navier-Stokes-Fokker-Planck model with centre- 
of-mass diffusion stated above, in the case oi K — 1. The anisotropic Friedrichs mollifiers, which 
naturally arise in the derivation of the model in the Kramers expression for the extra-stress tensor 
and in the drag term in the Fokker-Planck equation, were replaced by isotropic Friedrichs molli- 
fiers. We established the existence of global-in-time weak solutions to the model for a general class 
of spring-force-potentials including in particular the FENE potential. We justified also, through a 
rigorous limiting process, certain classical reductions of this model appearing in the literature that 
exclude the centre-of-mass diffusion term from the Fokker-Planck equation on the grounds that 
the diffusion coefficient is small relative to other coefficients featuring in the equation. In the case 
of a corotational drag term we performed a rigorous passage to the limit as the Helmholtz-Stokes 
mollifiers in the Kramers expression and the drag term converge to identity operators. 

In Barrett & Siili '9' we showed the existence of global-in-time weak solutions to general 
noncorotational FENE-type dumbbell models (including the standard FENE dumbbell model) 
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with centre-of-mass diffusion, in the case oi K ^ 1 , with microsropic cut-off in the drag term 



(1.11) 



In Barrett & Siih [10] we took that analysis further by showing the existence of global-in-time 
weak solutions to general noncorotational FENE-type dumbbell models (including the standard 
FENE dumbbell model) with centre-of-mass diffusion, in the general case K > 1, without cut-off 
or mollification. The weak solution was shown to satisfy an energy inequality, and in the absence 
of body forces it was shown to converge exponentially to the equilibrium solution of the problem. 

The present paper extends the analysis in Barrett & Siili to Hookean-type bead-spring 
models with centre-of-mass diffusion, in the general case K > 1, without cut-off or mollification. 
As was noted above our current analysis rules out the possibility of taking d = 1 in (jl.4al -c) , since 
we require that the potentials Ui have superlinear growth as \q\ — > oo. Since the argument is long 
and technical, we give a brief overview of the main steps of the proof here. 

Step 1. Following the approach in Barrett & Siili [5] and motivated by recent papers of 
Jourdain, Lelievre, Le Bris & Otto [25] and Lin, Liu & Zhang [28] (see also Arnold, Markowich, 
Toscani & Unterreiter [S], and Desvillettes & Villani [TS]) concerning the convergence of the 
probability density function to its equilibrium value ipooix,q) := M{q) (corresponding to the 
equilibrium value Uooix) := of the velocity field) in the absence of body forces /, we observe 
that if i/j/AI is bounded above then, for L E M>o sufficiently large, the drag term (jl.lip is equal 
to 

Yq- z{u)qMI3^^ 



(1.12) 



where £ C(K) is a cut-off function defined as 



/3^(s) := min(s,i) = 



s 
L 



for s < L, 
for L < s. 



(1.13) 



More generally, in the case of i^T > 1 , in analogy with (|1.12p , the drag term with cut-off is defined 
by 



K 

i=l 



a{u)q^Mp^ 



It then follows that, for L >• 1, any solution ip of (jl.lOp . such that ^/M is bounded above, also 
satisfies 



'dt 



+ (yy,)V'-H^y,, • (^a(u)g,M/3^(^^^^ 



1 



K K 



i=i j=i 

We impose the following decay/boundary and initial conditions: 



in n X D X (0,T]. 



(1.14) 



M 



1 



K 

^ A,; Vo 



a{u)qil3^ 



on fl X 



M 

K 

X 



eVxi^ • n = 



ondnxDx (0,T] 



V'(-,-,0) = Af(.)/3^ (V^o(-,-)/M(-))>0 



as \qi 



oo 



Dj X (0,r], for i = 1, 



on f2 X _D, 



(1.15a) 

(1.15b) 
(1.15c) 



6 



where n is the unit outward normal to dfl. 

Here tpo is a nonnegative function defined on 51 x 13, with J^ipoix^Q) = 1 for a.e. x G fi. 
We shall also assume that ipo has finite relative entropy with respect to the Maxwellian M; 
i.e. ^ il^o {x, q) log{ipQ {x, q)/M{q)) dq dx < oo. Obviously, if there exists i > such that 
0<'ijjo <LM, then M/3^(V'o/M) = ipo- Henceforth we shall assume that L > 1. 

Definition 1.2 The coupled problem (|l.la| -d). (fTSl) . (fTO)) . (|1.14p . (|1.15a| -c) wiH &e referred to as 
model (Pe^i), or as the Hookean-type (bead-spring chain) models with centre- of-mass diffusion and 
microscopic cut-off. 

In order to highlight the dependence on e and L, in subsequent sections the solution to (jl.l4p . 
(|1.15a| -c) will be labelled ^e,L- Due to the coupling of (I1.14p to (ll.lap through (|1.8I) . the velocity 
and the pressure will also depend on e and L and we shall therefore denote them in subsequent 
sections by u^^l and Pe,L- 

The cut-off has several attractive properties. We observe that the couple (mqo, V'oo), defined 
by Uoo{x) := and ipaoix, q) := M{q), is still an equilibrium solution of (|l.lal -d) with / = 0, (jl.Sp . 
(Il.yp . (|1.14l) . (ll.lSaf c'l for all L > 0. Thus, unlike the truncation of the (unbounded) potential 
proposed in El-Kareh & Leal [30], the introduction of the cut-off function into the Fokker- 
Planck equation (|1.10p does not alter the equilibrium solution {u^,ipoo) of the original Navier- 
Stokes-Fokker-Planck system. In addition, the boundary conditions for i/j on dft x _D x (0, T] and 
n X dD X {0,T] ensure that 

/ ■>jj{x,q,t) dq — / ip{x,q,0)dq for a.e. x € and a.e. t G IR>o. 
Jd ~ ~ ~ Jd ~ ~ 

Step 2. Ideally, one would like to pass to the limit i — > oo in problem (P^^l) to deduce the 
existence of solutions to (Pg). Unfortunately, such a direct attack at the problem is (except in 
the special case oi d — 2, or in the absence of convection terms from the model,) fraught with 
technical difficulties. Instead, we shall first (semi)discretize problem (Pe,L) by an implicit Euler 
scheme with respect to t, with step size At; this results in a time-discrete version (P^^) of (Pe^l). 
By using Schauder's fixed point theorem, we will show in Section |3] the existence of solutions to 
(PfX). In the course of the proof, for technical reasons, a further cut-off, now from below, is 
required, with a cut-off parameter 6 G (0, 1), which we shall let pass to to complete the proof 
of existence of solutions to (P^^) in the limit of 5 — )■ 0+ (cf. Section |3]). Ultimately, of course, 
our aim is to show existence of weak solutions to the Hookean-type models with centre-of-mass 
diffusion, (P^), and that demands passing to the limits At — > 0+ and i — > oo; this then brings us 
to the next step in our argument. 

Step 3. We shall link the time step At to the cut-off parameter i > 1 by demanding that 
At = o{L^^), as L — > oo, so that the only parameter in the problem (P^^) is the cut-off parameter 
(the centre-of-mass diffusion parameter e being fixed). By using special energy estimates, based 
on testing the Fokker-Planck equation in (P^^) with the derivative of the relative entropy with 
respect to the Maxwellian of the model, we show that jt^^ can be bounded, independent of L. 
Specifically is bounded in the norm of the classical Leray space, independent of L; also, 
the L°° norm in time of the relative entropy of tp^j^ and the L"^ norm in time of the Fisher 

information of ?/;^^ := ip^*j^/M are bounded, independent of L. We then use these i-independent 
bounds on the relative entropj^ and the Fisher information to derive L-independent bounds on the 
time-derivatives of and ip^l^ in very weak, negative-order Sobolev norms. 

Step 4- The collection of i-independent bounds from Step 3 then enables us to extract a 
weakly convergent subsequence of solutions to problem (P^^) as L — > oo. We then apply a 
general compactness result in seminormed sets due to Dubinskii fTS , which furnishes strong con- 
vergence of a subsequence of solutions (y^L j. , V'^l'^, ) to (P^£) with At — o{L^^) as L oo, in 
i^(0, T; L^{fl)) X LP{0, T; L^{fl x D)) for any p > 1. A crucial observation is that the set of func- 
tions with finite Fisher information is not a linear space; therefore, typical Aubin-Lions-Simon 
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type compactness results (see, for example, Simon 138' ) do not work in our context; however, Du- 
binskii's compactness theorem, which applies to seminormed sets in the sense of Dubinskii, does, 
enabling us to pass to the limit with the microscopic cut-off parameter L in the model (P^£), 
with At = o(L^^), as i — > oo, to finally deduce the existence of weak solutions to Hookean-type 
models with centre-of-mass diffusion, (Pe). 

The paper is structured as follows. We begin, in Section [2l by stating (Pg ^); the coupled 
Navier-Stokes-Fokker-Planck system with centre-of-mass difl^usion and microscopic cut-off for a 
general class of Hookean-type spring potentials. In Section[3]we show the existence of solutions to 
the time-discrete problem (P^£). In Section |4] we derive a set of L-independent bounds on 
in the classical Leray space, together with L-independent bounds on the relative entropy of 
and the Fisher information of ■ We then use these L-independent bounds on spatial norms to 

obtain L-independent bounds on very weak norms of time-derivatives of u^\^ and V'i^i • Section [5] 
is concerned with the application of Dubinskii's theorem to our problem; and the extraction of a 
strongly convergent subsequence, which we shall then use in Section [6] to pass to the limit with 
the cut-off parameter L in problem (P^^), with Ai = o(L~^), as L — > cx), to deduce the existence 

of a weak solution {u^jipg := M ipg) to problem (Pe). Finally, in Section [71 we show using a 
logarithmic Sobolev inequality and the Csiszar-Kullback inequality that, when / = 0, global weak 
solutions t M- {ue{t), ipe{t)) thus constructed decay exponentially in time to (0, M), at a rate that 
is independent of the choice of the initial data for the Navier-Stokes and Fokker-Planck equations 
and of e. We shall operate within Maxwellian-weighted Sobolev spaces that provide the natural 
functional-analytic framework for the problem. Our proofs require special density and embedding 
results in these spaces, which are proved in the Appendix. 



2 The polymer model (P 



Let f2 C K'^ be a bounded open set with a Lipschitz-continuous boundary dfl, and suppose that 
the set D := £>! X • • • X Dk of admissible elongation vectors q := {qi, . . . , qk) in (|1.10l) is such 

that A = K'', j 1, ■ • ^K. 

Collecting ()l.lal -d). (|1.8|) . and (|1.10|) . we then consider the following initial-boundary- value 
problem, dependent on the parameter L > I. As has been already emphasized in the Introduction, 
the centre-of-diffusion parameter e > is a physical parameter and is regarded as being fixed, 
although we systematically highlight its presence in the model through our subscript notation. 

(Pe,L) Find Ue,L : (x,t) e n X [0,r] ^ Ue,Lix,t) € and ■■ {x,t) € n X (0,T] 

Pe^L{x,t) € M such that 

du 

at ^ ~ 

inl7x(0,r], (2.1a) 

Yx-Ue,L = inl7x(0,r], (2.1b) 

u,^L = Q on917x(0,T], (2.1c) 

u^l{x,0) = uo{x) Vxefi, (2. Id) 

where v £ R>o is the given viscosity, f{x,t) is the given body force and t{i/j^^l) : {x,t) e 
ri X (0, T] K-s- T(-0e,L)(a;, t) e M''^'^ is the symmetric extra-stress tensor, dependent on a probability 
density function : {x,q,t) £ x D x [0,T] i-^ Ve,L(S;9;0 ^ defined as 



T{i;,,L) ^ k (j2c,{^P,^L)j -kp{i^,^L)i. (2.2) 
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Here k £ 



/ is the unit d x d tensor, 



JD ~ ~ ~ ~ ~ 



and 



(2.3a) 



p{lpe.L)ix,t) 



/ ^e,L{x,q,t)dq. 



The Fokker-Planck equation with microscopic cut-off satisfied by V'e.l is: 

K 



i=l 
K K 



a{us^L)qiM l3 



M 



i=l 3 = 1 



= £ A, Ve,L + Y?. ■ ( ^'^ Y 



in X £) X (0,r]. 



Here, for a given L> \, 13^ e C(K) is defined by (|1.13p . ^(t;) = Vx v, and 

A G R^^^ is symmetric positive definite with smallest eigenvalue ao G 
We impose the following decay/boundary and initial conditions: 

1 ^ 

— Aij Wq 



M 



M 



(T{Ue,L) qi P 



M 



— > as |gi| — )■ oo 



K 



on X 



X x(0,T], 

e\Ixil^e,L ■ n = Q on af7 X £) X (0,r], 

V'e.Ll-, -,0) = M(.)/3^(^o(-, OA^l-)) > on r! X 



1, 



(2.3b) 



(2.4) 



(2.5) 



(2.6a) 



(2.6b) 
(2.6c) 

where n is the unit outward normal to d^l. The boundary conditions for ijje.L on d^l x D x (0,T] 



and the decay conditions for tp^^L on fl x {^xf^ij^^ Djj x (0, T] as \qi\ —i'00,i = l,...,K, have 
been chosen so as to ensure that 



D 



^pe,L{x,q,t)dq ^ / ipe,L{x,q,0)dq V(x,t)eOT- 



(2.7) 



D 



Henceforth, we shall write V'e.L ■= '4'e,L/M, -tpo ■= ipo/M. Thus, for example, p.6cp in terms of 
this compact notation becomes: V'6,l(-, •, 0) = /3^(V'o(-, ■)) on il x D. 

The notation | • | will be used to signify one of the following. When applied to a real number 
X, \x\ will denote the absolute value of the number x; when applied to a vector u, \v\ will stand 
for the Euclidean norm of the vector w; and, when applied to a square matrix A, \A\ will signify 
the Frobenius norm, [tr(A"'" A)] ^ , of the matrix A, where, for a square matrix B, tr(i?) denotes the 
trace of B. 



3 Existence of a solution to a discrete-in-time problem 



Let 



H ■.= {weL^in):ya:-w = 0} and V :^ {w e H^ifl) :Vx ■ w = 0}, 



(3.1) 
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where the divergence operator V^; • is to be understood in the sense of vector- valued distributions 
on n. Let V' be the dual of V. Let S : V' V he such that S v is the unique solution to the 
Helmholtz-Stokes problem 

/ Sv-wdx+ / Yx (S v) : Yxwdx = {v,'w)y Vw^V, (3-2) 

where (•, •)y denotes the duality pairing between V' and V. We note that 

{u,Sv)y = \\Sv\\l,^^^ yveV'DiHlm', (3.3) 

and \\S ■ llifi(o) is a norm on V' . More generally, let Va denote the closure of the set of all 
divergence- free CJf (57) functions in the norm of HQ{il.)r\H'^{n), cr > 1, equipped with the Hilbert 
space norm, denoted by || • \\v„, inherited from H°'{il.), and let signify the dual space of Va-, 
with duality pairing (•, ■)v^- 

For later purposes, we recall the following well-known Gagliardo-Nirenberg inequality. Let 
r E [2,oo) if c? = 2, and r G [2,6] if c? = 3 and 9 = (5 ~ 7)- Then, there is a constant 
C = C(0,r,d), such that, for all 77 e H\n): 

||r/|U..(o) <C^||r/||i^(%||ry||?,,(^). (3.4) 

Let J" e C(R>o) be defined by J"(s) := s (logs - 1) -I- 1, s > 0. As lims_+o+ -^(s) = 1, the 
function J- can be considered to be defined and continuous on [0, 00), where it is a nonnegative, 
strictly convex function with J-(l) = 0. We then introduce the following assumptions on the data: 

dn e C°'^; uoeH; $0 ■= — > a.e. on fl x D with 

~ ~ M 

J"('0o) G X D) and / M{q) ijjo{x, q) dq = I for a.e. x & ft; 

J D ~ ~ ~ ~ 

and f L^{0,T-V'). (3.5) 

Here, i^(i7 x D), for p E [l,oo), denotes the Maxwellian- weighted space over il x D with 
norm 



^\\Ll,{nxD) ■■= < / Ml^l^dgd^ 



nxD 



Similarly, we introduce L\j{D), the Maxwellian- weighted space over D. 
On defining 

\MHlAi-l>,D) ■■= \ f M \\ip\^ + \Y,(p\^ + \Yg^\^] dqdx\\ (3.6) 
UnxD L J ~ J 

we then set 

X = Hl,{nxD) [^eLl^{nxD):\mHi^,,,D)<^}- (3-7) 

Similarly, we introduce HIj{D), the Maxwellian- weighted space over D. It is shown in 
[Appendix A| that 

C^iD) is dense in HIj{D) and hence C°°{n, C^{D)) is dense in X. (3.8) 

We have from Sobolev embedding that 

H\n-Ll,{D)) ^ L%Q;Ll{D)), (3.9) 

where s e [1, 00) if d = 2 or s e [1, 6] if d = 3. Similarly to p.4p we have, with r and 6 as defined 
there, that there exists a constant C, depending only on f2, r and d, such that 

Mwl^iD)) < Mm^n.,LUD)) y$eH\n;Ll,iD)). (3.10) 
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In addition, we note that the embeddings 

Hij{D) ^ LIj{D), (3.11a) 
HU^ xD) = L\n; Hl,{D)) n H\n; LIj{D)) ^ Llj{n xD) = L\n; LI,{D)) (3.11b) 

are compact; see [Appendix D| and [Appendix Y\ respectively. 

Let X' be the dual space of X with L\i{Q, x D) being the pivot space. Then, similarly to p.2p . 
lei Q : X' ^ X be such that Qrj\s the unique solution of 



M 



nxD 



dq dx 



= {M fj, if) 2 y^e X, 



(3.12) 



where (M ■,-)x denotes the duality pairing between X' and X. Then, similarly to (|3.3I) . we have 
that 



(3.13) 



and 11^ • II is a norm on X' . 

We recall the Aubin-Lions-Simon compactness theorem, see, e.g., Temam [39] and Simon [38] . 
Let Bq, B and Bi be Banach spaces, Bi, i = 0,1, reflexive, with a compact embedding Bq ^ B 
and a continuous embedding B ^ Bi. Then, for > 1, i = 0, 1, the embedding 



{?7 e L"o(0,r;Bo) : |f G (0, T; 6i) } L"''(0,T;B) 



(3.14) 



is compact. 

Throughout we will assume that (ll.4a| -c) and (|3.5p hold, so that (|1.7p and (|3.11al b) hold. We 
note for future reference that (|2.3al) and (11.71) yield that, for (p E L\j(fl x D), 



|C.(M(^)|2 dx = 



M(pUlq.,qJdq 



dx 



< M{U^r\q,\^dq 



D 



M\ip\^dq dx 



nxD 



< C 



[ M\ipfdqdx), 



i^\,...,K, 



(3.15) 



where C is a positive constant. 

We now establish a simple integration-by-parts formula. 

Lemma 3.1 Let (p G H\^{D) and suppose that B G K'*^'' is a square matrix such that tr(-B) — 0; 
then. 



/ M,{qi){Bqi)-Vq^(p{q)dq., 



AU{q^)^{q)U^{^\q,\^)q^qJ -.Bdq, 



K 



for a.e. {qi,...,q^-i,q^+i,...,qK) G X Dj \ , i = l,...,K, (3.16a) 

~ ~ ~ ~ \j = l,0^i I 

/ M(g) ^(Bg,)-V,,^(«)d<7= / M{q)^{q)Y,U[{\\1^?)q^qJ ■■ Bdq. (3.16b) 



Proof The set C^{D) is dense in H\^{D), see [Appendix A[ hence, there exists a sequence 
{<?n}n>o C C§°{D), converging to ip in H\i{D). The identity p.l6ap with (p replaced by (pn 
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is easily verified. First, on applying the classical divergence theorem for smooth functions, and 
noting (|1.6p and that tr(i?) = 0, we obtain that 



Br(0) 



M,(g,)(SgO • V,,^„(<7)dg, = / M,{q,) $„{q) q,qJU^{^\q,\^) : B dq. 



Br{0) 



M,^„{Bq,) -j^dSiq^) 
3r{0) ~ 



K 



for a.e. (gi,...,q,,_i,q,,+i,...,qi<-) G X Dj \ , i^l,...,K. (3.17) 



Here we have replaced Di by _Bi^(0), the ball of radius R centred at the origin. On noting (jl.4cp . 
we have that, as i? — > cxd. 



M,{q,) [pn{q) (Bq,) • ^ d5(qO 
j(0) ~ ^ ^ \qi\ 



\Vn\\L=°{D)- 



(3.18) 



Hence, on passing to the limit i? — )• cx) in p.l7p the boundary term vanishes and we obtain p.l6ap 
with if replaced by Then, (|3.16ap itself follows by letting n oo, recalling the definition of the 

norm in H\j{D) and (|1.7p . Finally, multiplying (I3.16ap by M/Mi, integrating over ^X^^ Dj^ 

and summing from i = 1 to K yields the desired result (j3.16bP . □ 

We now formulate our discrete-in-time approximation of problem (P^.l) for fixed parameters 
e G (0, 1] and L > 1. For any T > 0, let TV Ai = T and t„ = n At, n = 0, . . . , N. To prove 
existence of a solution under minimal smoothness requirements on the initial data, recall p.Sp . we 
introduce G V such that 



dx — / Uq ■ V dx 

Jn~ ~ ~ 



V-y e V: 



and so 



/ + Ai|y,u"p]dx < / \uo\^dx<C. 
Jn ~ Jn 



(3.19) 



(3.20) 



In addition, we have that u'^ converges to uq weakly in H in the limit of At — > 0+. For p G [1, oo), 
let 



{ip G X D) : ^ > a.e. on x D and 



/ M{q) (p{x, q) dq < 1 for a.e. x G fl}. 

Jd ~ ~ ~ ~ 



(3.21) 



Analogously to defining w*^ for a given initial velocity field Mqj shall assign a certain 
'smoothed' initial datum, t/j'^, to the initial datum i/jq- The definition of ■^^ is delicate; it will 
be given in Section HI All we need to know for now is that there exists a -0°, independent of the 
cut-off parameter L, such that: 

^P°eZ,; \ W ' / MTii,°)dqdx< MTi^o)dqdx. (3.22) 

[ iIjO e Hii{Q X D); JnxD ~ JnxD 

It follows from (jX^ and ([LTO]) that /3^(V^°) G Z2; in fact, /3^(^°) G L°°{n x D). 
Our discrete-in-time approximation of (Pe,L) is then defined as follows. 
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(P^*) Let ^ u° e V and 7/^^^ = I3^{^^°) e 2^2- Then, for n = l,...,iV, given 

te') e F X find (y^^i, Cx) e r X (X n Z2) 



At 



such that 
w dx + v / u" ^ : If dx 



70 = 



M 



= [ r-wdx-kV] [ C^{Miil\^) -.X^^wdx 



Vw e (3.23a) 



JnxD 



'— — - — (3 do dx 
At 1 ~ 



Vg;(^dq dx 



InxD 



Vx^dq dx = £ X; 



(3.23b) 



v. 



(3.24) 



where, for t G [t„„i,i„), and n = 1, . . . ,iV, 

/^*^+(-,i) = r(-):=^/*" /(•,i)dtG 

It follows from ([SS]) and ([3?24|) that 

/'^*'+ ^ / strongly m L^{Q, T; F') as A< ^ 0+. (3.25) 

We note here that since the test function w in p.23al) is chosen to be divergence-free, the term 
containing the density p in the definition of r (cf. (|2.2p ) is eliminated from p.23ap . and will play 
no role in the rest of the paper. 

In order to prove existence of a solution to we require the following convex regularization 



J. e.i^n 1 — 

of F defined, for any 5 G (0, 1) and L > 1, by 



Hence, 



^^ + s(log(5-l) + l fors<<5, 
Fl'{s):=l T{s) =s(logs-l) + l for 6 < s < L, 
+ s (log L - 1) + 1 for L < s. 

+ log (5—1 for s < 5, 
[^5]'{s)=\ logs for(5<s<L, 
f + log i - 1 for L < s, 

(5-1 for s < (5, 
[J'^fis) - { for ,5 < s < L, 

L-i for L<s. 



(3.26) 



We note that 



for s < 0, 



(3.27a) 
(3.27b) 

(3.28) 



- C'(L) for s > 

and that [J-^Ws) is bounded below by 1/L for all s e K. Finally, we set 

Phs) ■■= m]"r\s) = max{/?^(s),^}, 

and observe that I3g{s) is bounded above by L for all s G R. Note also that both and /3g ari 
Lipschitz continuous, and that their respective Lipschitz constants are equal to 1. 



(3.29) 
e 
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3.1 Existence of a solution to {Pf'D 

It is convenient to rewrite p.23ap as 

where, for all lUi £ i/g(il), i — 1,2, 
and, for all w e ff^fl) and € Lli{fl x D), 



(3.30) 



W2dx + Ati' / iL'i : Va; z«2 dx, (3.31a) 



We note that 



K 



^ - w-Atk ^C,(M^) : V^- 



dx. 



(3.31b) 



(W • V,j;)wi 



W2 dx 



iv-Vx)w2 ■ widx Wv<eV, ywi,W2 <E H^{fl), 



(3.32) 



and hence &(•, •) is a continuous nonsymmetric coercive bilinear functional on i?Q(f2) x Hq{Q). In 
addition, on recalling p.lSp . ib{'f){-) is a continuous linear functional on V for any if S L^(r2x Z?). 
For r > d, let 



:= e L'^(17) : y" « • V^wda; = e VK1^^(17)| 
It is also convenient to rewrite p.23b[) as 

where, for all ipi, (p2 G X, 



(3.33) 



(3.34) 



0(^1,^2) 



M \(pi^2 + At 



At 



2A 

and, for all v G ^ L°^{n x D) and ^ £ X, 



XI XI "^iM ■ '^9.'?2 dq dx, 



(3.35a) 



4(i',?7)(^) := 



M 



nxD 



K 



dq dx. 



It follows from ((333| and dSJ]) that, for r > d. 



(3.35b) 



(3.36) 



and hence a(-, •) is a continuous nonsymmetric coercive bilinear functional on X x X. In addition, 
we have that, for ah v e H^{^), rj e L°°(fi x £)) and tf e X, 

|4(w,5?)(^)| < II^'?,l^IIl2^(oxd) \MLl,(nxD) 

+ ^t[ M\q\'^dq] \\v\\L-^{nxD)\\Yxv\\L''{n)\\'^q'f\\Ll,inxD)- (3-37) 
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Hence, on noting that V'^l^ G -^2, (|l-4cp and (jl.Sp . £a{v,rf){-) is a continuous linear functional on 

X for all V e and ^7 G L°°(f2 x D). 

In order to prove existence of a solution to p.23a[ b). i.e. p.30p and (|3.34[) . we consider a 
regularized system for a given S € (0, 1): 

Find (u^^^, ■0^^ 5) e F X ^ such that 

K<L,5,Jf) = 4(^;U,5)H VweV, (3.38a) 

^) = V^e X. (3.38b) 

In order to prove existence of a solution to p.38al b). we consider a fixed-point argument. 
Given e Llj{n x D), let (u*, -0*) e F x X be such that 

b{u*, w) = 4(V')(w) e V^, (3.39a) 

a(V^*, 0) = iaif, PfWm e X. (3.39b) 

The Lax-Milgram theorem yields the existence of a unique solution to p.39al b). and so the overall 
procedure ()3.39a| b) is well defined. 

Lemma 3.2 Let G : L\j{VL x D) ^ X C i^/(f^ x D) denote the nonlinear map that takes ^jj to 
Tp* = G{'ip) via the procedure (|3.39al b'). Then G has a fixed point. Hence there exists a solution 

(Cl.^'Cl.^) eFxX to eSHib). 

Proof Clearly, a fixed point of G yields a solution of p.38a[ b). In order to show that G has 
a fixed point, we apply Schauder's fixed-point theorem; that is, we need to show that: (i) G : 
L\[{{1 X D) L,\,i{n X D) is continuous; (ii) G is compact; and (iii) there exists a C* G IR.>o such 
that 

\mLl,(n.D) < C. (3.40) 

for every 'tp E Lj^ti^ x D) and k E (0, 1] satisfying -ip = KG{ipi). 
Let {V'*^^-'}p>o be such that 

ip^P^ ^ ip strongly in i|f (17 X D) asp-J-oo. (3.41) 

It follows immediately from p.29p that, for any r E [2, 00), 

(^(P)) ^ ^f{ip) strongly in L''{n x D) as p ^ 00, (3.42a) 

and from (I3.15p . for i = 1, . . . ,K, 

C,{Mi)^P'^) ^ G,{Mi!) strongly in asp ^00. (3.42b) 

We need to show that 

^^(P) := G(^(P)) ^ G(V^) strongly in (17 X L>) as p ^ 00, (3.43) 
in order to prove (i) above. We have from the definition of G, see p.39al b). that, for all p > 0, 

aljfP\^)=Uvy\pU4'^P^)W) yipEX, (3.44a) 

where w'-^^ E V satisfies 

b{v'^P\w) = eb{i^^P^){w) Vw E V. (3.44b) 



15 



Choosing (p — j^'p^ in p.44ap yields, on noting the simple identity 

1{sx- S2)sy= s\ + [sx- s^f - s\ Vsi,S2e]R, (3.45) 

and that e V, (^3]) . the bound on the second term in (|3.35bp . as in ([333), ([TTic)) . (|L5)) and 
(P:^ tha"t,'for allp > 0, 



M 

nxD 



dq dx 



< I M IV'^'lT d? dx + C(i) At ( \\/^v'^P^\'^dx. (3.46) 
JnxD ' ~ ~ Jn ~ ~ 

Choosing w = v'-p'^ in (|3.44bp . and noting (|05l) . (|332|) . (|3?T5l) . (|3^ . a Poincare inequality and 
(|3:4T|) yields, for ah i > 0, that 

[ \\v'-P^\^ + \v'-P^ -u^~l\^] dx + Ati^ [ |V^f(f)pda; 

< I WrL\^ ^x + C M\\S rfm{n)+C At ( M \i^^P^ \^ dq dx < C . (3.47) 

Combining (|3.46p and p.47p , we have for all p > that 

\\v^P^h + \\vJP^\\HHn)<C{L,{At)''). (3.48) 

It follows from p.48p . p.9p and the compactness of the embedding (|3.11bp that there exists a 
subsequence {irf-^''\v''^''^)}pk>o and functions f} £ X and v£V such that, as pk — >■ oo, 

fj'-P'^^ fj weakly in L'iQ; Lli{D)), (3.49a) 

M^V^rv'-^-^^M^V^^ weakly in {n X D), (3.49b) 

M3 Vq^P"^ -> Vqf] weakly in L^{n x D), (3.49c) 
JyfP'-) ^ ?y strongly in i|f(17 x D), (3.49d) 

^ w weakly in V; (3.49e) 

where s e [1, oo) if d = 2 or s G [1, 6] if d = 3. We deduce from (|3.44bp . p.31al b). ((3:49e| and 
(I3.42b|) that v e V and G X satisfy 

b{v,w) =4(V^)(w) VtwGF. (3.50) 

It follows from p.44ap . p.35a[ b). p.49a| -e). and p.42ap that ff, ijj € X and i; G Jf, satisfy 

a(57,^)=4(t;,/3^(V^))(^) V^GC°°(17;Co°°(i^)). (3.51) 

Then, noting ([3371) . (fr4c|) and ([3^ yields that (|33T|) holds for all <^ G X. Combining this X 
version of (|3.5ip and p.SOp . we have that rj — G{tp) G X. Therefore the whole sequence 

ffp'> EE G(V^(P)) G{$) strongly in Llj{n x D), 

as p oo, and so (i) holds. 

As the embedding X ^ Lj^ii^ x D) is compact, it follows that (ii) holds. It therefore remains 
to show that (iii) holds. 
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As regards (iii), = kG{^) implies that {v,^j} £VxX satisfies 

b{v,w) = £b{4')iw) yw G V, 

Choosing w = v in p.52ap yields, similarly to p.47p . that 

dx + Ati' / \Yxv\'^dx 



(3.52a) 
(3.52b) 



ie,L I lHeX 



At 



(3.53) 



Choosing if — [J^f']'(V') in p.52bp . and noting the convexity J-g , p.29p and that v is divergence- 
free, yield 



At ^ ^ 

i=l 3 = 1 



nxD 

K 



Af V,.0-V,,([J-f]'(7A))dg dx 



< 



kA^ / Af (t(u) • Vg.V' dq dx 
j^-^ JnxD » ~ ~ ~ ~ ~ 

KAt V / C^iMiJj) : a{v)dx, 



(3.54) 



where in the transition to the final equality we applied ()3.16bp with B :— g;{v) (on account of 
it being independent of the variable q), together with the fact that tv{a{v)) = V^: ■ v — 0, and 

recalled p.3ap . Next, on noting p.29p and that y"^^ G F: i* follows that 



M<2V-V,([J-i^]'(^))d<zdx= / Mu 



nxD 



^1 ^ 



Vxip da; 



nxD ~ ■ iS^itp) ~ 

Af<^i • V,([Gf]'(^))dgdx = 0, (3.55) 



where Gf G C^'^iR) is defined by 



J_„2 , (S-L) 
2(5* ^ 2 



G^(.):=<; ,s-| 

2L^ 



if s < (5, 
if s e [(5, L] 
if s > L; 



(3.56) 



and so [G^]'is) = s//3|'(s). Combining (l333| and (IXKill . and noting (13351) . jS^ZEl, d!!]), (jM 
and a Poincare inequality yields that 



\v\ + \v — u 



e.L 



kL-^ At 



M 



nxD 



dx + KAtiy / \Vxv\ dx + k / MTs{ip)dqdx 
Jn " ~ ~ JnxD ~ ~ 

£|V.^|2 + |||V,V^P] dqdx 



< KAt{f\v)v + - / K^^rdx + Zc / Af J-5^(KV;';f:^)dgda; 



<'^Atiy [ \\/xv\^dx + KAtC{i^'^)\\S f 
2 « ~ ~ 



nil 2 



lu^'X^pdx + fc / A/J-f (KVj;'i^)dgda;. 



(3.57) 



17 



It is easy to show that J-^{s) is nonnegative for all s e R, with J-g{l) ~ 0. Furthermore, for 



any k e (0, 1], 



Thus we deduce that 



if s < or 1 < K s, 
if < KS < 1. 



(3.58) 



Hence, the bounds l|3.57p and p.58p . on noting p. 281) . give rise to the desired bound (|3.40p with 
C* dependent only on 6, L, At, k, v, /, y"^^ and Tp^^ ■ Therefore (iii) holds, and so G has a 
fixed point. Thus we have proved existence of a solution to p.38al b'). 



□ 



Choosing w = u"^^ in (|3.38aP and ^ = [^s]' ('^e l s) (|3.38bp . and combining, then yields, 
similarly to (|3.57p . that 



,n-l |2 



dx + k 



n-xD 



+ 



At 

k L^^ gq 
2A 



^ / \V.ulL,s\'dx + kL-'s [ M\vJl\s\'dqdx 
^ Jsi ~ ~ ~ JnxD ~ ~ ~ 



M\Vgi^lrs\'dqdx 



<AiCKi)ii5rii?,,(^) + i / \u:j}\'dx + k f MF^{r,:j})dqdx 

< C{L), (3.59) 

where, on recalling that V""!^ ^ Oj C'(L) is a positive constant, independent of 5 and At. We are 
now in a position to prove the following convergence result. 

Lemma 3.3 There exists a subsequence (not indicated) o/ {(w"^ i-, V'el 5)}(5>0; and functions 



g V and ip"^ e X n Z2, n e {1, . . . , N}, such that, as S 0+, 



!ie,L 

UeLS^HeL wcakly in V, 

'ii^,L,S-^ Ue.L strongly in L''{VL), 

where r € [1, 00) if d = 2 and r G [1, 6) if d = 3; and 

M5 ^-^^^ g M3 ipi^ weakly in L^{VL x D), 

M5 Vg?^"i^5 ^ Af5 VgV'"L weakly in L^{VL x L>), 
Vx'$e,L.s ^ Va;V'"L wcakly in L^{Q. x D), 
^l^^, i^l^ 



where s G [2, 00); and, for i — \, . . . , K , 



strongly in L^{Q x D), 
strongly in L'^{Q x D), 

strongly in L'^{Q). 



(3.60a) 
(3.60b) 

(3.61a) 
(3.61b) 

(3.61c) 

(3.61d) 
(3.61e) 

(3.61f) 



In addition, (u^ 1^,1^" j^) solves hS. 23a\ h] for n = 1,...,A^; consequently there exists a solution 
{(y"L: €L)}n=i to (Pfi), wtth yp^^ e V and ^ G Z2 for all n ^ 1, . . . , N . 
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Proof The weak convergence results p.60a|) and p.61a[) and the fact that "0"^ > a.e. on J7 x D 
fohow immediately from the first two bounds on the left-hand side of (|3.59l) , on noting p. 281) . The 
strong convergence result (I3.60b[) for u"^ g follows directly from (|3.60ap . on noting that V C Hq{^) 
is compactly embedded in L^{^) for the stated values of r. 

It follows immediately from the bound on the fifth term on the left-hand side of (|3.59p that 
(j3.61bp holds for some limit g G L^{n x D), which we need to identify. However, for any ?/ e 
C^ifl X D), it follows from (|1.6p and the compact support of on D that 

[Vq ■ {M^ v)]/M^ e L^{n X D), 
and hence the above convergence implies, noting (|3.61ap . that 



/ g -ridqdx ^ - i^^^ s ' , , i ~ M 

JnxD- ~ ~ ~ JU-kD M 2 ^ ~ 



/ ^ „ , ~ dqdx = - re,L V, • (M5 ^) dq dx 

JnxD M2 ~ ~ JnxD ~ ^ ^ ~ 



(3.62) 

as (5 —7- 0+. Equivalently, on dividing and multiplying by under the integral sign on the 
left-hand side, we have that 

/ M^^g ■ JVn-qdqdx^ - ' (^^^?7) dg dx yj]eCl(nxD). 

JnxD ~ ~ ~ ~ JnxD ' ~ ~ ~ ~ ^~ 

Observe that i] £ Cl{il. x D) i-> M^rj e x Z?) is a bijection of Cl{Q x D) onto itself; thus, 

the equality above is equivalent to 

/ M-^g-xdqdx = - f ^^^^ q ■ x) dq dx Vx e x i?). 

JnxD ~ ~ ~ JnxD 

Since C^iflxD) dCHflxD), the last identity also holds for all rj € C^iflxD). As Af5 g L°°{D) 
andM-5 g Lj^^p), it follows that M-^c,eL2^^(r]xD) and & L'{^^{i^x D). By identification 
of a locally integrable function with a distribution we deduce that M^^g is the distributional 
gradient of -0"^ w.r.t. q: 

M^^g ^Yq^e.L inV'irixD). 

As M^ig e Lfoci^ X -C)), whereby also Ygi'e.L ^ I'lod^ ^ follows that 

g = MiY,i'lL^Ll,{nxD). 

However, the left-hand side belongs to L^{fl x D), which then implies that the right-hand side 
also belongs to L'^{^ x D). Thus we have shown that 

g = M^Vq^lL e L^i^ X D), (3.63) 

and hence the desired result p.61bp as required. 

A similar argument proves p.61cp on noting p.61ap . and the fourth bound in p.59p . 

The strong convergence result (I3.61dp for -0"^ g follows immediately from (I3.61a| -c) and (I3.11bp . 
Finally, the desired results p.GlcI f) follow immediately from (|3.61dp . (|3?29| . ([23a|) and (|3?T5|) . 

It follows from (|3J0al b). (|33Tbl -f). (IS^Tal b). (j335al b). (|337l) and that we may pass to 
the limit, 5 0+, in (I3.38a[ b') to obtain that (u^^, 0^^) £V x X with 4)"^^ > a.e. on ^2 x D 
solve ([3:301 and ([3341) . i.e. ((SfSSal b). 
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Next we prove the integral constraint on ip" i^- First we introduce, for to = n — 1, n, 

C7:l{x):= f M{q)i;^^ix,q)dq. (3.64) 
Jd ~ - ~ ~ 

As "0"^ G X and '(/'"^^ G 2^2: we deduce from the Cauchy-Schwarz inequahty and Fubini's theorem 
that C"l ^ H^{fl) and C^l^ ^ i^(il). We introduce also the following closed linear subspace of 
X = Hlj{n X D): 



H\n)<»l{D) -.^ i^lfe Hliinx D) : (f{-,q*) = (f{-,q**) for ah g*, g** G 



(3.65) 



Then, on choosing (p = £ H^{il,) (g) 1(1?) in p.23b[) . we deduce from (I3.64p and Fubini's theorem 
that 



Q.L-C 



At 



■\'^ipdx = yipeH^n). (3.66) 



By introducing the function z™^ := 1 — ^™^, m = n—1, n, we deduce from (|3.66p . and as u"^^ 
is divergence-free on with zero trace on dfl, that 



At 

Let us now define by 



if dx 



■V^ipdx = yipeH'^ifl). (3.67) 



[x]± :=i(x±|x|) 



the positive and negative parts, [x]j^- and [x]^, of a real number x, respectively. As 4'"~l ^ ■^s, 
we then have that [^"2^1- ~ ^ ^-^^ ^- Taking = [^^ lI- ^ ^^s^ function in p.67p . noting 
that this is a legitimate choice since [z^ G H^{Vl), decomposing z™^, m = n — l,n, into their 
positive and negative parts, and noting that w"^^ is divergence-free on V, and has zero trace on 
dft, we deduce that 

\\[zl^U\' + Ate\\Y..KL]-r=0, 
where || • || denotes the L'^{^l) norm. Hence, [2:"^]- = a.e. on ft. In other words, 2;"^ > a.e. 
on ri, which then gives that Qj^ < 1 a.e. on fi, i.e. ■(/'"^ G Z2 as required. 

As (u^^,'0°^) G V X Z2, performing the above existence proof at each time level tn, n = 
1, . . . , iV, yields a solution {(u?,^, V^?,i)}^^i to {P^X). □ 



Having shown in Lemma [331 that problem (P^^) has a solution, we shall next develop suitable 
estimates on this solution, independent of the cut-off parameter L. 



4 Entropy estimates 

Our starting point for the analysis here is the final result of the previous section, stated in Lemma 
13. 3[ concerning the existence of a solution to the discrete- in-time problem (Pf^^). The model (P^^) 
includes 'microscopic cut-off' in the drag term of the Fokker-Planck equation, where L > 1 is a 
(fixed, but otherwise arbitrary,) cut-off parameter. Our ultimate objective is to pass to the limits 
L 00 and At — !• 0+ in the model (P^^), with L and At linked by the condition At = o{L^^), as 
L ~^ 00. To that end, we need to develop various bounds on sequences of weak solutions of (P^^) 
that are uniform in the cut-off parameter L and thus permit the extraction of weakly convergent 
subsequences, as L — >■ 00, through the use of a weak-compactness argument. The derivation of 
such bounds, based on the use of the relative entropy associated with the Maxwellian M, is our 
main task in this section. 
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Let us introduce the following definitions, in line with f|3.24p 



At 



■y",L(-) + 



At 



ie,L 



{■), te[tn^i,tn], n^l,...,N, 



and 



At,+ / 



(•,t) n^l,...,N. 



(4.1a) 



(4.1b) 



We shall adopt y^^^'^"* as a collective symbol for u^^, u^^'*- The corresponding notations V'^l 
V'^^, and ■0^1'^ are defined analogously; recall p.20p and p.22p . 



We note for future reference that 

At „.At,± 



i e (t„_i,t„), n = 1, 



.,7V, 



(4.2) 



where t+ := t„ and t„ := t„_i, with an analogous relationship in the case of V'i^i- 

Using the above notation, p.23a[ b') summed for n = 1, . . . , can be restated in a form that 
is reminiscent of a weak formulation of (|l.lal -d): Find {uf*£^{t),'ipf'*£'^{t)) G V x {X Ci Z2) such 



that 



Jn 



'E,L 



dt 



Jo. 

T 



■ w dxAt 



I At.- v7 \ At.- 



At. J 
X "■e,L 



Ax dt 



if 



At,+ 



K 



dt 



\/w e L\0,T;V), 



(4.3a) 



M 



JilxD 
T 



dt 



(pdq dx dt 



TAt,+ At,- 7At 



Vx^P dq dx dt 



I I M {eVx^^i^ 
Jo JnxD L ~ 

I fT r K K 

^'^ E E ^'^^ Y-^. V^e'^r • V,, ^ dq dx dt 

/ / M 5][a(7.fl'+) q,] /3^(^fl'+) • V,,^ dg dxdt = 
Jo JnxD -1 » ~ ' ~ ' ~ ~ ~ 



V^eLi(0,T;X); 



(4.3b) 



subject to the initial condition uf^i(-,0) = m" G F and ijj^l{-,-,0) = /3^(-0''(-, •)) G Z2. We 
emphasize that (|4.3a[ b) is an equivalent restatement of problem (P^^), for which existence of a 
solution has been established (cf. Lemma [ 



Similarly, with analogous notation for {C"l}^=oi P-661) summed for n = 1,...,N can be 
restated as follows: Given uf*£'{t) G V solving B3al b). find Ct,L^{t) G /C := {t? G H\n) : t] G 
[0, 1] a.e. on ft} such that 



Jn 



At 

'£,L 



dt 



ipdxdt 



Jn 



At,+ At,- /-At,+ 



U 



6,L '^E,L 



Vxy^dx dt — 

\f(p e L\0,T;H\n)), (4.4) 
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subject to the initial condition C^i(-,0) = f^M{q)(3^{ip°{-,q))dq; cf. ((XM)) and recall that 

Once again, on recalling (|3.64p and (I3.66p . we have established the existence of a solution to 
(gJl) and that 

Cl{x,t) = / M(g)V^fi(x,g,i)d(Z for a.e. {x,t) e O x (0,T). (4.5) 

In conjunction with , defined by (|l.I3p . we consider the following cut-off version J^^ of the 
entropy function J" : s e M>o J^{s) = s(log s — I) + I e R>o: 

f°fr'^^'' . '-'-'^ (4-6) 

I Sr" + s(logL-l) + l, L<s. 



Note that 



and 



Hence, 



, f log., 0<,s<L, 

^ ' \ f +logL-l, L<s, ^ ^ 



r (4.8) 

1 L' 



/3^(s) =min(s,i) = [(J-^)"(s)]-\ s S R>o, (4.9) 
with the convention l/oo :~ when s — 0, and 

> J-"(s) = s-i, s e R>o. (4.10) 

We shall also require the following inequality, relating to J^: 

T^{s)>T{s), seR>o. (4.11) 

For s > 1, this follows from (j4.10l) . with s replaced by a dummy variable a, after integrating twice 
over a G [l,s], and noting that (J^^)'(l) = J"'(l) and (J"^)(l) = For s e [0,1], we have 

J^^(s) = J'is) of course, by definition. 

4.1 L-independent bounds on the spatial derivatives 

We are now ready to embark on the derivation of the required bounds, uniform in the cut-off 
parameter L, on norms of u^^'^, V-'f^i'^ ^-nd C^j^^ ■ As far as yf^^'^ is concerned, this is a relatively 
straightforward exercise. We select w = X[o,t] uti,'^ ^ test function in (|4.3al) . with t chosen as 
tn, n € {1,...,A^}, and X[o,t] denoting the characteristic function of the interval [0,t]. We then 
deduce, with t ^ tn, that 



At Jo 

1 



<IImo||^ + - / ||/"''+(s)|lt.,ds 

K 

-2k ' ' 



f I M{q) f]g.<7?(7;(i|g.H V7fi'+ : V . ufl+ dq dx ds , (4.12) 
Jo JnxD ~ j^^i ~ ~ ~ 

where, again, \\ ■ \\ denotes the norm over fl and we have noted p.20p . We intentionally did 
not bound the final term on the right-hand side of (|4.12p . As we shall see in what follows, this 
simple trick will prove helpful: our bounds on ip^L^ below will furnish an identical term with 
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the opposite sign, so then by combining the bounds on u^^"*" and "0^1' this pair of, otherwise 
dangerous, terms wih be removed. This fortuitous cancellation reflects the balance of total energy 
in the system. 

Having dealt with uf''^£'^ , we now embark on the less straightforward task of deriving bounds 

on norms of ?Af^|^'^ that are uniform in the cut-off parameter L. The appropriate choice of test 

function in (j4.3bp for this purpose is ^ = X[o,t] i.^^)' {.i^f]^^) with t = tn, n G {1, . . . , N}; this can 
be seen by noting that with such a ip, at least formally, the final term on the left-hand side of 
(|4.3b[) can be manipulated to become identical to the final term in (I4.12p . but with opposite sign; 
and this will then result in the crucial cancellation of terms mentioned in the previous paragraph. 
While Lemma [3.31 guarantees that ip'^*£'^ {■ , ■ , t) belongs to Z2 for all t G [0,r], and is therefore 

nonnegative a.e. on 51 x 13 x [0,T], there is unfortunately no reason why i/'i^i'^ should be strictly 

positive onQx D X [0, T], and therefore the expression (J^^)' {ip^*£'^) may in general be undefined; 

the same is true of which also appears in the algebraic manipulations. We shall 

circumvent this problem by wor king with ( J"i')'(^^[^+ + a) instead of ( J"^)'(V^^['+), where a > 0; 

since 1/'^^'''' is known to be nonnegative from Lemma [331 {J^^)' {ipf^^ + a) and {J^^)" {^fl^'^ + a) 

are well-defined. After deriving the relevant bounds, which will involve J^^{'il!^*£^ + a) only, we 

shall pass to the limit a — > 0+, noting that, unlike {J^^Y {tp^'£^) and (J^^)" {'ip'^'^^), the function 

is well-defined for any nonnegative i'f^'^- 
Thus, we now take any a G (0, 1), whereby < a < 1 < L, and we choose 

^=X[o,t](-^^)'(V^i^L+ + a), witht = t„, ne {l,...,iV}, 

as test function in (j4.3bp . As the calculations are quite involved, we shall, for the sake of clarity of 
exposition, manipulate the terms in (I4.3bp one at a time and will then merge the resulting bounds 
on the individual terms with (|4.3ap to obtain a single energy inequality for the pair {u'^*£'^ , V'^l^)- 
We start by considering the first term in (j4.3bp . Clearly + a) is twice continuously 

differentiable on the interval (—a, 00) for any a > 0. Thus, by Taylor series expansion of s G 
[0, 00) I— > J-^{s + a) ^ [0, 00) with remainder, and c G [0, 00), 

(s - c) + a)= ^^(s + a) - ^^(c + a) + i(s - c)^ {J^^)"{es + (1 - 6)0 + a), 

with e G (0, 1). Hence, on noting that t G [0,r] 1-^ ^^X^'' ■'^) ^ ^ 

is piecewise linear relative to 

the partition {Q ~ to,ti, . . . ,t]\; — T} of the interval [0, T], 

"At 



Ti := / / Af — ^X[o,t](-F^)'(V^fr+«)dgdsds 
Ja JnxD OS 

M^ iijjfX + a) + a) dq Ax As 



Jn/.D 9s 

MF^ {4,fl+ {t) + a) Aq Ax - [ MF^ {fi^ + a) Aq Ax 
nxD ' ~ JnxD 



1 



+ 9a7 / / MiT^ne^p^r + (1 - ^)v^s"r + «) (cr - d^ as. 

Jo JnxD 

Noting from (|4.8I) that {F^)"{s + a) > l/L for all s G [0,oo) and all a > 0, this then implies, 
with t = t,i, n G {1, . . . , N}, that 

Ti > / MF^ (^/^f [■+ (t) +a)AqAx- [ MF^ {/3^ (V^" )+a)Aq Ax 
nxD ' ~ JnxD 

Mii^f'£+ -ijfl^-fAqAxAs. (4.13) 



1 



2AtL Jq 



nxD 
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The denominator in the prefactor of the last integral motivates us to link Ai to L so that 

AiL = o(l) asAt^O 

(or, equivalently, A< = o{L^^) as i — > cx)), in order to drive the integral multiplied by the prefactor 
to in the limit of i — > oo, once the product of the two has been bounded above by a constant, 
independent of L. 

Next we consider the second term in (j4.3bp . using repeatedly that V^^ • = and that 

uf*£'^ has zero trace on dfl: 



M 



JnxD 
t 



JnxD 



lo JnxD 

where in the last line we added in the form of 



a 

lo JflxD 

Hence, similarly to p.55p . 



A^yfi'" • Y^J^'^Yi^Ti'^ + a) dg dxds = 0. 



L2 := e 



M(J-^)"(V^^i'+ + a)|y,(^,^i^+ + dgdxds 

At.-/TAt. 



JnxD 



JnxD 



A/(J-^)"(Vfi'+ + a)|y.^,l^+P dg dxds 



J f2 X r> 



JnxD 



l/L 



if ?Af + a > L 



dq Ax ds 



M(j-^)"(V'rr + «)iy.crr m dxd^ 

MufX ■ y-[G^(^'i'L^ + ")] d£ dx ds, 



/ ^ O X D 

where denotes the (locally Lipschitz continuous) function defined on R by 



G\s) 



L 

ifs>i. 



On noting that the integral involving vanishes, (|4.10p then yields the lower bound 

T2>e/ / M(^f['+ + a)-i|y,(^;f['+ + a)|2d9dxds. 
Jo JnxD 



(4.14) 



Next, we consider the third term in (|4.3bl) . Thanks to (j2.5p we have, again with t = i„ and 
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ne {1,...,N}: 

K K 



= Vlf I i^'y'i^t}:^ + ") E E y^.^^'r • y,.^,^r a. a. 

^■^ Jo JflxD .^^ 

f I M(-^^)"(V^fr + «) Eiy,.v^fl^+p agaxa^ 



> — 
- 2A 



^ f f M{T^)"{$fl:+ + a)\Y,i^tL^\' dqdxds. (4.15) 
JnxD 



2A 



We emphasize here that, unhkc (|4.14p above, in (|4.15p we refrainea from using (|4.10p to further 
bouna {J-^y {■iljf'*£^ +a) from below by J^"(?/'^^'^+a) = (V'^^'^+a)"^. Performing this aaaitional 
lower bouna will be postponea until later, after a term similar to T3 that arises as a byproauct of 
manipulating term T4 below has been absorbea in term T3. 

We now consiaer the final term in (j4.3bp . with t = tn, n E {1, . . . , N}: 

T4:=- r [ M V[a(uf['+)g,]/3^(v^fi'+) • v,.x[o.t] + a^a^as 

Jo JnxD -T ~ ~ ~ ' ~ ' ~ ~ 



i=l 
K 



A/E[(Y-!f^L^)9.]/3^(V^S''') • (-^i)"(V?I'+ + «)v,.v^,^r aqa^a. 
=1 ~ ~ ~ ~ ~ ~ 

JaJnxD j=i ~ ~ ~ P (V'e,L +a) ~ ' ~ ~ 
M Et (Y- !f^L^) 9^ ] ■ Y^.^^'l^ dx ds 



+ 11 M y[(V,Mf*/+)q, 
fj M j2<l^qJU[{\W)^^f'^ ■.V^uf'y' dqdxds 

JoJnxD ._i ~ ~ ~ a ~ ^ ~ 



^9i^£^L^ aq dx ds 



Jo JnxD ~ ~ ' 



"^q.^tL^ dqdxds, (4.16) 



where in the transition to the final equality we appliea (I3.16b[) with B := Yx U^l^ {on account 
of it being inaepenacnt of the variable q), together with the fact that 

Summing (|i?T5|) ana (|i?TBl) yieUs, with t = <„ ana n G {1, . . . , iV}, the following 

inequality: 



MTH^^t'L^it) + a)dqdx+ I I M(^fi'+-v^fr)'dgaxas 

flxD ~ ~ ^^I-t^Jo JnxD ~ ~ 

M dq dx ds 

lo JnxD V'e.L ~^ ^ ~ ~ 

"° / * ^ M (J-^)"(^f + a) |V,^f I'+p dq dx ds 



2A 



JnxD 
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< 



MJ-^(/3^(V'°) +a)dg da: 



JO JnxD j^]^ ~ ~ ~ » ~ ^ ~ 

JO JnxD ,_i ~ ~ 



Vq^Vfi'"^ dgdaids. (4.17) 



Comparing (I4.17P with (|4.12l) we see that after muhiplying (|4.17l) by 2k and adding the resuhing 
inequahty to (|4.12p the final term in (I4.12p is cancehed by 2k times the second term on the 
right-hand side of (|4.17p . Hence, for any t — tn, with n E {!,..., N}, we deduce that 



At^,+ _^At,-„2, 



2k 



MJ-^(V^f['+(i) + a) dq dx + 



~ 
k 



nxD 



2ks 



AtL Jo 



M(^;^i'+-7Afi'-)2d(7dxds 



nxD 



M ~ dq dx ds 



'o JnxD ip^ £ + a 



+ 



aok /■* 

T 



JnxD 



<\\uor + - \\f^'^+{s)rv,ds + 2k MJ^^{P^{^°) + a)dqdx 



nxD 



K 



2k / Af^[(V: 



JnxD 



At, + \ 



Vg.^^['+ dqdxds. (4.18) 



It remains to bound the fourth term on the right-hand side of (I4.18p . Noting that is Lipschitz 
continuous, with Lipschitz constant equal to 1, and /3^(s -t- a) > a for s > (recall that < a < 
1 < iy), we have that 



< 1 - 



+ ^(J-L)"(V,At,+ +a) 



/3^(^^^i-+ + a) 



< 



^/a when V'^^'^ < 
when ^ L. 



< 



(4.19) 



With this bound and (ll.4cl) . for t = t„, n G {1, . . . , N}, wc then have that 



K 



2k / Mj2[{V,uf^)q, 



10 JnxD 

rt 



< 2k 



M\\7,uf^l+\\q\ 



10 JnxD 



^qiip'^^£~^ dq dxds 



|VqVe L^l dg dxds 



Af |V.<n |g| J(J-i)"«r +a) |V,<ri d<zdxd, 



/o JnxD 



At,+ 



2k^ / |V, 



Jn 



M\q\JiT^r{^f;£+ + a) |V,^,"['+| dg) dxds 



At, 
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IV. uf'^+\ / M(J-^)"(^^^*'+ + a) |V,V^^l^+p dq ] dxds 



< 2k^/^ ( [ ||V, ufl+W^ds) Y / / M(^^)"(^fl^+ + a) |V,^fl,'+|2 dg da:ds 



I I M(^-)"«r +a) |V,<n^ dgd.d. 



2XkCM f r UT-, _ At,+ ||2 



+ a^-^ U ||V,u,Xlrdsj, (4.20) 
where 

Ca/ / M\q\^dq. (4.21) 

Substitution of into (|i?T5)) and using (|i?TU)) to further bound (J^^)" {■ipf'^^ + a) from below 

by ■F"{^^t'£'^ + a) = (V^^['+ + a)-i and BUj) to bound J^^(^,^l,'+ + a) by :^(^,'^i+ + a) from 
below finally yield, for all f = t„, n e {1, . . . , N}, that 



i-^r wip + ^ - yf:Mp + ny. (.)f d. 



- 2fc / J-(^f (t) +a)dqdx + -^ [ [ M(^f - V^f )2 d<7 dg; ds 



<llyo||' + -/ \\f^''+{s)\\l,ds + 2k I AfJ-^(/3^(V^")+a)d<?ds 



2Afc Cm 
a 



riiy.y^r(s)|pd.. (4.22) 
flo Jo 

The only restriction we have imposed on a so far is that it belongs to the open interval (0, 1); 
let us now restrict the range of a further by demanding that, in fact, 

«<"<-^K^'2TI^)- ^'-''^ 

Then, the last term on the right-hand side of (|4.22p can be absorbed into the third term on the 
left-hand side, giving, for i = t„ and n € {1, . . . , N}, 

II At,-|-/,\||2 I 1 /" II At.+ At.-||2 1 

IKX (Oil + "-^■■f' ""t^'L " 

2AfcCA/\ /■* II At,+ , 



v — a 



ao 



/ \\v,uf'i+{s)rds 

Jo ~ ~ 



k 



2k I MT{ijf^{^{t) + a)dqdx+-— / Mi^f}:+ - V^^l'-)^ dq dx ds 



nxD ~ ~ '^tLjQ jQxD 



2k e I I M ~ dqdxds-l- 177 / / -At + dgdxds 



< hoir + - / ||/''*^+(s)||(^,ds + 2A: / Af J-^(/3^(0'J) + a) dg dx. (4.24) 

JnxD ~ ~ 
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Let us now focus our attention on the final integral on the right-hand side of (I4.24[) . which we 
label T^la) and express as follows: 

T5(a) := / MJ'^{P^{ip°) + a)dqdx 
JnxD 

= [ MJ'^{p^{$°) + a)dqdx+ [ MJ'^{p^{^°) + a)dqdx, 

where 

'&L,c. {{x,q)enxD : L-a<(3'^{$°ix,q))<L}. 

We begin by noting that 

/ MT^{(3^{ii°)+a)dqdx^ [ MT{(3^{i;^) + a)dqdx. 

For the integral over ^L,a we have 

j MF^{(3^{^°)+a)dqdx 

^ I M 

< / M 

= "(1 + 777) / Mdqdx+ I MJ'{P^{^°) + a)dqdx 
<^a\n\+ [ MJ'{l3^{$") + a)dqdx. 

Thus we have shown that 

T5{a) < la\fl\ + [ MJ'{P^{$°) + a)dqdx. (4.25) 
2 JnxD 

Now, there are two possibilities: 

Case 1. If + a < 1, then < /3^(V'°) < 1 - a. Since L > 1 it follows that < /3^(s) < 1 

if, and only if, /3^(s) = s. Thus we deduce that in this case /3^(7/''^) = V'^i ^-iid therefore 
< + a) = + a). 

Case 2. Alternatively, if 1 < + a, then, on noting that I3^{s) < s for all s e [0, 00), it follows 

that 1 < /?^(^'') + a < 'tp'^ + a. However the function is strictly monotonic increasing on 
the interval [1, 00), which then implies that = J'(l) < F{l3^{%p°) + a) < + a). 

The conclusion we draw is that, either way, 

< J'(/3^(V^°) + a)< + a). 

Hence, 

T5{a) <^a\n\+ [ MJ^i^P" + a)dqdx. (4.26) 
2 JnxD 



(/3^(V^°)+a)2-L2 



2L 



+ (/3^(V")+a)(logL-l) + l 



dq dx 



{L + aY 
2L 



+ (/3^(f°) + a)(log(/3^(V^°) + a) - 1) + 1 



dq dx 
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Substituting (|4.26p into (|4.24p thus yields, for i = t„ and n e {1, . . . , A^}, 



1 





2\kC 



M 



Jo ~ ~ 



2fc / Af^(V^f (<) + a) dg dx + — - / / M(V^f - )2 dq dx ds 

2k e / M ~ d(jda;ds+ 77V / / ^'^ X., , dgdxds 



/ ; At,+ , ^o. I . . 4. 



< ho||' + 



-/ ||/'^*^+(s)||^,ds + 3afc|f7| +2fc /" MJ-(V^" + a) dq dx. (4.27) 

The key observation at this point is that the right-hand side of (|4.27|) is completely independent 
of the cut-off parameter L. 

We shall tidy up the bound (|4.27p by passing to the limit a ^ 0+. The first a-dependent 
term on the right-hand side of (14.271) trivially converges to as a — >■ 0+; concerning the second 
a-dependent term, Lebesgue's dominated convergence theorem implies that 

lim / MJ^i^P^ + a)dqdx^ [ Af dg dx. 

Similarly, we can easily pass to the limit on the left-hand side of (|4.27p . By applying Fatou's 
lemma to the fourth, sixth and seventh term on the left-hand side of (I4.27|) we get, for t = tn, 
11 e {1,...,A^}, that 

lim inf„^o+ / MF{i^f'j]+{t) + a) > [ M^(vi^f['+(i)) dq dx, 
JnxD ' JnxD 



lim info,. 



JnxD 



M ^ ■ ' dq dx ds 



> 



M 



JnxD 



dq dx ds 



JnxD 



' dq dx ds, 



lim infQ,_s.o_|_ / / M _/ dgdxds > M ~ — dgdajds 

Jq JnxD t/j^ £ + a ~ ~ Ja JnxD £ ~ ~ 



= 4. j I M\^qJijf^£+\^ dqdxds. 



JnxD 



Thus, after passage to the limit a — > 0+, on recalling p.22p . we have, for alH = i„, n e {1, . . . , N^, 
that 



iiy.rwii^ + xT u 



+ 2k / MT{i^fl:+{t)) dqdx + —- / Af (^f - V^f dq dx ds 
JnxD ' ~ ^r^Jo JnxD 

t 



M\Yx\/^t'£^\''dq dxds 



'0 JnxD 

+ ^ /7 M\V,J^fl^\'dqdxds 
^ Jo JnxD 
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<ll«o||' + - f \\f^''+is)\\l,,ds + 2k f M^(^")dgds 
Jo ~ JnxD 

<llyo||' + - r \\f{s)\\l,ds + 2k I M^(7Ao)dgdx=: [B(yo,/,V^o)]' 



(4.28a) 



(4.28b) 



SlxD 



where, in the last hne, we used p.22p to bound the third term in (|4.28ap . and that t G [0,T] and 
the definition (I3.24p of /^*'^ to bound the second term. 

We select ip = X[o.t] C^l^ a-s test function in (|4.4p . with t chosen as i„, n £ {1, . . . , N}. Then, 
similarly to (|4.12p . we deduce, with i = t„, that 



iicf^rwii' + TTT / iicr(^)-cr(s)ii'ds+2£ / iiv<r(s)fd 



At Jo 



At,- 



< 



r(^")dg||'< (4.29) 



where we have noted ((05)) . ([332]) and that /3^(V'°) € Z2. 

Next, we develop L-independcnt bounds on the time-derivatives of Tp^j^ and C^}^- 

4.2 L-independent bounds on the time-derivatives 

We begin by bounding the time-derivative of -^^^ using (j4.28bp : we shall then bound the time- 
derivative of u^l^ in a similar manner. 

4.2.1 L-independent bound on the time-derivative of -^^^ 

It follows from (|4.3bp that 



I M — --— Lp dq dx dt 
JOxD dt ^ ~ 







< 






JO JnxD 



M Vx^f*!^^ ■ ^xpdq dx dt 



JnxD 

T 



M uf'l- (V^^['+) • dq dx dt 



K K 



Txl I E E ^<?. ^e""!' + ■ V,. V dq dx dt 

^^Jo JnxD i=i j=i ~ ~ ^ ~ 

/ / M Y.[a{ufl+) q. ] /3^(^fl'+) • V,,^ d<7 dx dt 
Jo JnxD ,_i ~ ~ ' ~ ' ~ ~ ~ 



=:Si+S2 + S3 + S4 ypeL\0,T;X) 



(4.30) 



We proceed to bound each of the terms Si, . . . , S4, bearing in mind (cf. the last sentence in 
the statement of Lemma 13.31) that 



V'f 1'^ > a.e. on n X D X [0,T], 



Miq)dq^l, 



0< / A/(q)V'f['+(x,g,t)d(7 < 1 for a.e. (x,t) e X D. 



(4.31a) 
(4.31b) 



D 

In addition, we make use of the following result 

Lemma 4.1 For i ^ 1, . . . , K and for allt — tn, n = 1, . . . ,N, 

1? 



nxD 



M$f '£^{t)dqdx < — 

' ^ Cil 

1 



nxD 



MT{^fl^{t))dqdx+\n\ I Me 



< 



[B(uo,/,V'o)]' + a,< 



dq 
(4.32) 
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where and cn are as defined in (|1.4ap . and 

Proof First we recall the logarithmic Young's inequality 

r s < r log r ~~ r + e'^ for all r, s E IR>o. (4.33) 
This follows from the Fenchel- Young inequality: 

r s < g*{r) + g{s) for all r, s G R, 

involving the convex function g : s G M i— > g{s) G (— oo,+cx)] and its convex conjugate g*, with 
17(5) = e* and 

{+00 if r < 0; 

if r = 0; 

r (logr — 1) if r > 0, 

with the resulting inequality then restricted to M>o- It immediately follows from (|4.33p that 
r s < J-{r) + e* for all s E K>o. Hence we have that 



Oil 



2 JnxD 2 ^ ^ ~ jnxD 



M 



dq dx. (4.34) 



The desired result (|4.32p follows immediately from (|4.34l) . on noting (j4.28bp and (|1.4cp . □ 
We shall use throughout the rest of this section tests function if such that 

^ e L^{0,T]W^'°°inx D)). (4.35) 
We begin by considering Si, noting (|4.31a[ b'): 

i-T 



Si = 2e 



< 2e 



JnxD 

T 



M y,VV^fi'+ • V^^dq dxdt 



Jn 



M 



D 



?At.- 



d£ lly:r(^|lLoo(£,) 



dxdt 



< J-\8ke 



M 



dq dx dt 



lo JflxD 

Hence, by (I4.28bp with t = tN = T, 

lo Jfl 

Next, on noting (I4.31al b') and (I4.9p . we consider term S2 

At,- ||2 



y2;^|li=o(D)dx dd 



S2 < 



Si < \/— B(yo,/,V'o) ( / / iiy.^||io.p)dxdt 



lyf"!'"! / M/3^(V^i'r)dg ) \\Y.ip\\L^iD)dxdt 



(4.36) 



/o "'n 



< 



u 







/o Jo 



y2:^lli«=(£,) dx dt 



(4.37) 



31 



where Cp{ft) denotes the (positive) constant appearing in the Poincare inequahty < Cp(il) 
\\Vx v\\ ^ foi' ^'^y ?^ S F C Hl{Q). On recaUing the definitions of uf*£^ from (j4.1bp . and noting 
(feO)) and (|4.28bl) we have that 



T-At 



ly.y^rip di 



< llyoir+ / WY.ufl^W'dt 



< 



(l + i) [B(uo,/,^o)]' 



Therefore, 



S2 < Cpin) {l + iy B(yo, /, / lly.^lli^p) dx dt 



^Sl 



■) 



(4.38) 



(4.39) 



Alternatively, without the use of the Poincare inequality, directly from the second line of (|4.37p . 
we have that 



S2 < Vf 



At, 



eSS.SUPtgfo^T] WUeX 



Jn 



(4.40) 



Similarly as above, 

ess.suptg[o -rj||w^^"(i)||2 = max (^||uo|p,ess.sup4g(o -r-AtjIk^L^WI 

< max (1 1 uo II ^ ess.suptg [o,T] W^tl'^ (*) I 

< [B(zio,/,^o)]'. 

Combining (|4?39l) . (|440)) and (|44T|) . we have that 

S2<min(Cp(r!) (l + i)^Vr)B(uo,/,V^o)( / / || V.^H^ dx di 



(4.41) 



(4.42) 



Jn 



We are ready to consider S3; we have that 
'0 JnxD 



S3 = 



1 

2~A 



i=i j=i 



< 



/ / ^^2v/^EE^..y,^y^.y,,^dgd3;di 



12. 



/o JnxD 



/o JnxD 



K 



lYg^^l"^ dgdjdi 



/o JnxD 

i-T 



(^ElY^'^l'j dgdxdt 
lYg^l dq dx dt 



Jo \Jd 



M^f'j:+dq 



M 



d 



?At,- 



dq \\Yq'f\\L-^{D} dxdt 



< 



\A\ 2aok 



■\/2aokX \ ^ J 



M 



nxD 



V 2 / />T 



dqdxdt] I / / WYqipWl^^jj-^dxdt 



Jn 
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Thus, by (l4:28bl) . 



S, < 



1^1 



-\/2aofcA 



B(wo,/,'0o) 



Jo 



(4.43) 



Finally, for term S4, recalling the inequality fi^{s) < s for s e M>o, (|4.31al b') and (|4.32p , we 
have that 



Jo JnxD ~ ~ ' ^ ~ ^ ~ 

M\q\ |a(<i'+)|/3^(^,^i'+) |V,^| dgdxdt 



< 



< 



< 



< 



< 



/o JnxD 



l'^(<L^)l ( / Mk|/3nCr)dg) ||V,^|U^(,,)dxdi 



/o Jsi 



jL /TAt 



JSl 
T 



2 TAt,- 



AfV^^^[^+dqj ||Vg(?|U^(^)d^di 



D 



.1 



l^(^^L^)l^da;dij |^ess.suptg[o,T] ^ M \q\'^ ijj^]^^ (t) dq dx 



"'O 



1 2 „'tAt,+ / 



< 



\ 



2 



i=l 
T 



kcii 



[BK,/,VoP + C,, 







lY9^lli~(nxD)dt , 



if 



where we have noted, on applying a Young's inequality, that |gp < ■^[{"d ~ 1)K + y~^(^|gtP)^]- 

1=1 

Hence, by (I4.28bp . 



S4 < 



2 



K r 



il 



-[B(uo,/,^oP + a, 



cxp 



B(uo,/, V'o) 



(4.44) 



Upon substituting the bounds on Si to S4 into (li^Ul) . with ^ e L'^{Q,T]W^^°°{n x D)), and 
noting that the latter space is contained in L^{Q^T\X) we deduce from (|4.30p that 



'0 JnxD 



M —^Sdq dx dt 
at ^ ~ 



<C, max [B(wo,/,V'o)]^B(wo,/,V'o) / W\\w^,^inxD)<^t ] (4-45) 
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for any e L^{0,T; W^'°°{il x D)), where C* denotes a positive constant (that can be computed 
by tracking the constants in (I4.36l) - (|4.44l) ). which depends solely on e, v, Cp(f2), T, \A\^ aq, fc, A, 
M, K and {c,i}f£i. 

We now consider the time derivative of It follows from (|4.4p . (|4.29p and (j4.4ip that 



Jo. 



< 



dt 



if) Ax At 



< 



Jn 



At.- 



At,- /-At,-i 



+ess.sup,e[o,T]llC'+llL==(o) 



Ax At 



< 



e(H) +B(uo,/,V'o) 



|V:,^||2dt 



l|V.(^f dtj 
e L\0,T;H\n)). 



(4.46) 



4.2.2 i-independent bound on the time-derivative of uf]^ 

In this section we shall derive an L- independent bound on the time-derivative of u^^. Our starting 
point is (I4.3ap . from which we deduce that 



'0 Jn 



< 



dt 



Jn 

T 



• w dec dt 



At,+ 

6,L 



{f^''+,w)vAt 

Ui + U2 + U3 + U4 







+ k 



w Ax At 
Jn 



T 

Jn 



V:r u^i^ : w Ax At 



V / I C, (M V^f -.^xwAx At 

,_i Jo Jn~ ' » ~ ~ 

e L\0,T;V). 



(4.47) 



On recalling from the discussion following (|3.3I) the definition of Fcr, we shall assume henceforth 



that 

Clearly with this choice of a, L'^{0,T;Va) C L^{0,T;V) and we will exploit the embedding 
H-^in) W^^°°{n). Using (|4;4T]) and (|4~28bp . we have 



ui < ess.sup,g[o.T]iiMf:rii I / iiy.w,^i^+f dtj (y 



lji'llL~(a) dt 



< 



"B(uo,/,i/'o)] 



L°=(n) 



dt 



(4.48) 



For term U2 we have, 



U2 < V^U J^\\Yxufl+fAt\ U iiy. 



wfAt 



< V^BiuoJ,$o) / WYxwfAt 



(4.49) 
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Concerning the term U3, on noting the definition of the norm || • and that thanks to (|3.24p we 
have 

ll/'^*'^llL2(o,T;y') < II/IIl2(o, T;y')' 

it follows that 



U3 < V^BiuoJ,$o) / WY.wfdt 



(4.50) 



We now bound the term U4. On noting (jl.4bp . Young's inequality, (|4.31b[) and (|4.32p . we have 
that 



U4 = fc 



K 







dx dt 









2kf: r 1 1 M 

Jo Jn Jd 



< 



< 



Ci2 akil + Cj3 2l9»l 



ipf*£~^ I Vj; wldqdx dt 



dx 



dt 



< 2k 



^Ec»2 + V(^+C,3) T— [BK,/,V'0)]' + a,exp / ||V, 



w|lioo(f2) dt. 

(4.51) 



Collecting the bounds on the terms Ui to U4 and inserting them into (|4.47l) yields 



T , dufX 



Jn 



dt 



• w dx dt 



<C,, max [B(mo,/,V'o)]',B(uo,/>o),1 / \\w\\ldt\ (4.52) 



for any w G i'^(0, T; V^o-), cr > l + jrf, where C** denotes a positive constant (that can be computed 
by tracking the constants in (|4.48p - (|4.5ip ). which depends solely on SI, d, k, X, uq, M, d, K, 
{c.2}^£i and {c.3}^i. 



5 Dubinskii's compactness theorem 

Having developed a collection of L- independent bounds in Sections 14.11 and 14.21 we now describe 
the theoretical tool that will be used to set up a weak compactness argument using these bounds, 
— Dubinskii's compactness theorem in seminormed sets. 

Let ^ be a linear space over the field M of real numbers, and suppose that is a subset of A 
such that 

(V^ G X) (Vc e R>o) opeM. (5.1) 

In other words, whenever ip is contained in A4, the ray through ip from the origin of the linear 
space A is also contained in A4 . Note in particular that while any set A4 with property (|5.ip must 
contain the zero element of the linear space A, the set J\4 need not be closed under summation. 
The linear space A will be referred to as the ambient space for M . 

Suppose further that each element of a set Ai with property (|5.1I) is assigned a certain real 
number, denoted [95]^^; such that: 

(i) [(p]m > 0; and [pIm = if, and only if, (p = 0; and 

(ii) (Vc e M>o) [cp>]m = c [p]m- 
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We shall then say that is a seminormed set. 

A subset S of a seminormed set Ai is said to be bounded if there exists a positive constant Kq 
such that [ip]M < Kq for all (p & B. 

A seminormed set Ai contained in a normed linear space A with norm || • ||^ is said to be 
embedded in A, and we write A4 ^ A, if: 

(3i^o G M>o) (V^ e X) MA<Ko[ip]M- 

Thus, bounded subsets of a seminormed set are also bounded subsets of the ambient normed linear 
space the seminormed set is embedded in. 

The embedding of a seminormed set Ai into a normed linear space A is said to be compact if 
from any bounded, infinite set of elements of A4 one can extract a subsequence that converges in 
A; we shall write A4 '-H- A to denote that A4 is compactly embedded in A. 

Suppose that T is a positive real number, if maps the nonempty closed interval [0, T] into 
a seminormed set Ai, and p E R, p > 1. We denote by LP{0,T; Ai) the set of all functions 
(p : t €[0,T]^ ip{t) e M such that 



i/p 



Mt)rMdt 



< oo. 



The set LP{0,T; Ai) is then a seminormed set in the ambient linear space LP{0,T; A), equipped 
with 



'\LP{0,T;M} 



We shall denote by L°°{0,T; Al) and ['p]l'^{o,T:M) the usual modifications of these definitions 
when p = oo. The following theorem is due to Dubinskii |18) . 

Theorem 5.1 Suppose that Aq and Ai are normed linear spaces, Aq ^ Ai, and At is a semi- 
normed subset of Ao such that Ai ^-H- Aq. Consider the set 



y 



ip : [0,T] ^ M : [p]lp(o,T:M 



dp 



dt 



< oo 



LPi{0,T-Ai) 



where 1 < p < co, I < pi < oo, ||-mi is the norm of Ai, and dp/dt is understood in the sense 
of Ai-valued distributions on the open interval (0,T). 
Then, y is a seminormed set with seminorm 



'\y 



LP{a,T-M) 



dp 



dt 



LPi(0,T;Ai) 

in the ambient linear space LP{0,T;Ao) D W^^p^{0,T;Ai), and y ^ LP{0,T;Ao). 
In the next section, we shall apply Dubinskii's theorem by selecting 
Aq — i(i+|(7|)2fljv,/(^ ^ ^) 

M{q){l + \q\f''\p{x,q)\dxdq<^ 



■.= \peLl,{nxD) : \\p\U := / . 



and 



M 



p G Aq '■ P > with 

^lig) ( Yx^J'?ix,q) + Vq^p{x,q) 



flxD 



dxdq < oo 
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and, for (p £ A4, we define 



• M 



dx dq. 



Note tliat is a seminormed subset of the ambient normed linear space ^o- Finally, we put 

Ai := M-^W{n X D)' -.^ {ip : M(p eH'{nx £>)'}, 
equipped with the norm 

and take s > 1 + ^{K + l)d. With such s it then follows from the Sobolev embedding theorem on 

ll^ll \{Mp,x)\ ^ W\\Li,{QxD)\\x\\L^inxD) 
mlA^ = sup — < sup T—r 

X&H'{nxD)\\X\\H=(nxD) x^H^iUxD) II XII H= (f2x_D) 

< Ko\\p\\Ll^(nxD) < -ft^oll^lUo, 

where Kq is any positive constant that is greater than or equal to the constant K^, the norm of 
the continuous linear operator corresponding to the Sobolev embedding {H^{Q x D) ^)H''^^{Q x 
D) ^ L°°(f7 X D), s > 1 + ^{K + l)d. Hence, we have that Ao ^ Ai. 

Trivially, M ^ Aq. Wc shall show that in fact M ^-H' ^q. Suppose to this end that B is an 
infinite, bounded subset of A4. We can assume without loss of generality that B is the infinite 
sequence {pn}n>i C Ai with [(^njvvt < Kq for all n > 1, where Kq is a fixed positive constant. 
We define pn := and note that p„ > and pn € Hl,i{Q x D) D L'^^^^^^^^ x D) for all 

71 > 1, with 

Since H\j{^l x D) Ci Lj^^^q^-^2t) ^ ^) compactly embedded in L^^^|g|^2o^,/(i^ x D) (see 

[Appendix F| at the end of the paper for a proof of this), we deduce that the sequence {p„}„>i 
has a subsequence {p„j.}fc>i that is convergent in x D); denote the limit of this 

subsequence by p; p G L^-|^^|^|^2fljv,/(^ x D). Then, since a subsequence of the sequence {pnk}k>i 

also converges to p a.e. on x D and each p„j. is nonnegative on fl x D, the same is true of p. 
Now, define ip := p^, and note that ip G L^^^|g|^2i)jv/(f^ x D). Clearly, 

Wn^-'fU^ „ 29 inxD)= Mil + lqlf"^ (pn^+p)\pn^-p\dxdq 



<\\Pnk+P\\L^ „ (OxD) II - P IIl^ (nxD) 

< llPnfclU^ „ (axD) + II pIIl^ ^ (nxD)] 

X II Pnfc - P||l2 „ 2 5 (f2xD)- 

As {pnk}k>i converges to p in L^^^|^|^2t)^/(^^ x £)), and is therefore also a bounded sequence 
in L^-^^|^|j2t)^/(f^ X D), it follows from the last inequality that {(pnk}k>i converges to ip in 
^}i+\q\)^^ X D) ~ Aq. This in turn implies that the seminormed set A4 is compactly embedded 
in ^0 • Thus we have shown that the triple A4 Aq =-> Ai satisfies the conditions of Theorem 
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Remark 5.1 There is a deep connection between A4 and the set of functions with finite relative 
entropy on D, exhibited by the logarithmic Sobolev inequality: 



Af(g)|p(<z)|2log 



\m\' 



dq<- f M{q)\Y,m\^dq e HIj{D), 

~ 1^ J D ~ 



(5.2) 



with log-Sobolev constant k > 0; the inequahty (|5.2p is known to hold whenever M satisfies the 
Bakry-Emery condition: Hess(— log > k Id (in the sense of symmetric Kd x Kd matrices) 

on D — DiX ■ ■ ■ X Dk, asserting the logarithmic concavity of the Maxwellian on D. The inequality 
(|5.2[) follows from inequality (1.3) in Arnold, Bartier & Dolbeault [3]. 

The validity of the Bakry-Emery condition for the Hookean Maxwellian, Ui{s) — s for i — 
1, . . . ,K, is an easy consequence of the fact that 



Hess(-logM(g)) = Hess (j^U, 

= diag (hbss (^[/i(i|gi|2)^ 



,Hess[UK{l\qK\^)]] =ld. 



(5.3) 



. , qx) E Di X ■ ■ ■ X Dk — D, so k — 1. If we replace Ui{s) = s by the potential in 
. ,K, it is easy to show, on noting that dct(/ + a fy^) = 1 + a • 6 for all a, b E R'^, 



for all g = (q i , 
(lO)) for i = 1, 
that AC > 1. 

More generally, we see from ()5.3|) that if qi E Di i-^ Ui{^\qi\'^) is strongly convex (see, e.g., 
Hiriart-Urruty & Lemarechal [12], p. 73) on Di with modulus of convexity Ki > 0, i — 1, . . . ,K, 
then M satisfies the Bakry-Emery condition on D with n — min{Ki, . . . , kk}- 

On writing ^p{q) := |p(g)P ( > 0) in (|5.2p . we have that 



M{q)ip{q)log 



D 



^(q) , 2 

— '^1 - ~ 



M{q) 



D 



d(j. 



(5.4) 



for all tp such that > on and y/if e H\[[D). Taking if — {p/M where is a probability 
density function on Z?, we have that ||<^||ii (jj) = ||</'||li(_d) = 1; thus, on denoting by v the Gibbs 
measure, defined by dv = M{q)dq, the left-hand side of (|5.4p becomes 



S{cp\M) 



M 



referred to as the relative entropy of ip with respect to M. The expression appearing on the 
right-hand side of (|5.4p is 1/(2k) times the Fisher information I{(p) of if, where 



nip) := E 



Vg log ip{q) 



Vqlogifiq) ip{q)dv 



D 



dv. 



where, E is the expectation with respect to the Gibbs measure. o 
The following simple lemma will be helpful in the next section. 

Lemma 5.1 Suppose that a sequence {'Pn\^=i converges in L^(Q,T] L^^-^^i^q^y^ x D)) to a 
function ip G ^-'-(0, T; i|--|^^|^|j2tfjy,_f (i^ x D)), and is hounded in L°°(0, T; L^j^_|_|g|pfl^^(fi x D)), 
i.e. there exists Kq > such that \\p>n\\L'^(Q t-l'^ (nxD)) — ^^o for all n > 1. Then, 

ip E L^(0, T; L^j^^|q|^2fl^f (f^ X D)) for all p £ [l,oo), and the sequence {i^Ti}ri>i converges to ip in 
LP{0,T;Ll^^q^y,,ljin X D)) for all p e [l,oo). 
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Proof Since {(^„}„>i converges in L^{0,T;L^^_^^q^^2ii]^j{^ x D)), it follows that it is a Cauchy 
sequence in L^{0, T; L^^^^q^y2^ x D)); thus, for any p e [1, oo), there exists uq — nQ{e,p) G N 
such that for all m,n> rio(e,p) we have 



(1+|(7|)2i5m 



Hence, for all m,n> no{e,p), 



i/p 



(i+igi)2fM^ ' 



< eSS.SUPtg[o.T]ll^n 



(i+igi)2fjf^ ' 



1/p 



\\^n - 'fvaWh^ (nxD) 



< e. 



This in turn implies that {'fn}n>i is a Cauchy sequence in the function space LP{0, T; L^i^^q^y^ m 
{ft X I?)), for eachp e [l,oo). Since LP(0, T; Ljj^_|_|g|pi,^^(il x D j) is complete, {^„}„>i converges 
in LP(0, T; Lj-^^|g|^2tf jv/(^ ^ -^)) ^'^ ^ limit, which we denote by <^(p), say. Since, by assump- 
tion, {^„}„>i converges in £^(0, T; Lj^^|^|)2i»M(^ ^ -^))' ^^(^'^'Hi+m^^Aii^ ^ ^)) ^ 
L^(0, T; L^-^^^q^~^2i) jy.ji^ x ^)) for each p G [1, oo), it follows by uniqueness of the limit that ^(p) = if 
for aU p e [1, oo). Hence, also, ip G ^^(0, T; L|^^|^|^2«m(^ ^ -^)) '^l^ P ^ t-*^' °^)- '-' 

6 Passage to the limit L — > oo: existence of weak solutions 
to Hookean-type bead-spring chain models with centre- 
of-mass diffusion 



The bounds (|4.28b[) , (|4.45p and (|4.52p imply the existence of a positive constant C* , which depends 
only on B(uo, fjipo) E^nd the constants and C**, which in turn depend only on e, v, Cp(r2), T, 
1^1, flo, fc. A, SI, d, K and Af , but not on L or At, such that: 



A^.^ 



ess.suptg[o_j.]|ju^_;f; ''{t)\ 
+ ess.suptgfo^T] 



1 

At 



At,- 



nxD 



AtL 



JflxD 



JnxD 

T 



+ 







5<i 



2 



^e^i'^r '^'^ '^'^ 



JnxD 



^e,L I 



dq dx ds 



dt 



dt 



V' 



M 



e.L 



dt 



dt < a. 



3.1) 



H={nxD)' 



where || • \\v^ denotes the norm of the dual space V'^ of Va with a > 1 + ^d, (cf. the paragraph 
following (|4.47p ): and || • \\H={nxD)' is the norm of the dual space i/*(Sl x D)' of i/*(Sl x D), with 
s > 1 + ^{K + l)d. The bounds on the time-derivatives stated in the last line of (|6.ip follow 
directly from (|4?52)) . and (|445t using the Sobolev embedding H^fl x D) ^ x D). 

By virtue of (14.411) . (|4.38p . the definitions (|4.1al b). and with an argument completely analogous 
to (|4.38p on noting (|4.9p , (I3.22p and (13.51) in the case of the fourth term in (|6.ip , we have (with a 
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possible adjustment of the constant C*, if necessary,) that 



1 /"^ 

ess-supigfo^Tillw^L'^^^Wir + A< io 



At.+ At,-||2 
K.L -^e,L II 



Jo ~ ~ 



ess.suptg[oy] / MJ'{^f^^''^\t)) dq dx + 



nxD 



AtL 



JnxD 



Mi^fl:+ -$f^'£~fdqdxds 



M\y^Mf'£+\^dqdxds+ / M\\7gMfl;+f dqdxds 



at 



dt- 



M 



9i 



dt < a. 



On noting (|4.31al b'). (|lAa[ bV and (I53)) . we also have that 

^M,±) ^ p on X Z? X [0, T] 



and 



M{q)^if^j}''^\x,q,t)dq < 1 for a.e. {x,t) Gftx [0,T]. 



Henceforth, we shall assume that 



At = o{L' 



as L 



3.2) 



(6.3) 



(6.4) 



(6.5) 



Requiring, for example, that < At < Co/{L logL), L > 1, with an arbitrary (but fixed) constant 
Co will suffice to ensure that (|6.5I) holds. The sequences 



L>1 



as well as all sequences of spatial and temporal derivatives of the entries of these two sequences, 
will thus be, indirectly, indexed by L alone, although for reasons of consistency with our previous 
notation we shall not introduce new, compressed, notation with At omitted from the superscripts. 
Instead, whenever L ^ cx) in the rest of this section, it will be understood that At tends to 
according to (|6.5p . 

We are now almost ready to pass to the limit with L oo. Before doing so, however, we first 
need to state the definition of the function ip'^ that obeys (|3.22p . for a given -00 satisfying p.Sp . 
We emphasize that up to this point we simply accepted without proof the existence of a function 
■0° obeying (jX^ for a given ■00 ■ The reason we have been evading to state the precise choice 
of ■0° was for the sake of clarity of exposition. The definition of and the verification of the 
properties listed under p.22p rely on mathematical tools that were not in place at the start of 
Section[3]where the notation ■0° was introduced, but were developed later, in the last two sections. 
The details of 'lifting' 0o hi^to a 'smoother' function are technical; they are discussed in the 
next subsection. 

A second remark is in order. One might wonder whether one could simply choose 0^^ as 0^0; 
indeed, with such a choice all of the properties listed in (|3.5p would be automatically satisfied, bar 
one: there is no guarantee that [0o]^^^ G Hj^^flx D). Although the property [0^"]^/^ e H\j{Q'><D) 
has not yet been used, it will play a crucial role in our passage to the limit with L oo in Section 
16.21 In fact, in the light of the logarithmic Sobolev inequality (15. 4p . on comparing the requirements 
on 0^0 in p.5l) with those on 0° in p. 221) . one can clearly see that the role of the condition 
|-^oji/2 g H\j{Vl X D) in (13.221) is to 'lift' the initial datum ^00 with finite relative entropy into a 
'smoother' initial datum 0° that also has finite Fisher information, in analogy with the process 
of 'lifting' the initial velocity uq from H into m° in V . That the choice of 0o as ■0° is not a good 
one can be seen by noting the mismatch between the third term in (|6.2p arising from the Navier- 
Stokes equation on the one hand, and the sixth and seventh term in (|6.2p that stem from the 
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Fokker-Planck equation. The absence of bounds at this stage on -0^*' and -ip^* in those terms 
is entirely due to the fact that, to derive (|6.2p . we did not use that [-0"]^/^ g Hlj{fl x D). This 
shortcoming of (|6.2p wiU be rectified as soon as we have defined '0" and shown that it possesses 
all of the properties listed in (I3.22p . 

6.1 The definition of ^/^° 

Given 0o satisfying the conditions in p.Sp and A > 1, we consider the following discrete-in-time 
problem in weak form: find C^'^ G H\^{VL x D) such that 

/ M- --^ ^dgda;+ / M V^^C^'^ • V:^^ + V^C^^'^ • Vg^ dg da; = (6.6) 

for all if e Hli{n X D), with C^'^ P^ii^o) & LIj{VL x D). Here {3^ is defined by (fTTSl) . with 
L replaced by A. The function which we shall encounter below, is defined by ()4.6p . with L 
replaced by A. ^ 

The existence of a unique solution Q^'^ £ x D) to (|6.6p . for each At > and A > 1, 

follows immediately by applying the Lax-Milgram theorem. The parameter A plays an analogous 
role to the cut-off parameter L; however since we shall let A — > oo in this subsection while, for 
the moment at least, the parameter L is kept fixed, we had to use a symbol other than L in 
(|6.6p in order to avoid confusion; we chose the letter A for this purpose in order to emphasize the 
connection with L. 

Lemma 6.1 Let C,^'^ be defined by (|6.6p . and consider --y^'" defined by 

7^'"fe) := / Af(g)C^'"fe,g)dg, n = 0, 1. (6.7) 

Jd ~ ~ ~ 

Then, (^^'^ is nonnegative a.e. on nx D, andO< j^''^ < 1 a.e. on fl. 

Proof The proof of nonnegativity of C^'^ is straightforward (cf. the discussion following (|3.67p ). 
We have that [C^^°]- = a.e. on fl x D, thanks to (^3)) and the definition of l3^. We take 
f — [C^' ^ ] - as a test function in (|6.6p , noting that this is a legitimate choice since C^'^ G H\j (O x D) 
and therefore [C^'^]_ € Hli{Vt x D) also (cf. Lemma 3.3 in Barrett, Schwab & Siili [7]). On 
decomposing C,^'^ — [C^'^]+ + [C^'^]-, and using that 

[^A,i]^ [^A,i]_ ^ 0, y,[c^-i]+.y,[c^^i]_ = o and yjc^'i]+.y,[c^-i]_ = o 

a.e. on O X £), we deduce that 

= X7 / ^""'^ t^""'']- dx = ^J M [C'^'°]+ [C^'i]_ dq dx < 0, 

JnxD ~ JOxD 

where || • || denotes the L^{il x D) norm. This then implies that 

[fA.l]_||2 < 0. 

Hence, [C^'"'^]_ = a.e. on J7 x _D. In other words, C^'^ > a.e. on £7 x £), as claimed. 

In order to prove the upper bound in the statement of the lemma, we proceed as follows. 
With 7^'" as defined in (|6.7p . we deduce from the definition of C,^'^ and Fubini's theorem that 
^A,i ^ H^{V,). Furthermore, on selecting (p — ip £ H^{^1) ® l(-D) in (|6.6p . recall (I3.65p . we have 
that 

/ 1^'^ -1^'° f V:,7'^'i- V^(pdx = \lLpeH^{n). (6.8) 

At ~ Jn~ ~ ~ 
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As C^^" = 13^ i^o), and < /3^(s) < s for all s G R>o, we also have by (|33]) that 

0<7'^'°= / MP^{ipo)dq< I MiPodq^l on n. (6.9) 



D J D 



Consider 

^A,n _ _ ^A,n^ n = 0, 1. 

On substitutmg 7^^^" = 1 — z'^'", n = 0, 1, into (|6.8p . we have that 

/ T- V3da;+ / V^z^'^ ■ V:^ipdx \/ipeH^{VL). (6.10) 

Also, by (16. 9p . we have that < z^'^ < 1. By using an identical procedure to the one in the first 
part of the proof, we then deduce that [z^'^]_ = a.e. on n. Thus, z^'i > a.e. on which 
then implies that 7^^'^ < 1 a.e. on f2, as claimed. □ 

Next, we shall pass to the limit A — > cxi; as we shall see in the final part of Lemma [6.21 below, 
this will require the use of smoother test functions in problem (j6.6p , as otherwise the term involving 
^A,o _ /?^(^q) is not defined in the limit. In any case, our objective is to use the limit of the 
sequence {C^'^}a>i, once it has been shown to exist, as our definition of the function ■0°. We shall 
then show that ■0° thus defined has all the properties listed in (|3.22p . 

To this end, we need to derive A-indepcndent bounds on norms of C^'^, very similar to the L- 
independent bounds discussed in Section 21 Since the argument is almost identical to (but simpler 
than) the one there (viz. (|6.6p can be viewed as a special case of (j3.23bP . with /", ^^'^ HeL 

taken to be identically zero, X = ^, e = 1, N ^ 1, and A chosen as the K x K identity matrix), 

we shall not include the details here. It suffices to say that, on testing (|6.6I) with J^'{(^'^ + a) and 
passing to the limit a 0^, analogously as in the derivation of (I4.28ap in Section |4l we obtain 
that 

f MT(C^'^)dqdx+ 4At f M \\'^\/^\^ dqdx + AAt f Af | V, \/^| ^ dg dx 

JuxD ~ ~ JnxD ~ ~ ~ JnxD ~ ~ ~ 

< [ MJ"(V^o) dg da;. (6.11) 

Our passage to the limit A — > cxi in (|6.6p is based on a weak-compactness argument, using (16. lip , 
and is discussed below. ^ 

We have from Lemma [6.11 that {[C^'"'^] ^ }a>i is a bounded sequence in L\.j{il, x D). Using 
this in conjunction with the second and third bound in (|6.1ip we deduce that, for At > fixed, 
{[C'^'^]^}a>i is a bounded sequence in Hl^j{il x D). Thanks to the compact embedding of H\,^{Vl x 
D) into L\,i(p, X D) (cf. [Appendix FP , we deduce that {[C^'^]Ha>i has a strongly convergent 
subsequence in L\[{Q. x D), whose limit we label by Z, and we then let Z^. 

For future reference we note that, upon extraction of a subsequence (not indicated), C^^'^ then 
converges to C,^ a.e. on x D; and C^'^{x, ■) converges to C,^{x, ■) a.e. on D, for a.e. x ^VL. 

By definition, we have that > 0; furthermore, thanks to the upper bound on 7"^'^ stated in 
Lemma |6. 11 the remark in the previous paragraph, and Fatou's lemma, we also have that 

[ M{q)(\x,q)dq <1 for a.e. x G 17. (6.12) 

Jd ~ ~ ~ ~ 

Further, again as a direct consequence of the definition of C^, we have that 



■a,i^./^i strongly in (17 X D). (6.13) 
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Application of the factorization ci — C2 = (-v/ci ~ a/^J) (\/^^" ^i^h ci, C2 G M>o, the Cauchy- 
Schwarz inequahty and (j6.13l) . then yields that 

C^'^ strongly in Llj{n x D). (6.14) 

Finally, we define 

^° (6.15) 
It follows from the nonnegativity of and ()6.12p that 

<>0 a.e. on X L> and 0< / M{q)ti;°{x,q)dq < 1 for a.e. x G rj. (6.16) 

Jd ~ ~ ~ ~ 

Further, from the bound on the first term in (|6.1ip and Fatou's lemma, together with the fact 
that, thanks to the continuity of J^, (a subsequence, not indicated, of) {J^{C^'^)}a>o converges to 
J-{C^) — J'iw) a-e. on f2 X Z?, we also have that 

/ MF{^°)dqAx< I MJ"(?/^o) dg dx. (6.17) 
Next, we note that from (|6.13p we have that, as A — > cxo. 



M2 y'^-a,! ^ \JC strongly in L^{n x D). (6.18) 

We shall use (|6.18l) to deduce weak convergence of the sequences of x and q gradients of d^^'^ . We 
proceed as in the proof of Lemma [5751 The bound on the third term on the left-hand side of (|6.1ip 
implies the existence of a subsequence (not indicated) and an element g £ L?{Q, x D), such that 



V q\J C,^'^ ^ g weakly in L2(f7 X D). (6.19) 
Proceeding as in (j3.62p - p.63l) in the proof of Lemma [3.31 with ip" s, ip" and (5^0+ replaced 



by y C^'^, y Sind A ~> oo, respectively; we obtain the weak convergence result: 



and similarly for the a; gradient 



Mi\/q\/(^-^ ^ M^Vq^JC weakly in 2.2(^7 X D), (6.20a) 



M^V^.yC^^^ ^ M^V^yC weakly in L2 (17 X £)), (6.20b) 

as A — > oo. Then inequality (|6.1ip . (|6.20a[ b) and the weak lower-semicontinuity of the L^{n x D) 
norm imply that 

AAt [ M\V^\f^\^dqdx+ AAt [ M \V g\/f^\^ dq dx < [ MT{ipQ) dq dx. (6.21) 
JnxD ~ ~ ~ JnxD ~ ~ ~ JnxD ~ ~ 

After these preparations, we are now ready to state the central result of this subsection. Before 
we do so, a comment is in order. Strictly speaking, we should have written ip'^^ instead of ip'^ in 
our definition (16.151) . as "0*^ depends on the choice of At. For notational simplicity, we prefer the 
more compact notation, V'", with the dependence of V''^ on At implicitly understood; we shall only 
write V'At' when it is necessary to emphasize the dependence of At. Of course, "00 is independent 
of At. 

Next we shall show that, with our definition of ■0'^, the properties listed under (I3.22p hold, 
together with additional properties that we extract from (|6.15p , which were not stated in (|3.22p as 
they will only be required later, in our passage to the limit with L — > oo (and thereby At — > 0+, 
according to At = o(L~^),) in the next subsection. 
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Lemma 6.2 The function ■i/;^ = -0^^ defined by (j6.15l) has the following properties: 



MF{i>°) dq Ax < 



nxD 



MT{ipa)dq dx; 



@ 4 At 



A/J"(-0o)d(7 dx; 



nxD 



O limAt-)-o+ "0" = V'O; weakly in L\f(Q, x _D); 
© limAt->o+ = "00, weakly in L\.j(n x £)). 

Proof 

O This property is an immediate consequence of (|6.16p and the definition p.2ip of Zi. 
© This property was estabhshed in (|6.17p above. 

© The inequahty fohows by using (j6.15p in the left-hand side of (j6.2ip . 

O We begin by noting that an argument, completely analogous to (but simpler than) the one 
in Section l4.2.1l that resulted in (|4.45p . applied to (|6.6p now, yields 



M S dq dx 

At 1 ~ 









/ M 




JnxD 



dq dx 



2 



iYx<?iIl-(d) + iiy?^ 

for all (p e H'^in; L°°{D)) n L^{n; W^^°°{D)). On noting (lOTI) we deduce that 
Af (^A4_^A,0)^^g dx 



dx 



nxD 



(IxD 



MJ"(^o)dg dx 



L--{D) 



|V,(?| 



dx 



(6.22) 



for all ip e H^{n; L°°{D))nL^{n; W^'°°{D)). As the right-hand side of (O^ is independent 
of A, we can pass to the limit A — > cx) on both sides of (|6.22[) . using the strong convergence 
of C^'^ to in L\^{Vl X £)) as A ^ oo (see (|04l) and the definition of ([OS]) ) together 
with the strong convergence of — ^^(iPq) to 0o in -^m(^ x D), as A — oo, with At kept 
fixed. We deduce that 



M (-0° - -00) "^dq dx 



nxD 



< {My 



nxD 



MT{ijjo)dq dx 



2 

L~(D) 



2 



dx 



(6.23) 



for all^e H^{n:L°"{D))nL^{n;W^'°°{D)) and therefore in particular for all ^ e H^iflxD) 
with s > 1 + ^{K + l)d. 

As the last two factors on the right-hand side of (j6.23p are independent of At, we can pass to 
the limit At — > 0+ on both sides of (|6.23l) to deduce that 0" = converges to 0o weakly 
in M-^H^i^n X D)' for s > 1 + i(X + l)d, as At^O+. 

Noting (|6.17p and the fact that T{r)/r — oo as r — >■ oo , we deduce from de le Vallee Poussin's 
theorem that the family {0At}At>o is uniformly integrable in L\^{n x D). Hence, by the 
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Dunford-Pettis theorem, the family {'0^j}At>o is weakly relatively compact in L\j{il x D). 
Consequently, one can extract a subsequence {V'a* Ife^i ^^^^ converges weakly in L\f{il. x D); 
however the uniqueness of the weak limit together with the weak convergence of the (entire) 
sequence -0° = V'a* ^ ipo in M^^H''{^ x D)' , s > 1 + \{K + as Ai 0+ then implies 
that the (entire) sequence 

ifp = converges to V'o weakly in L\i{Q, x Z?), as — > 0+, 

on noting that L\i{Vl x D) is continuously embedded into M^^H^iVt x £>)' for, again, 
s > 1 + ^{K + \)d (cf. the discussion following Theorem 15. ip . 



© It follows from G Zi and pITSl that 



< 



M Ldqdx < 



nxD 



M (i^{i}°)dqAx < I M^j^ dqdx <\Vl\. 



(6.24) 



SlxD 



On noting that J- is non-negative and monotonically increasing on [l,cx)), and that J'{s) G 
[0, 1] for s G [0, 1], we deduce that 



/ M - L]+)Aqdx 
JnxD ~ ~ 

M T{ [ip° - L]+) dq dx + 



< 



M dq dx 



flxD 



MT{i}°)dq dx < C. 



MT{[ij" - L] + )dq dx 



nxD 



Applying the logarithmic Young's inequality ()4.33p . we have that 

(logi) gj" - L]+ < :F{[i/j" - L] + ) + L. 
The bounds (p?^ . (|05l) and (p?^ then imply 



(6.25) 



(6.26) 



nxD 



M[tl;° - L]+dqdx ^ M[il;° - L]+dqdx 



< 



< 



1 

logL 

C 
logi' 



MTiiiJj" - L] + )dqdx- 



L 

JthO 



M L dq dx 



ipo>L 



(6.27) 



Hence for any (p E L°° {fix D) , we have from (|6.27p . on recalling the relationship At = o{L ^), 
that = V'At satisfies 



lim 

At^O+ 



A/(V'° -/3^(V'°))<?d(j dx 



nxD 



= lim 

At-!.0+ 

< I lim 

At^04 



= 0. 



M[i;° - L]+[pdq dx 
nxD ~ ~ 

M[ip" ~L]+dqdx^ ||<?||l~(OxD) 

(6.28) 



+ JnxD 



Therefore, similarly to (|6.23l) . we have that the sequence {V'a* ~ /3^(V''At)}Ai>o converges to 
zero weakly in M-^H^Q. x D)' for s > ^{K + l)d, as At 0+. 

Noting (|6.25p and the fact that F{r)/r — > oo as r — > oo, we deduce from de le Vallee Poussin's 
theorem that the family {V'a* ~ /3^(V''At)}At>o = {[V-'a* ~ -^]+}At>o is uniformly integrable 



45 



m Lj^in X D). Hence, we can proceed as for the sequence {i'%f}At>Q in the proof of O to 
show that the (entire) sequence 

- I3^iip°) = ip%f - P^{ii\t) converges to weakly in L\j{n x D), as Ai 0+, 

on noting that L\.[{Q. x D) is continuously embedded into M^^H'^{Vl x D)' for s > ^{K +l)d 
(cf. the discussion following Theorem IS.ip . Hence, we have proved the desired result. 

That completes the proof of the lemma. □ 

With the information contained under item 3 in Lemma f6.2[ we can now return to the inequality 
()6.2|) . and supplement it with additional bounds, in the sixth and seventh term on the left-hand 
side. The first additional bound can be seen as the analogue of ()4.38|) : 



4/ / M 

10 JnxD 



4At 



M 



dq dx dt 



nxD 

T-At 



|V.V/3^(V'°)l' + |V,V/3^(^o)P 



dq dx 







M 



nxD 



< iAt / M 
JnxD 

Jo JnxD 



dq dx 



M 



T 



"At,+ |2 



dq dx dt 



dq dx dt 



< / MT{ipo)dqdx+ 4 / / M 
JnxD ~ ~ Jo JnxD 



/:At,+|2 



dq dx dt 



(6.29) 



where in the transition to the last line we used p.22p and the bounds on the sixth and seventh 
term in (j6.2p : here and henceforth signifies a generic positive constant, independent of L and 
At. On combining (16.29^ with our previous bounds on the sixth and seventh term in (|6.2p . we 
deduce that 



4/ / M 

Jo JnxD 



^At,±|2 



dq dx dt < Ci,. 



(6.30) 



It remains to derive an analogous bound on V'^)^. To this end, let n G {1, • . • ,N} and consider 
t G {tn-i,tn)', we recall that 



t - t 



At 



(6.31) 



For ease of exposition we shall write 

t — tn-l 



7+ 



At 



and 7_ := 



tn t 

At 



in the argument that follows, noting that 7+ + 7_ = 1 and both 74- and 7_ are positive. The 



functions t G (tri-i,^n) I— >■ i^^L^i'i 'lO ^I'S constant in time and ^Z'"]^'*^^'' (a;, t) > on x 
{tn-i,tn), n E {1, . . . ,N}. For any a G (0, 1) we have that 



-At,(±), 
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< 2 



< 2 



At,+ 



Hence, on bounding 7^ by 1, we deduce that 



< 2- 



7At.- , ■ 
T. + a 



for all {x,q,t) enxDx (t„-i, t„), n = 1, . . . , iV, and all a G (0, 1). On multiplying (|6:32|) by M, 
integrating over 57 x £) x (i„_i, summing over n = 1, . . . , A^, and passing to the limit a — s- 0+ 
using the monotone convergence theorem, we deduce that 



JnxD 



MlV^Jip^ll dgdxdt < 2 



/O JdxD 



M\Y,JiPfl:+\^dq dxdt 



JnxD 



MlYxJ^ti'l^dq dxdt 



Analogously 
4 



JnxD 



M\Yq^/i'fX\ dqdxdt < 2 



/o JnxD 



JnxD 



M\Vg\/Tpf'£-\'^dqdxdt 



Summing the last two inequalities and recalling (|6.30p . we then deduce that 

'lYxA/v^^iP + iy^Jv^^iP dqdxdt<c,, 



M 



JnxD 



(6.33) 



where, again, C* denotes a generic positive constant independent of L and At. 

On noting that (1 + \q\Y^ < + \q\^^) < 2^^'^K'^-\1 + Y.f=i \q^?'^), it follows from 

([Q]) and (p:^ that for a.e. t € [0, T] ~ 



SlxD 



M{\ + \q\f''^f'j^{t)dqdx<C.. 



(6.34) 



Finally, on combining (|6.30l) . (|6.33p and (|6.34l) with (|6.2I) we arrive at the following bound, 
which represents the starting point for the convergence analysis that will be developed in the next 
subsection. 

With a > l + \d and s > 1 + \{K + \)d, we have that: 



ess.suptg[o^T]||M^[^'^^(t)||2 + 



At.- 



ds+ r \\V,uft^\s)fds 
Jo « ~ ' 



-ess.suptg[oy] 



nxD 



M F{i}f'^'^\t))dq dx + 



1 



nxD 



■ ip^L )^dgdxds 



■ess.supjgfo-T] 

fT 



Mil + \q\r'>i^t[^'^\t)dqdx 



nxD 



'0 JnxD 



M\V^\hpf'J;-^^f dq dxds 



JnxD 



rAt(,±)|2 



dq dx ds 



+ 



'0 



dufX 
dt 



dt 



M- 



dt 



dt < a. 



(6.35) 
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6.2 Passage to the limit L — )• oo 

We are now ready to pass to the limit and prove the central result of the paper. In what follows, 
(■j ■)_ff=(f2) denotes the duality pairing between H'^^n)' and H''{n,) relative to the pivot space 
i^(r2) with inner product (•,•); similarly, •)j;fs(QxD) denotes the duality pairing between 
M^^H'^{n X D)' and x D) relative to the pivot space L\.j{Q x D) with inner product 

(01, 02)m := / M 01 02 dg dx; 
JnxD 

and (•, ■)v^ denotes the duality pairing between the spaces V'^ and Va relative to the pivot space 

Theorem 6.1 Suppose that the assumptions p.Sp and the condition (|6.5p . relating At to L, hold. 
Then, there exists a subsequence o/{(m^^, '0^^)}l>i (not indicated) with At — o{L~^), and a pair 
of functions {ug^^il);,) such that 

Ue e L°°(0, T; L^{n)) n L^{0, T; V) n H\0, T;V'„), ct > 1 + id, 

and 

e Li(0,r;Lii^|g|)2.M(^ X D))nH\0,T;M-^H''{n X D)'), s > 1 + ^{K + l)d, 
with il'e>0 a.e. on Q. x D x [0,T], 

M{q)i)^{x,q,t)dq = l for a.e. {x,t) e Vl x [{),T], (6.36) 

and hence ipe G L°°{Q,T\L\j(Vt x D)); and finite relative entropy and Fisher information, with 

Tii^e) e L°°{0,T;Lli{nx D)) and e L^{Q,T;HIj{VL x D)), (6.37) 

whereby ips G -^""(O; ^(i+|g|)2f m(^ ^ ^))/ such that, as L —> oo fand thereby At — ^ 

^'f L ^"^^ ^ weafc* m L°°(0, T; L^(f7)), (6.38a) 

L ^"^^ ^ weakly in L^{0, T; F), (6.38b) 

uf.L'^^ ^ strong//?/ in L^{0, T; L''(r2)), (6.38c) 

^ ^ weakly m L^{Q, T; \/;,), (6.38d) 
where r G [1, oo) if d = 2 and r G [1, 6) if d = 3; and 



D 



M*V:r\/i/'^L'^^ ->M5 V:r\/^e weakly in L\Q,T;L\Vtx D)), (6.39a) 



M^\Iq^ij'^2 ' ^ M-^VqSji^e weakly m L\0,T; {fix D)), (6.39b) 

M^^^M^ weakly m L^(0,T;H%nx D)'), (6.39c) 

V^e'^L '^^ ^ V^e sirori3/2/znLf(0,T;Lii+lqlp,^,(r!xD)), (6.39d) 
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for all p ^ [1,00); and, 



K 



K 



■Y,C^{Mi,f^^■'^'^) ^ V, ■Y,C^{M4,e) weakly m L^{0,T;V',). 



The pair of functions {u^,ipg) is a global weak solution to problem (Pe), in the sense that 



.39e) 



\~dt 



■,w ) dt 



T 

Jn 



dx dt 



f {f,w)v~kY, I C^{M^^,):^,wdx 



dt 



(6.40) 



and 



dt 



„T « K K 



M 



JnxD 



Vj;(^ dq dx dt 



+ TT f I ^^Y. IZ^^J ^9.^^ • ^9'^ dqdxdt 
^■^Jq JnxD i=i j=i ~ ~ ~ ~ 

r"^ f ^ 

-/ / My][a{u^)q.,]^j^ ■ Wq^ip dqdxdt = G L^(0, T; x D)). (6.41) 

Jo JnxD j^]^ » ~ ^ ~ ^ ~ 

The initial conditions We(-, 0) = yo(') md tpe{', 0) — V'o(', •) are satisfied in the sense of weakly 
continuous functions, in the function spaces Cw{[0,T]; H) and Ciui[0,T]; L\f(nx D)), respectively. 
The weak solution {ug,il)g) satisfies the following energy inequality for t & [0:^]-' 



\\u,{t)f + vj \\V^u,{s)fds+ 2k f MT{i^,{t))dqdx 

JnxD ~ ~ 

nt 







+ 8ke 



JnxD 



A/|Va;V V-el dg " 



2aofc 



JnxD 



M |V,Y?Aer dqdxds 



<lkof + 



1 



2k 



MJ-(^o)dgdx< [B(wo,/,V'o)]' 



(6.42) 



nxD 



with J^{s) — s(logs — 1) + 1, s > 0, and [B(uo, /, V'o)]^ as defined in (|4.28bp . 



Proof Since the proof is long, we have broken it up into a number of steps. 

Step 1. On noting the weak* compactness of bounded balls in the Banach space X°°(0, T; 
and recalling the bound on the first term on the left-hand side of (|6.35l) , upon three successive ex- 
tractions of subsequences we deduce the existence of an unbounded index set C C (1, 00) such that 
each of the three sequences {u^^^'^^}l&c converges to its respective weak* limit in i°°(0, T; L?{VL)) 
as i ^ 00 with L £ C. Thanks to (|4.1al b'). 



Wv^fTi^) - yf:r(^)lPd.s < iC.Ai, (6.43) 



where the last inequality is a consequence of the second bound in (|6.35p . On passing to the 
limit L — )■ 00 with L e £ and using (16.51) we thus deduce that the weak* limits of the se- 
quences {u^^l^'^^^Lec coincide. We label this common limit by Ue] by construction then, e 
L°°(0,r;L2(f7)). Thus we have shown (|6.38ap . 
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Upon further successive extraction of subsequences from {u^ ^ }lg£i ^-nd noting the bounds 
on the third and eighth term on the left-hand side of (|6.35l) the hmits (j6.38bl d') follow directly 
from the weak compactness of bounded balls in the Hilbert spaces L^{0,T;V) and L^{0,T;V'^) 
and (j6.38al) thanks to the uniqueness of limits of sequences in the weak topology of L^{0,T; V) 
and L^{0, T; V'^), respectively. 

By the Aubin-Lions-Simon compactness theorem (cf. (|3.14p '). we then deduce (|6.38cp in the 
case of on noting the compact embedding of V into L^{^) H H, with the values of r as in 
the statement of the theorem. In particular, with r ~ 2, — > u^, strongly in L'^{0,T; L^{i})). 
Then, by the bound on the left- most term in (|6.43l) . we deduce that uf^£^ also converges to u^, 
strongly in L^(0, T; as L — > oo (and thereby At — )• 0+). Further, by the bound on the 

middle term in (j6.43l) we have that the same is true of u^'j^ ■ Thus we have shown that the three 
sequences u^^'^'' all converge to u^, strongly and i^(0, T; ^^{Vt)). Since the sequences ^'''^^ are 
bounded in L^(0, T; H^{Vt)) (cf. the bound on the third term in (16.351) *) and strongly convergent in 
i^(0, T; we deduce from p.4p that (|6.38cp holds, with the values of r as in the statement 

of the theorem. Thus we have proved (|6.38al -d). 

Step 2. Dubinskii's theorem, with ^q, Ai and M as in the discussion following the statement 
of Theorem 15. 1[ and selecting p = 1 and pi — 2, imply that 

: [Q,T]^M : Mmn.T-M) + 

^ L\0,T;Ao) = ii(0,T;Li,^|^|).,j,,(r! x D)). 

Using this compact embedding, together with the bounds on the sixth, the seventh, the eighth and 
the last term on the left-hand side of (|6.35p . in conjunction with ()6.3p . we deduce (upon extraction 
of a subsequence) strong convergence of {V'^|^}l>i in L^{0,T; L^^^^q^-^2ii x D)) to an element 

i^e e 2.1(0, T;Li^^|^|)2^^,^(f7 X £>)), as L ^ oo. 

Thanks to the bound on the fifth term in (j6.35p . by the Cauchy-Schwarz inequality and an 
argument similar to the one in (I6.43p . we have 

n M (1 + \q\Y^ - V^f I dq dx At 

Jo JnxD ~ ~ ~ 

<IM( f M(l + \q\rdq)n M(V^fi'+-V^fi'-)2dgdxdt 
■J \Jd ~ ~/Ja JnxD ' ' ~ ~ 

< \n\ AtL (^j^ M (1 + \q\f'> dq^ . (6.44) 

On noting from (|1.4cl) that the integral over D is finite and recalling (16.51) , and using the triangle in- 
equality in the 2.1(0, T; L^^^j^^^^yf, x D)) norm, together with (|6.44p and the strong convergence 

of {V-'f^i}L>i to tpe in 2.1(0, T; iJ^^|^|^2tfA,/(f^ x D)), we deduce, as L — ^ oo, strong convergence of 

{V'^L'^}i>i ill ^HO,T;iJi+lgl)2^M(^><-0)) to tl^e same element V^e € 2.1(0, T; 2.i^^|^|,,^^^(l]x2:>)). 

This completes the proof of (|6.39dp for p = 1. 

From the sixth bound in (|4.18p we have that the sequences {"0^^ '^''}l>i are bounded in 

L°°(0, T; 2.i-i^^l^lj2fl^,j(f^ X 23)). By Lemma 15.11 the strong convergence of these to V'e in 2.1(0, T; 

2.Jj^lg|-,2Sj,/(f^x 2?)), shown above, then implies strong convergence in 2.^(0, T; 2.i^^|^|j2tf j^,/(f^ x 2))) 

to the same limit for all values of p G [1, oo). That now completes the proof of (I6.39dp . 

Strong convergence in 2.^^(0, T; 2. J-|^^|^|j2tfj(,f(i^ x 2))), p > 1, implies strong convergence in 

2.^(0, T; L\j{VL X D)), p > 1, and convergence almost everywhere on f2 x 2) x [0, T] of a subsequence. 



dip 



< oo 



L^{0,T-Ai) 
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Hence it follows from (|6.3p that -ip^ ^ on fl x D x [0,T]. Furthermore, by Fubini's theorem, 
strong convergence of {ipf*]^''^''}L>i to ipi. in L^{0,T; L\j{n x D)) implies that 

/ M{q)\^}f^j;}^\x,q,t)-^e{x,q,t)\dq^Q as L ^ oo 

J D ~ ' ~ ~ ~ 

for a.e. {x,t) £ fix [0,r]. Hence we have that 

/ M{q)i)f]^'^\x,q,t)dq / M{q)%jje{x,q,t) dq as L oo, 
Jd ~ ~ ~ jd ~ ~ ~ 

for a.e. {x,t) <E fl x [0,T], and then (16.41) implies that 

/ M{q)i^eix,q,t)dq <1 for a.e. (x, i) G O x [0, T]. (6.45) 

We will show later that the inequality here can in fact be sharpened to an equality. 

As the sequences {V-'^1'''^''}l>i converge to V'e strongly in L^{G, T; L\j{Vl x D)), it follows that 

(upon extraction of suitable subsequences) they converge to V'e a.e. onilx D x [0,T]. That then, 

in turn, implies that the sequences {J^{^f^£^^^)}L>i converge to T{ipe) a.e. on Q. x D x [0,T]; in 

particular, for a.e. t £ [0,T], the sequences {•7-"(^^^'^''^''(-, ■,i))}L>i converge to J^(V'£(-, ■, t)) a.e. 
onftx D. Since is nonnegative, Fatou's lemma then implies that, for a.e. t e [0,T], 

/ M{q)J^{tpe{x,q,t))dxdq<\\mmfL^oo M {q) {tpf'^'^ {x , q , t)) dx dq 

JnxD ~ ~ ~ ~ ~ JnxD ~ ' ~ ~ ~ ~ 

< ^[B(uo,/,^o)]', (6.46) 

where the second inequality in (|6.46p stems from the bound on the fourth term on the left-hand 
side of (|4.28b|) . As the expression on the left-hand side of (j6.46p is nonnegative, we deduce that 
T{i^e) belongs to L°°{Q,T]L\j{Q. x D)), as asserted in the statement of the theorem. 

We observe in passing that since | — y/c2 \ < ■\/|ci — C2I for any two nonnegative real 
numbers ci and C2, (j6.39dl) directly implies that, as L 00 (and At — !■ 0+), 

y^J^^y^ strongly in LP(0,r;L2^(l]xZ?)) Vpe[l,oo), (6.47) 
and therefore, as L ^ cx), 

M5 strongly in iP(0,r;i2(f} X D)) Wpe[l,(X)). (6.48) 

By proceeding in exactly the same way as in the previous subsection, between equations (|6.18p 
and (|6.20al b) with C'^'^j = "0° and A replaced by ips and L, respectively, but now using 

the seventh and the eighth bounds in (|6.35p . we deduce that (|6.39a[ b) hold. 

The convergence result (|6.39cl) follows from the bound on the last term on the left-hand side 
of (|6.35|) and the weak compactness of bounded balls in the Hilbert space L^{Q,T; H'^{il x D)), 
s> l + ^{K + l) d 

The proof of (|6.39ep is considerably more complicated, and will be given below. 

After all these technical preparations we are now ready to return to (|4.3a[ b') and pass to the 
limit L 00 (and thereby also At 0+); we shall also prove (|6.39eP and will also pass to the 
limit on the initial conditions for (I4.3al b) . Since there are quite a few terms to deal with, we shall 
discuss them one at a time, starting with equation (j4.3b[) . and followed by equation (|4.3ap . 

Step 3. We begin by passing to the limit L -> 00 (and At — > 0+) on equation (j4.3bp . In 
what follows, we shah take (p G C{[G,T]]C°°(rt;C^{D))). Note that C°°(II; (7^(13)) is dense in 
L^{n;H^{D)) n H'{n;L^{D)) = H^iil x D), and so C([0, T]; C°°(n; C^(Z)))) is a dense linear 
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subspace of L^(0, T; iJ^(rj x D)) for any s > 0. As each of the terms in (j4.3bp has been shown to 
be a continuous linear functional with respect to cp on L^(0, T; 7J*(f2 x D)) for s > 1 + ^{K + l)d, 
the replacement of L^{0, T; H'{n x D)) by its dense linear subspace C([0, T]; C°°(n; C^{D)) for 
the purposes of the argument below is fully justified. 

Step 3.1. Passing to the limit on the first term in ()4.3bp is easy: using (|6.39cp we immediately 
have that 



dt ~ \ dt \ dt 

as L -> oo (and At 0+), for all If e C{[0,T];C°°{TI,C^{D))), as required. That completes 
Step 3.1. 

Step 3.2. The second term will be dealt with by decomposing it into two further terms, the 
first of which tends to 0, while the second converges to the expected limiting value. We proceed 
as follows: 



/ / MVjf'^+-V.,^dqdxdt 
Jo JnxD ~ ' ~ ^ ~ 

= 2e [ I M J V, J ^^['+ • V,^ dg dx dt 

Jo JnxD V ' ^ V ^ ^ ^ 

= 2e [ [ M (Ji^f'£+ - vJ^Tl'^ ■ V,(? dq dxdt 

Jo JnxD \ J ~ 1 ' ^ ^ ^ 



+ 2e [ I M \[^e^ x\l ^^f^L^^ ■ Vx"^ dqdxdt 
Jo JnxD ~ V ' ^ ^ ^ 

=: V1+V2. 

We shall show that Vi converges to and that V2 converges to the expected limit. 



T 



|Vi|<2£ / / ( / M|J^;^i'+- V^el^dg 
Jn \Jd * 



'y^M|V,^/J5+pdgy ||V,^|Uo.(o) dxdt 

<2e r ( f MlJ^-^.l^dqdxY 
Jo \JnxD ~ ~J 

V M\V ^J^^f dqdX \\V ML^{nxD)dt 
\JnxD ~ V ^ ^ 



2e / ll\A^^5^-\A^IU^a^xi3) 



X V2:y>flJ^||L2(J2x_D) l|Va;^||L~(nx_D) dt 

iyM^5+-VM^irL.(^,xi5)di) ( riiv.^iiS(oxD)d^ 



were r £ (2, 00). Using the bound on the sixth term in (|6.2p together with the Sobolev embedding 
theorem, we then have (with C* now denoting a possibly different constant than in (|6.2p . but one 
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that is still independent of L and At) that 



|Vi| < 2CH\\jM^f^l+ - ^Mi,, ||l.(o,t;L^(oxd)) ||V,(?|| 



L'^(0,T;L°=(axD)) 



where we also used the elementary inequality \y/ci ~ \fca\ < \/|ci — C2I with ci,C2 G K>o. The 
norm of the difference in the last displayed line is known to converge to as L — ?> 00 (and Ai -> 0+), 
by (|6.39dp . This then implies that the term Vi converges to as i — > 00 (and At — > 0+). 
Concerning the term V2, we have that 



V2 = 2£ 



JnxD 



M t Vx^ipf^l:'^ ■ V M V^(p dq dx dt. 



Once we have verified that y M -ip^Yx'P belongs to L^(0, T; _L^(f] x D)), the weak convergence 
result (|6.39al) will imply that 



V2 ^ 2e 



M^V^VV'e- y M^p^V^ip dqdxdt = e 



M Vx^s ■ x'P d(j da; dt 



JnxD 



JnxD 

as L — > CXI (and At — > 0+), and we will have completed Step 3.2. Let us therefore show that 

\J M ip^ Yx^ belongs to L^{0, T; L^{^ x D))\ the justification is quite straightforward: using (|6.45l) 
we have that 



JnxD 



M |Vj,(^pdg dxdt 



JnxD 

T 



< 



< 



Jn 

T 



I'^2:<?IIl~(D) ( / Mip^dq] dxdt 



Va;(^||?oo,r)^ dx dt 



Jn 



L^iO,T-L^{n-L°^(Dy)} 



< 00. 



That now completes Step 3.2. 

Step 3.3. The third term in (j4.3bp is dealt with as follows: 



JnxD 

T 



AI uf*£ i^f\P^ ■ xW dqdxdt 



M Ue ipe • Va;(^ dq dx dt 



JnxD 



JnxD 



M {ue — uf*£ ) 'ipf]]^ ■ Va;^ dq dx dt 



+ I I Mu,{^,- ■ V,^ dq dx dt. 



JnxD 



We label the last two terms V3 and V4 and we show that each of them converges to as i — > 
(and At — >■ 0+). We start with term V3; below, we apply Holder's inequality with r £ (l,oo) in 
the case of d = 2 and r G (1, 6) when d = 3: 



|V3| = 



< 



Jn 

T 



) dg 


da; dt 







Jn 



We - U 



At,- 

'6,L 



AfV^fi'+dg 



\'^x'f\\L-°{D) dx dt 



53 



< 



/ / lUe-ufi I ||Va;^||Loo(D) dxdt 

Jo Jn ~ ~ 



< 



T 

II ^ 

L 

At,- 



11^^ - ^e.L \\LHO,T;Lrin)) II Y-"?!! i2(o,T;L^ (0;L- (D))) ' 



where in the transition from the second hnc to the third line we made use of (j4.31bp . Thanks to 
(|6.38cp the first factor in the last line converges to 0, and hence V3 converges to also, as L — >■ 00 
(and At ^ 0+). 

For V4, we have, by using Fubini's theorem, the factorization 

M I 



- i'tr) = (\/^- V^) (\/^+ V^) ' (6.49) 

together with the Cauchy-Schwarz inequality, (j4.31bp . ()6.45p and the elementary inequality 
I y/ci - ^/cj I < a/|ci - C2I with ci , C2 € M>o , that 

IV4I < f f M \u,\ \$, - |V,^| dq dxdt 

Jo JnxD ~ ' ~ ^ ~ 

rp 1 

<2^ j^\u,\ (^^A/|V^s- v/^S^l'dgV ||V:,^|U~(Z5)djdt 

rp 1 1 

<2^ (^Iju.fdxY ^A/IV^-^^g+l^dgd^y \\Yx^\\L^in>,D)dt 

<2/ \\ue\\L^n) ( Af - dgdx ) W'^x'fWL^inxD) <it 

Jo ~ \JnxD ' ~ -/ ~ 

< 2||Ue||Loo(o,T;L2(o))||4 -'?i^L^II^(0,T;Lj,,(f^xI)))llY^^II^'"(0,T;L-(axD))• 
By (I6.38ap the first factor in the last line is finite while, according to (j6.39dp (with p — 1), the 
middle factor converges to as L 00 (and At 0+). This proves that V4 converges to as 
L ^ 00 (and At — ^ 0+), also. That completes Step 3.3. 
Step 3.4. Thanks to (|6.39bp . 



^q^ji^f^i'^^ A/5 Vg-y/i^e weakly in L^{0,T; L'^{n x D)). 

This, in turn, implies that, componentwise. 



^<i,\J'^f!L^^ ^q,\l^e weakly in L'^{0,T; L^{n x D)), 

for each j = 1, . . . , K, whereby also, 

K K I — 

Af 5 ^ A.jVq^ V^S^ M^- J2 V weakly in 1^(0, T; L^{fl x i?)). 

J = l ~ 3 = 1 

That places us in a very similar position as in the case of Step 3.2, and we can argue in an identical 
manner as there to show that 



Jo JnxD 
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K K 



Jo JnxD .^-^ ~ ~ ~ ~ 



as L ^ oo and At — > 0+, for all (p e L^-^ (0, T; H^^'°°(r2 x D)), where r £ (2, oo), and in particular 
for all if e C{[0,T];C°°{Ti; C^{D))). That completes Step 3.4. 

Step 3.5. The final term in (|4.3bl) . the drag term, is the one in the equation that is the most 
difficult to deal with. We shall break it up into four subterms, three of which will be shown to 
converge to in the limit of L — >■ cxd (and At — >■ 0+), leaving the fourth term as the (expected) 
limiting value: 



/ / MJ2 ['^("f L^) 1^] (V^f L + ) • V,, ^ dq dx dt 

Jo JnxD j^]^ « ~ ^ ' ~ ^ ~ 

JO JnxD j^]^ L ^~ ~ 

Jo JnxD ,_i 



,At,+ 
' .L 



/3^(V^f['+) • Vg^fidqdxdt 



/ / ^'^E (Y-!f^i^) 9' (/3''(^s)-V^e) • Vg,(?dqdjd^ 



dt 



1=1 



Jo JnxD .^^^ ~ ^ ~ 



ipe • Vg.^ dg dx di 



• Vq.i^ dg da; dt. 



(6.50) 



We label the first three terms on the right-hand side by V5, Vg, V7, respectively, and we proceed 
to bound each of them. We shall show that each of the three terms converges to 0, leaving the 
fourth term as the limit of the left- most expression in the chain, as L — >■ 00 (and At — >■ 0+). 

We begin by bounding the term V5, noting that cp is fixed with compact support in D, 
is Lipschitz continuous, using the factorization (|6.49p together with (j4.31bl) and (|6.45p . and then 
proceeding as in the case of term V4 in Step 3.3: 



|V5|< 



M \q\ |V,u^i'+| - |V,^| dq dx dt 



JnxD 

T 



< 



'E,L 



M l^t.L^-y^el dq ] dx 



n ~ ~ \JD 

< 2 ||Va;W^^'^||L2(o^T;L2(o)) 



IqI Vg^llioo(oxD) dt 



TAt,+ 



A\\h{0,T:L\,inxD)) II I9I Y9'?IU°°((0,T)xOxD))- 



On noting the bound on the third term on the left-hand side of (|6.35p and the convergence result 
(|6.39d[) that was proved in Step 2, we deduce that term V5 converges to as i -> 00 (and 
At 0+). 

We move on to term Vg, using arguments similar to those used for terms V5: 



IVel < 



M \q\ \V^uf 'j-+\ |/3^(V'e) - A\ |V,^| dq dx dt 



JnxD 

T 



< 



I / M|/3^(^,)-^,| dq dx 



\q\yq(p\\L-^(nxD) dt 



< 2 llW^u 



At,+ | 



2(0,r;L2(O)) ll/3^(V'e) - V'e|lii(o_T;L]„(OxL>)) H I? I Y^^ll L~((0,T) xfJxD)) • 
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Note that < t/i^ — I3^{tpe) < i^e and that -ip^ — f3^{'tp^) converges to almost everywhere on 
X D X (0, T) as L ^ oo. Note further that, thanks to ()6.39dp with p = 1, ■0^ e L^{0, T; Llj{n x 
D)). Thus, Lebesgue's dominated convergence theorem imphes that, the middle factor in the last 
displayed line converges to as L — > oo. Hence, recalling the bound on the third term on the 
left-hand side of (|6.35|) . we thus deduce that Vg converges to as L — > oo (and At — > 0+). 
Finally, we consider term Vy: 



V7 



/7 '''^[(^ 

Jo JnxD 



At,+ 



tp^ ■ Vq-(p dq dx dt. 



We observe that, before starting to bound V7, we should perform an integration by parts in order 
to transfer the x gradient from the difference 



At.H 

e,L 



onto the other factors under the integral sign, as we only have weak, but not strong, convergence 
of Yxuf*!^^ — YxUe to 0, (cf. (|6.38b|) ) whereas the difference uf]^^ — converges to strongly, 
by virtue of (|6.38cl) . 

We note is this respect that the function a; e f2 ipe{x,q,t) has a well-defined trace on dfl 
for a.e. (g,t) £ D x (0,T), since, thanks to (|6.39ap . 



''tpe{-,q,t) e and therefore 

for a.e. e D x (0,r), implying that 



an 



for a.e. {q, t) £ D x (0, T), with 2p e [1, 00), when d — 2 and 2p £ [1, 4] when d = 3, whereby 

M-,q,t) eLP{dn) 
an 

for a.e. {q,t) G Dx (0, T), with p € [1, 00), when d = 2 and p G [1, 2] when d = 3. As the functions 

Mg, uf'^£~^ have zero trace on d^, the boundary integral that arises in the course of integration by 
parts then vanishes. First, we write 



V7 



M 



JnxD 



K d 

E E , 

i—l m,n— 1 



d_ 

, dXrr 



V^e (Vg,^)„ dq dxdt. 



Here, {u^*j^^)n and {ue)n denote the nth among the d components of the vectors u^^i^ and Ue, 
1 < n < d, respectively, and (yg.^)„ denotes the nth among the d components of the vector V q^(p, 
1 < n < d, for each i e {1, . . . , K}. Similarly, {qi)m denotes the mth component, 1 < m < d, of 
the d-component vector qi for i e {1, . . . , K}. Now, on integrating by parts, and cancelling the 
boundary integral terms, with the justification given above, we have that 



V7 = 



T !• K d r- 

^^E E (i^fX')r.-{Ue)n) (?.),: 



d 



T p K 

M 



JnxD 



i—l 7n,Ti— 1 '- 

K d r 



9a;„ 
9x„ 



(^V^j (V5i(^)„^ dg da:; di 
(Vq;<^)„ dg da; dt 



„T !• K d 

'0 JnxD 



-1 m,n— 1 



t/'e ( Vg, if)n ] dq dx dt 
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= : V7.1+V7,2- 

For term Vy^i, we have that 



|Vr,i|< / / |(<Dn-("e)n \{q^U 

Jo JnxD i^irn,n=l ' ~ ~ ~ 

/ K d 

E E |(! 



dx„ 



dq dx dt 



1 

lo JnxD 



4^) 



{Ue)n \{qi)rn\" 



?n,n— 1 



E E 

i—l m,7i—l 



dXr, 



dq dx dt 



M - ue\ IgllV.T/'el |V,(^| dg dxdt 



JnxD 

T 



< 



M\u 



nxD 



At,- 
e.L 



, dt 



At,+ 



< 2 



I At,+ 
K.L -"e 



I \y dq dxj II \q\yq(p\\L'^{nxD}' 

M\/^|V:r\/^|dg^ da; \\\q\Yq'f\\L=-{nxD)dt 

dt, 



where in the transition to the last hne we used the Cauchy-Schwarz inequahty in conjunction with 
the upper bound (|6.45p . Hence, 



|V7,l| < 2 ||uf^^'+ - Ue||i2(o,T;L2(0)) l|Vj;\/?Ae||i2(o,T;L|j(OxD)) II kl ||i<»(o,T;L~(nxD))- 

Noting ()6.38cp with r = 2 and (|6.39ap we then deduce that Vy^i converges to as L — (and 
Ai-^0+). 

Let us now consider term V7.2. We proceed similarly as in the case of term YtX- 



7,2 



< 



< 



< 



< 



lo JnxD 



K 



i—l m,n— 1 



d 



dXn 



lo Jn 

i-T 



lo Jn 

i-T 



.^■^E E \i^f!^^)n-iUe)n\\iqjUt 

M\q\iie |V^Vq^|dg^ dajdi 
Mijj^dq] II |g| Va;Vq^|lLoo(£,) dxdt 



dq dx dt 



At,+ 
'e,L 



At,+ 



D 



Jn 



At,+ 
'■e.L 



Ue\ II kl Va;Vq^||L=o(£,) dxdt, 



where in the transition to the last line we used (|6.45l) . Hence, 



|V. 



7,2 



< llw 



At,+ 



"ellL2(0,T;L2(o)) II \q\y x"^ q^\\L■^iO,T■.L^{n■.L•^{D)))■ 



lS^otmg (|6.38cp with r = 2, we deduce that Vg_2 converges to as i — > (and At — > 0+). 

Having shown that both ¥7^1 and ¥7^2 converge to as L — >• (and At — ^ 0+), it follows 
that the same is true of V7 = ¥7^1 + ¥7^2- We have already shown that V5 and Vg converge 
to as i — 7> (and At — ^ 0+). Since the sum of the first three terms on the left-hand side of 
(|6.50p converges to 0, it follows that the left-most expression in the chain (|6.50p converges to the 
right-most term, in the limit of L — cx) (and At — >■ 0+). That completes Step 3.5. 
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Having dealt with (I4.3b[) . we now turn our attention to (I4.3a|) . with the aim to pass to the 
hmit with L (and /Si). In Steps 3.6 and 3.7 below we shall choose as our test function 

w e C([0,r];C'g°(17)), with ■ w = on rj, for ah i e [0,T]. 

Clearly, any such w belongs to i^(0,T; V) and is therefore a legitimate choice of test function in 
(|4.3ap . Furthermore, the set of such smooth functions w is dense in L'^{0,T; Va), a > 1 + ^d. As 
each term in (|4.3al) has been shown before to be a continuous linear functional on L^(0,T; Va), 
a > 1 + the replacement of L^{0, T; Va), cr > 1 + ^d, with such smooth test functions for the 
purposes of the argument below is fully justified. 

Step 3.6. The terms on the left-hand side of (|4.3ap are handled routinely, using (|6.35l) and, 
respectively, (|6.38dl) . (j6.38cl) with r — 2 and (j6.38b|) . In particular, the second (nonlinear) term 
on the left-hand side of (|4.3al) is quite simple to deal with on rewriting it as 



and then considering the difference 



(We Ue - (Ki uf*£ , Yx w) dt, 



which is bounded by 




Ue (g) - (g) uf ]^ h^in) di I lly^r ?«|| L~(0,T;L~ (H)) • 



By adding and subtracting (8) uf*£~ inside the first norm sign, using the triangle inequality, 
followed by the Cauchy-Schwarz inequality in each of the resulting terms, and then applying the 
first bound in (|6.35|) . and (|6.38cp with r = 2, we deduce that the above expression converges to 
as L — ?> oo (and At — > 0+). The convergence of the first term on the right-hand side of (|4.3a|) 
to the correct limit, as L — > cx) (and At 0+), is an immediate consequence of (|3.25p . We refer 
the reader, for a similar argument in the case of the incompressible Navier-Stokes equations, to 
Chapter 3, Section 4 of Temam [39^. That completes Step 3.6. 

Step 3.7. The extra-stress tensor appearing on the right-hand side of (|4.3ap is dealt with as 
follows. We have from (|1.4bp . and similarly to (|4.5ip . that 



k y2 CiiMi^^l^^) -.X^^wdxdt-k / V / C.,{M^Pe) -"^xwdxdt 
Jo Jn~ '' a, ~ ~ Jo Jn ~ ~ ~ ~ 



T K 

E 

i=l 



flxD 



M U- q^ (?f (V'f 1'"*" - ^e) -"^xwdq dx dt 



< 2fc||V,«;|Uc 



= (0,T;L~(n)) / E / 

'0 JnxD 



c»2 1 hkt? 



M{l + \q\)^'>\i^^'£+ -$,\dqdxdt. 



JnxD 



At,- 



t/j;, I dq dx dt 



By (j6.39dp with p = 1 the integral in the last line converges to zero as L — ?> oo (and At — >■ 0+). 
The fact that Vs converges to then directly implies (|6.39ep . thanks to the denseness of the 
set of divergence-free functions contained in C([0,T]; C'^{fl)) in the function space L^(0,T; Va), 
a > I + ^d. That completes Step 3.7, and the proof of (|6.39ep . 

Step 3.8. Steps 3.1-3.7 allow us to pass to the limits L — oo and At 0+, with L and 
At linked by the condition At — o{L^^) as L — > oo, to deduce the existence of a pair {ue,^e) 
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satisiying (|6.40l) . (|6.4ip for smooth test functions <f and w, as above. On noting the denseness of 
the set of divergence-free functions contained in C([0, T]\C'^{^)) in L?{Q, T; V^), a > 1 + ^d, and 
the denseness o{C{[0,T]]C'^{Ti;C^^{D))) in L'^{0,T; H^^nx D)), s > l + l{K+l)d, respectively. 
That completes Step 3.8. 

Step 3.9. The weak continuity of and stated in the theorem is shown as follows. First 
we consjxler 

As e H^iO,T;M-'^H'{n x D)'), s > 1 + ^{K + l)d, it follows by the Sobolev embedding 
theorem that ij^^ G C{[0,T]; M-'^ H" {n x D)'). In fact, since G L°° {0,T; L\j{n x D)), it is 
possible to show that 'tpe € C'u){[0,T'];i]^^(^2 x D j), i.e. V'e is weakly continuous as a function from 
[0,r] into i]v/(^ ^ -^); particular, (■i/'e(-, ■,t)-ipo{-, •), i^) -J> as i 0+ for all if G ^mC^ ^ -^)'' 
where (•, •) denotes the duality pairing between L\f{il, x D) and L\j{n x D)'. In order to prove 
the statements in the last sentence we invoke Lemma 1.4 in Chapter 3 of Temam [39] . which we 
quote here. 

Lemma 6.3 Suppose that X and Y are Banach spaces, with X continuously embedded into Y . 
If ip Cz L°°{0,T; X) and is weakly continuous as a function with values in Y , then ip is weakly 
continuous as a function with values in X . 

By taking X = Llj{n x D) and Y = M^^H^in x D)' with, once again, s > 1 + \{K + l)d, the 
continuous embedding of X into Y follows (cf. the discussion following Theorem 15. ll for a proof); 
Lemma [^751 then implies that ip^ G Cu,([0,T]; L\i{^ x D)). On recalling point © of Lemma l6.2l and 
(|6.39c[ d) and following the abstract framework in Temam j39j. Chapter 3, Section 4, we deduce 
that Tp satisfies the initial condition in the sense of weakly continuous functions from [0,r] into 
L\j{il X D). The proof for is similar, by taking X ~ H and Y ~ V'^ and noting that u'^ — > uq 
weakly in H as At —?' 0+. That completes Step 3.9. 

Step 3.10. The energy inequality (j6.42p is a direct consequence of (j6.38al -c') and (|6.39al b.d'). on 
noting the (weak) lower-semicontinuity of the terms on the left-hand side of (I4.28bl) . and (|6.46p . 
That completes Step 3.10. 

Step 3.11. It remains to prove (|06l) . To this end, it follows from (|429|) . (|446l) and on 
recalling that G L^(0,r;/C) that there exists a subsequence of {C^|,}l>i (not indicated) with 
At = o(L-i) and a function Q G L'^{Q,T-K.)rM'^{0,T;H'^{n)') such that, as L ^ oo (and thereby 
At 0+) 

C^'i'^^'^Ce weak* in L°°(0,T;L°°(f})), (6.51a) 

Cfi+ ^ Ce weakly in L^{Q, T; iJ^(17)), (6.51b) 



weakly in L^{Q, T; H^{VL)'). (6.51c) 



dt dt 

Noting (j6.51aL c) and (I6.38c[) . we can pass to the limit as L oo in (|4.4I) to obtain that 

/ (%,^)ffi(0)d<+ / / [eV.C.-^ed -V.^dxdt^O V(^G L2(0,T;i/i(r!)). (6.52) 
Jq at Jq Jn ~ ~ J ~ 

Noting Fubini's theorem, point © of Lemma 16.21 see also (|6.28l) . and (13. 5p . we have for any 
(p G C°°(n) c L°°{n X D) that as L oo 

Ctl{x,0)f(x)Ax= [ ( f M{q)l3^{ip\x,q))dq]ipix)dx 
n ~ ~ ~ Jn \Jd ~ ~ ^ ^/ ~ ~ 

M{q) I3^{i>^{x, q)) Lp{x) dg dx 
nxD ~ ~ ~ ~ ~ ~ 

/ M{q)ii)oix,q))(p{x)dqdx 
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M{q) ipo{x, q)) dq ] ip{x) dx — / ip{x)dx. (6.53) 
n \J D ~ ~ ~ ~J ~ ~ Jn ~ ~ 

Hence we deduce from (|6.51a| -c) and (I6.53P that (^{x,0) ~ 1 for a.e. x E il. We note also from 
Fubini's theorem, and (|6.39dp that 



I I I Mi^,dq\dxdt^ f f \ f M{i;^,l-i^,)dq\dxdt 

Jo Jn JD ~ ~ Jq Jn Jd ~ ~ 

< [ [ M\ii^l-ii^\dqdxdt^O as L ^ oo. (6.54) 

Jo JqxD ' ~ ~ 

Hence (^1 fj^M^p^dq strongly in L^{0,T; L^{il)). Comparing this with (|6.51ap . we deduce 
that 



Ceix,t)^ M{q)ij,{x,q,t)dq for a.e. (a;, i) € X (0, T). (6.55) 

~ Jd ~ ~ ~ ~ 

Clearly the linear parabolic problem (|6.52p with initial data Ce(a:, 0) = 1 for a.e. x E ft has the 
unique solution = 1 on x [0,T]. This implies (|6.36p . and completes Step 3.11 and the proof. 

□ 

7 Exponential decay to the equilibrium solution 

We shall show that, in the absence of a body force (i.e. with / = 0), weak solutions (ue,^e) to 
(Pg), whose existence we have proved in the previous section, decay exponentially in time to the 
trivial solution of the steady counterpart of problem (P^) at a rate that is independent of the 
specific choice of the initial data for the Navier-Stokes and Fokker-Planck equations. Our result 
is similar to the one derived by Jourdain, Lelievre, Le Bris & Otto except that the arguments 
there were partially formal in the sense that the existence of a unique global-in-time solution, 
which was required to be regular enough, was assumed; in fact, the probability density function 
was supposed to be a classical solution to the Fokker-Planck equation; (cf. p. 105, therein; as well 
as the recent paper of Arnold, Carrillo & Manzini 0] for refinements and extensions). In contrast, 
we require no additional regularity hypotheses here. 

Theorem 7.1 Suppose that the assumptions of Theorem \6.1\ hold and M satisfies the Bakry- 
Emery condition (cf. Remark \5.1\) with k > 0; then, for any T > 0, 

ll«e(r)|p + A||^^(r)-i|||,^(^^^) 



|-ito|P + 2fc/ MT{ipo)dqdx 
JnxD ~ ~. 



where 70 := min , ^Iy'^ • In particular if f = 0, the following inequality holds: 

|mo||^ -l-2fc / MI-{i^a)dqdx 



+ - / 11/11^, d., (7.1) 



(7.2) 



Proof We take t ~ ti = At and write ~ tQ in (I4.24p . and we replace the function T on the 
left-hand side of (|4.24p by T^, noting that, prior to (|4.24p . in (|4.18p we in fact had on the 
left-hand side of the inequality, which was subsequently bounded below by J^; thus we reinstate 
the J"^ we previously had. We recall that m° = y^L~(^i) and I3^{ip°) ~ ''pf^'iti) and adopt the 
notational convention t-i :— —00 (say), which allows us to write jt^|^'^(to) instead of uf*£~{ti) 
and %l)f*£^{ta) instead of V'^^ (^i). Hence we have that 



At,-h/, Mi2 ( 1 /" \| At.+ At,-|i2 1 
teX (^l)ll + At / ""-^ " " 
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— a 



Jtn ~ ~ ' 



ao / Jto 



+ 2k f M (^f r (ti ) + a)dqdx+-^ r I M{^f'^ - ^f]:^-f dg dx ds 
JQxD ~ ~ ^T^^Jto JnxD ' ~ ~ 



/ / M^^— — dgdxds+ Yx 7A^+' dgd:rd5 



<\\uf^+{to)f + - f]\f^'^+{s)\\l,ds + 2k f MT'^{i^fl:+{to)+a)dqdx. (7.3) 



ti 

II 

to ~ JnxD 

Closer inspection of the procedure that resuhed in inequahty (|4.24p reveals that (j4.24p could have 
been, equivalently, arrived at by repeating the argument that gave us (j7.3p on each time interval 
[t„_i,i„], n = 1, ... ,7V; viz.. 



2k / M J-^ (i„) + a) dg dx + ^ / Af - tA^^^- )2 dg dx ds 

A^*.+ |2 , |Y7 „?Af,+ i2 



+ 2ke I / ■ dgdxd. + — / / M - dqdxds 



<||ufr(t„_i)|P + - /"||/^*^+(.s)||?.,d. 

+ 2fc /" AfJ-^(?Af['+(i„_i) +a)dgdx, n = l,...,iV, (7.4) 
JnxD ' ~ ~ 

summing these through n and then bounding on the left-hand side below by T . 

Here we proceed differently: we shall retain on both sides of (|7.4p . and omit the second, 
fifth and sixth term from the left-hand side of (|7.4I) . Thus we have that 

/"jiy..f:r(.)ipd. 

+ 2k I MF^{i)f^£+{tn) + a) dq dx 



nxD 

t 



2af)k 



I I M |V,7^^i'+ + a|2 dqdxds 
Jt„_i JnxD ~ V ' ^ 



<ii<i'+(i„_i)f +- rii/^*^+(.)fv"d5 
jt„_i ~ 

-h2fc /" A.fJ"^(7^fL + (t„_i) -|-a)d9dx, n = l,...,A^. (7.5) 

JnxD ' ~ ~ 

Thanks to Poincare's inequality, recall (|4.37p . there exists a positive constant Cp = Cp(il), such 
that 

\\uf^{;s)\\ <Cpin)\\V^,uf;£+{;s)\\ (7.6) 
for s e (t„_i, n = 1, . . . , iV. Also, by the logarithmic Sobolev inequality (|5.4p . we have for a.e. 
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X € Q that 

/ ^(g)[<r(;.g.^) + "]iog ^TTT^'^ 

< - / 

for s G = 1, . . . , -/V. Hence, for a.e. x € fi, 



2 

dq, 



{x, q, s) + a] dq 
1 

dq 



< - f M{q) 

K Jd 



^M{q)[ii^;^ + {x,q,s)+a]dqj log i^j ^M{q)[^^}+ {x, q, s) + a] dqj , (7.7) 

for s G {tn-i,tn], n — 1,...,N. Note that, thanks to (I4.31b[) and the monotonicity of the 
mapping s e M>o logs G M, the second factor in the second term on the right-hand side of 
(I7.7P is < log(l + a). Since a G (0, 1), we have log(l + a) > 0; also, the first factor in the second 
term on the right-hand side of (|7.7p is positive thanks to (I4.31ap and by ()4.31bp it is bounded 
above by {1 + a). Hence the second term on the right-hand side of (|7.7p is bounded above by the 
product {1 + a) log(l -I- a). We integrate the resulting inequality over V, to deduce that 

Miq)[$f*£^{x,q,s) + a]\og[ipf*j^'^{x,q,s) + a] dq dx 
nxD ~ ' ~ ~ ~ ^ ^ ~ 



< - f M{q)\V l^f'£+ix, q, s) + dq dx + \n\ (1 + a) log(l + a), 
JnxD ~ - V ~ ~ ~ ~ 

for s G n = 1,...^N. Equivalently, on noting that slogs = J'{s) ^ (1 ~ s), we can 

rewrite the last inequality in the following form: 



nxD 



M{q)Tiiif^'£+ix,q,.s) +a)dqdx 



2 

< - 

K 



Miq)\V H,^'£^{x,q,s) + a\' dq dx 
nxD ~ ~ V ' ~ ~ ~ ~ 

-I- [ M{q){l-^^f*£+{x,q,s)-a)dqdx+\n\{l + a)\og{l + a), (7.8) 

for s G {tn-i, tn], n = 1, . . . , N . This then in turn implies, thanks to the fact that q, •) is 

constant on (i„_i,i„] for all {x,q) GflxD, that 

^At I M{q)F{^f'^{t^) + a)dqdx 
^ JnxD ~ ~ ~ 

<^/*" / M|vVv?l'+ + «Pd<zdxds 



A 

K an fc 

Kar,k 

+^ 

for n= l,...,iV. 



i„_i JnxD 

t 



I I Af (1 - - a) dq dx ds 

Jtri^i JnxD ' ~ ~ 

At\n\{l + a)\og{l + a), (7.9) 
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Using (frB)) and ((7^ in (|73)) yields 



' JnxD ~ ~ 

+ 2k f M [J-^(^fl^^+(t„) + a) - ^(^f['+(i„) + a)] dg da: 
JnxD ~ ~ 

< ||wfl'+(i„-i)|P + 2fc / MT{ipfl:+{U,-i)+a)dqdx 
~ ' JnxD ~ ~ 

+ 2k f M [J-^(^^f[^+(t„_i) +a)- J-(V^f['+(t„_i) + a)] dq dx 
JnxD ' ~ ~ 



naok 
A 



M(l - ^^f - a) dg da;di 



e,L 

tn-i JnxD 



''"°^.Ai|r!|(l + a)log(l + a) + i T" ds, forn = 1, . . . , A^. (7.10) 



A 



We now introduce, for n = 1, . . . , N, the following notation: 

, , ■ f ^ f 2A A; Cm \ koq 
7(q!) :~ mm — ^ \ v — a 



Cp \ clq J 2A 

^„(a) := + 2fc / Af ['+(t„) + a) dq dx, 

~ ' JnxD ' ~ ~ 

Bnia) := 2k [ M [^^(^f ['+(t„) + a) - ['+(t„) + a)] dg da;, 
JnxD ' ' ~ ~ 



naak 



C„(a):=:^ / / M{l-^pf'£+{x,q,s)-a)dqdxds 
^ Jtn-i JnxD ' ~ ~ ~ ~ 

+ 1^At\n\{l + a)\og{l + a) + - r \\f^''+\\l,ds. 
^ ^ Jt„-i ~ 

We shall assume henceforth that a is sufficiently small in the sense that (|4.23p holds. For all such 
a, 7(ck) > 0; further, trivially, An{a) is nonnegative; by (14. lip , we have that -B„(a) is nonnegative, 
and by (j6.4p . C„(a) is also nonnegative. In terms of this notation (17.101) can be rewritten in the 
following compact form: 

(1 + 7(a) At) An{a) + B„(a) < A„_i(a) + B„_i(a) + C„(a), n = 1, . . . , iV. 

Equivalent ly, we can write this as follows: 

(l + 7(a)At)A„(a) < ^„__i(a) + D„(a), n^l,...,N, 

where D„{a) :— C„(a) — (-B„(a) — i3„_i(a)). It then follows by induction that 

n 

An{a) < (l+7(a)Ai)-'Mo(a)+^i^,(a), n=l,...,N. 

That is, 

A„(a) + B„(a) < (1 + 7(a) Ai)-'Mo(a) + | Bo{a) + Q(a) I , n^l,...,N. 
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In particular, with n = iV, by omitting the nonnegative term Biq{a) from the left-hand side of 
the resulting inequality, and recalling that T — — NAt, we have that 



\\ufl^{T)r + 2k I MF{i,^^l^{T)+a)Aqdx 



At. 



< 



2M1\ 

N J 



flxD 
N 



nxD ~ ~. 



+ 



2k 

KQok 

A 

naok 
A 



M 



nxD 



J-^(^,^i^+(0) +a)- J-(V7^r (0) + a) 



dq dx 



lo JnxD 



M(l — ^pf'i^~^{x, q, s) — a) dq dx ds 



T \n\ (1 + a) log(l + a) + - r\\f^''+fv, ds. 

Jo ~ 



(7.11) 



Using that 



y^l;+(0)=y° and ^^^1^+ = /3^(V^«), 
we then obtain from (I7.1ip that 



\\uf*l^{TW + 2k I AfJ-(^,-r(r)+a)dgdx 



nxD 



< 



N J 



-N 



2k 



MT{P^{^j°) +a)dq dx 



2k 



M 



flxD 



inxD 



dq dx 



KQok 

A 

Kao k 
A 



M(l — L~^{x, q, s) — a) dq dx di 



JnxD 



T\n\il + a)\og{l + a) + - r\\f^'^+\\l,ds. (7.12) 

Jo ~ 



Applying (14. lip and (I4.26P in the second factor in the first term on the right-hand side of (|7.12l) 
and applying (j4.25l) in the square brackets in the second term on the right-hand side, we have that 



\\uf^+iT)r + 2k I M^i^^'^^iT)+a)dqdx 



At,+ / 



nxD 



< 1 



N 

+ 3ak\n\ 

KGok 



-N 



,0||2 



3ak\n\+2k A/J"(V^" + a) dq dx 

JnxD ~ ~. 



inxD 



JnxD 



M{1 — tjj^ £ {x,q, s) — a) dq dx ds 



^ T\n\{l + a)\og{l + a) + - 
A f 



11/ 



(7.13) 



We now pass to the limit a — )■ 0+, with L and At fixed, in much the same way as in the previous 
section. Noting that limQ,_j.o+ tI") = 7o, we thus obtain from (|7.13p . p.20p and p. 221) . that 



uf]^iT)r + 2k I Mn^:i^{T))dqdx 



At. 



nxD 



< 



Kaok 



||uo|r + 2fc / MT{'ipo)dqdx 
nxD ~ ~. 



M{l-ip^*£+)dq dxds + i 



JnxD 



11/ 



(7.14) 
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Using (|6.38cp . (|OS)) . (j6.39d|) . (IPS)) and (13211), we pass to the limit with L ^ oo (whereby 
At — 0+ according to At — o{L^^) and therefore N — T/ At — >■ oo), to deduce from (|7.14[) that 



\\u,{T)Y + 2k I MF{iPe{T))Aqdx 

\uo\\^ + 2kf MJ^{^^o)dqdx]+- [ \\f\\^y,ds. (7.15) 



nxD 



Now, the Csiszar-Kullback inequality (cf., for example, (1.1) and (1.2) in the work of Un- 
terreiter, Arnold, Markowich & Toscani 40 ) with respect to the Gibbs measure /i defined by 
d/i = M{q) dq, yields, on noting (|6.36p , for a.e. x G il, that 



MF{i;,{x,q,T))dq 



which, after integration over 17 implies, by the Cauchy-Schwarz inequality and subsequent squaring 
of both sides, that 



^XT) - lllii (oxD) < 2|^^l / AfJ-(4(T))dqdx 

JnxD 



By combining this inequality with (|7.15l) . we deduce (|7.ip . Taking / = 0, the stated exponential 
decay in time of (u^, i/je) to (0, 1) in the L'^{^) x L\.j{V. x D) norm then follows from ()7.ip . □ 



Remark 7.1 By recalling our notational convention — Mip^, we see from (|7.2p that, in the 
absence of an external force, from any initial datum (uq, V'o) with initial velocity uq Cz H and initial 
probability density function tJjq that has finite relative entropy with respect to the Maxwellian M, 
the solution will, as i — )■ oo, evolve to the trivial solution (0, M) of the steady counterpart of the 
unsteady problem at an exponential rate that is independent of the choice of the initial datum. 
By introducing the free energy (the sum of the kinetic energy and the relative entropy): 

m-^lheit)f + k f MT{Mt))<iqdx, 

^ JnxD 
we deduce from (|7.15p that, for any T > 0, 

Jo ~ 

Thus in particular when / = 0, the free energy decays to as a function of time. o 
Remark 7.2 It is interesting to note the dependence of 

V KUo 

a 
p 



the rate at which the fluid relaxes to equilibrium, on the dimensionless viscosity coefficient v 
of the solvent, the minimum eigenvalue ag of the Rouse matrix A, the geometry of the flow 
domain encoded in the Poincare constant Cp(r2), the Weissenberg number A, and the Bakry- 
Emery constant k for the Maxwellian M of the model; as we have shown in Remark 15.11 in the 
case of a Hookean chain k = 1. We observe in particular that the right-hand side of the energy 
inequality (|6.42p and the rate 70 at which the fluid relaxes to equilibrium are independent of the 
centre-of-mass diffusion coefficient e appearing in ()1.10|) . o 
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Appendix A Density of C^{D) in Ll{D) and Hl{D) 

Suppose that w is a positive measurable function defined on D := 'R.^'^ such that w, G L°°{B) 
for any bounded open ball B C D centred at the origin of D. 

Consider, for j ■ = 1, 2, . . . , the sequence of concentric open balls Bj := B{0,j) in D. Consider 
further the sequence of functions {Cj}j>i C C^{D) such that < < 1 on D; Q = 1 on 
Bj] (j = on £) \ -Bj+i; \Cj\w^-^{Bj.f.i\Bj) ^ C, where C is a fixed positive constant, C > 1, 
independent of j. 

For V e L^{D), define vj := (j v. Clearly, 

ll^-^.lli^(D)= / (l-C,)'|^;|^^g)dg= / {1-Qr\v\'w{q)dq< [ \v\'w{q)dq. 

Jd ~ ~ JD\Bj ~ ~ JD\Bj ~ ~ 

Hence, given S > 0, there exists j >1 such that the right-most expression in this chain is bounded 
above by (5/2)^; therefore, also. 

Note further that vj € L^{Bj+i); this follows by observing that 

ll^^jlli2(B,.+i) = / \vjf dq < max w-\q) \vjfw{q)dq 

- '^'^^'^n l^l^«^(«)d«< max u;-i(g)||?;|||2(£,) <oo. 

As C^(i3j+i) is dense in L'^{Bj^i), there exists a function (pj € C^{Bj^i) such that 

II II ^ 

- V3\\L^Bi+^) < — — TT- 

2[msKq^Bj+^ w{q)\'^ 

On extending ^pj by from Sj+i to the whole of D and noting that ^pj G Co°(-D), we have that 

Vj — (pj\'^w{q) dq — / \vj — (pj\'^w{q) dq < max w{q) / \vj ~ ipj]"^ dq, 



from which we deduce that 

whereby, using the triangle inequality. 

Thus we have shown that for each 6 > there exists ipj G Cq°{D) such that ||w — <^j||i2^(£)) < d. 
This then implies that Cg°{D) is dense in L'^{D). 
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An analogous argument yields the density of C^{D) in Hl^{D). Given v G H^{D), we define 
Vj := Qv, with Cj as above, and note that, since -Bj+i \ Bj C D \ Bj and niax(f + 2C^, 2) < 3C^, 
we have 



W^-^jWhHD) < / \v\'w{q)dq + 2 \Vgv\'w{q)dq + 2C' K w{q) dq 

J D\Bj ~ ~ JD\Bj ~ ~ J Bj + i\Bj 

< 3C' [ {\v\' + \Y,v\')w{q)dq. 

J D\Bj 

Hence, given 5 > 0, there exists j >1 such that the right-most expression in this chain is bounded 
above by (5/2)^; therefore, also, 

h-'Vi\\Hi(D) < ^S. 

As Vj e iJg(i3j+i) and C^(i?j+i) is dense in iJg(_Bj+i), there exists a function £ C^{Bj+i) 
such that 

hj - V3\\m(B,+,) < — -— X- 

2 [maxggSj+i w{q)\2 

On extending (pj by from Bj+i to the whole of D and noting that (pj E C^{D), we have that 



\h - VjWhkd) = / - + lYqi-^j - '/'jOn w(g)dg < max wiq) \\vj - 

from which we deduce that 

whereby, using the triangle inequality, 

W^-'PjWhkd) < \\v - VjWhud) + \\vj - 'PjWhud) < + = S. 
This then implies that 0^(0) is dense in Hl{D). 
Remark Appendix A.l 

1. An argument identical to the one above shows that C^{D) is dense in H^{D) for any k > 0. 

2. In the special case when w is the Gaussian weight function, the density ofC^{D) in H^(D) 
follows by noting Definition 1.5.1 on p. 13 and Proposition 1.5.2 on p. 14 in Bogachev 



3. For the purposes of the present paper, the relevant choices of w are M and (1 + |g|)'^''M. 
4- Thanks to Theorem 8.10.2 on p. 4-18 in the monograph of Kufner, John & Fucik f26'^, 

LliiD), Ll,{nxD), Ll^^qj,,^,{D), Ll^^qj,,j,,{nxD), HIj{D), Hii{n x D) 

are separable Hilbert spaces for any ?9 > 0. 

Appendix B Fl,^(M^)nL2^^,^^,),,^,^(R^) ^ ^ > 1 

We begin by introducing our notational conventions. For r > 0, let 

D(^r)-^{q(^^'^ ■■ \q\<r} and D'-'''^ := {q e R"^ : \q\ > r} . 

Suppose that G is an unbounded domain in M'' such that GnL'(r) and GHD'-'''^ satisfy the segment 
property and the cone property (see 4.5 and 4.6 in Chapter 4, p. 82, of Adams & Fournier [T]). 
Let /i be a positive finite measure on G defined by d/i — w{q)dq, such that w S L°°{G) and w 

is locally bounded away from zero on G; i.e. for any compact set K C G, w{q) > 6k > a.e. on 
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K. We consider /i to be a measure on all of by defining //([/) = /x(C/ n G) for all Borel subsets 
U or R''. For 1 < p < cx) and any domain E contained in R'' we denote by Lf^iE) := LP{E; fi) 
the set of all (equivalence classes of /i-almost everywhere equal) functions defined on E whose pth 
power is Lebesgue-integrable on E with respect to the measure /z. 

Proposition Appendix B.l (Hooton |23] ) Let A4 be a subset of L^{G), 1 < p < oo, such 
that: 

(i) for each m G Z>o, the set of restrictions of functions in M to G D -D(m) is a relatively 
compact subset of LPj{G D 

(ii) given S > 0, there exists m e Z>o such that 

[ \^\Pdii<6 \f$eM. 

Then Ai is relatively compact in L^{G). 

The proof of this proposition is a simple modification of the proof of Theorem 2.33 in [I]. 

Theorem Appendix B.l For •§ > 1 the following compact embeddings hold: 

HM-i^"^) n L(^^|^^|)2iiMi(^'') '"^ ^fl+\q^\y^^Mii^'^)^ I = 1,. . . ,K. 

Proof Our argument follows closely the proof of Theorem 3.1 in the work of Hooton ^23". 

Since Mi{s) = Ci4 exp(— cas'') as s — +cx), with > 1, where qi, Ci4, i = 1,...,K, are 
positive constants whose actual values are of no relevance in the argument that follows, we shall, 
for the sake of clarity of exposition and ease of writing, set these constants to 1, and assume below 
that ^ 

AU{qJ :=e"(''2'l') , i^l,...,K. 

Suppose that if e (R'') n L^^^^^^ (R'^) , define ^{qi) :— [(p{qi)]'^, and consider the 
integral 

hAr,t):= f e"(^l2-l')'(l + |g,|)2^$L-iil ) dg„ 

where < t < r, d > 0, and 

D'-M := {g, e R'' : \q^\>r}, i^l,...,K. 
By performing the change of variable 




we deduce that {l — (t/\qi\)) \qi\ — \qi\, and therefore \qi\ = withO <t< rand \qi\ > r. We 

observe that in terms of the radial co-ordinate \qi\ and the d—1 angular co-ordinates (jiis, ■ ■ ■ , 4>i,d-i 
the d components of qi can be expressed as follows: 

qi,i = ki| cos(0i,i), 

qi,2 = \qi\ sin((/)j,i) cos((/>i^2), 

qi,3. = \qi\ sin(0j,i) sin(0j,2) cos((/)i^3), 

qi,d-i = \qi\ sin(0j,i) ••• sin(0j,d-2) cos((/)j,d-i), 
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qi,d = ISil sin(^i,i) ••• sin(^i,d_2) sin((^i,d_i), 
and thereby, letting Cd{4>i,i, • ■ • , '/'i.d-i) sin''^^(0i4) sin''~^(0i,2) • • • sm{4>i^d-2) , we have that 

dgi = \qi\'^~^Cd{(pi,i,...,(pi4-i) d\qi\ dcpi^i d(pi^2 ■■■ d(pi^d-i- 
With this notation, and observing that d\qi\ = d\qi\, we have that 

I ,-i d-l 

\qi\+t 



\qi\<^-^ 



dqi. 



Thus, 



= / 



{l + \q,\+tf^^qi) 



\qi\+t 



dqi. 



We deduce by bounding the integrand of this integral and renaming the dummy integration variable 
qi into qi that 

(^lg'l')''(l + |g,|)2''%,)dg,. 



where, independent of i, 



'){r,t) := sup 

qi££)i,(T-t) 



£,i,(r-t) 



l + \q^+t 



2i9 



Taking in particular t = 1 yields that 
hAr,^) = [ 



< 7(r,l) 



£)i,(r-l) 



(l + |g,|r$(q,)dg,. 



(B.l) 



Next, observe that 



/ 



(i + k.l) 



2i5 



< 



< 



Jo Jd*'( 
Jo iJd''( 



'o dt 

{^ + \q^\r 



dqi 



do/ qi 
V~ qi 



dt 



dqi 



ii + \qi\f\Y,Mqi)k dqi 



dt 
dt, 



(B.2) 



where | • j^^ denotes the l-norm on R'*, defined as the sum of the moduli of the components of the 
vector in whose norm is taken, and 



7i(r,t) := sup 



The Newton-Leibnitz formula implies that 



\qi\+t 



^i,i)<^[q^-Q^ 
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and we thus deduce from (jB.ll) and (IB.2[) that 



{l + \q,\r^^q,)dq 



<7(r,l) 



< S{r) 







(l + |g.ir<i>(gOdg, 

(l + |g,iny,,$(g,)ki dg, 



At 



(1 + 19.1)2" a>(Q.)dg 



where 



3 (^'2'''^ (l + |g,iny,.$(g.)kd<?, 



(5(r) := max 7(r, 1), / 7i(r,i)dtL 



(B.3) 



(B.4) 



We recall that $ := [i^J^; hence, Vg^^ = 2tf V qi^f- Use of Young's inequality 2ab < + under 
the second integral sign in (|B.3p with a = {1 + \qi\)''^\(p\ and b — \Yqi'?\ei yields that 



(1 + 19,1)2" [0{q,)fdq, 

(1 + 19,1)2" [^(90]'d9, 



\^q^viqt)\i, dg. 



(B.5) 



The inequality (|B.5|) is the analogue of inequality (8) on p. 575 in the work of Hooton [23] . 
Next, we show that if i? > 1 then limr_j.oo ^{r) = 0. Let us suppose to this end that < t < 1, 



■& > ^ and 2 < r; then. 



i\1^\+t? 



\q^r + 2M\q, 



|2i5-l 



Hence we deduce that 

0<7(r,t) <e-(^) " '5*('-*)^''-^2^-M |j = (2'^-i-2"32'') 
This implies that, for any i) > ^, 

lim 7(r, 1) = 0. 

r— J-oo 

Similarly, for < t < 1, 2 < r, we have that 



(B.6) 



< (2'^-i-''3''l 



sup e 

qieD^.^'—t') 



sup e 



(l + k^l+i)" 



We shall assume in the rest of this section that "i? > 1. 

For r > 2, consider the function G defined on [r — 1, oo) x [0, 1] by 



Gix,t) 



il + x + tf 



where x> r — 1, Q < t < 1, 
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with c := 2^ ^ d > and x :— \qi\. Our aim is to show that, in the hmit of r — > oo, 

/ sup G{x,t)dt 

Jo x>r~l 

converges to 0. As 71(7-, t) is positive for all r > 2 and t G [0, 1], this will then imply that 

lim / 7i(r,i)dt = 0. 

Clearly, for any t E [0, 1] and any r > 2 fixed, x E [r — l,oo) !—> G{x,t) e (0, 00) is a C°° 
function, and, for any t E (0, 1] we have that limaj^+oo G{x, t) = 0. Thus, for t E (0, 1] the function 
X I— ?> G{x,t) either attains its maximum value over the interval x E [r — l,oo) at the left-hand 
endpoint x = r — 1 of the interval [r — l,oo), or at an internal point xo{t) E (r — l,c») whose 
location is potentially i-dependent. 

In the former case, G(r - 1, t) := e"^* ('^"i' (r + tf. If on the other hand the maximum of 
the function x 1— G{x, t) is attained at an internal point xoit) E (r — 1, 00), then xo{t) annihilates 
the x-partial derivative of G: 

G^{x,t) = -ct{2d - 1) x2''-2e-^*"''"'(l + x + tY + i^e-^*""'"^! + x + tf-^ 

Now, Gx{x, i) = if, and only if, 

„2i?-2 - ^ 



tx'^-'{l + x + t) 



c(2i?- 1) 



As, by hypothesis, > 1, it follows that > 0; on the other hand, for each t E (0, 1], the 

function x 1— )■ tx'^^~'^ (1 + 2; + is a strictly monotonic increasing bijection of [0,oo) onto itself. 
Hence, for each t E (0, 1] there exists a unique positive real number xo(t) such that 

t [X„it)r-' (1 + [.TO W] +t)= ^^^J_^y (B.7) 

i.e. Gx{xo{t),t) = 0, which means that xo(i) E (0, 00) is the unique point of maximum of the 
function x i-> G{x,t) in the interval (0,oo), with t E (0, 1]. We emphasize that, at this stage, it is 
not yet clear whether or not, for a fixed value of r > 2, xo{t) E (r — l,c»). Differentiating both 
sides of (IB.7|) with respect to t e (0, 1], we deduce that 

xo{t){l+xo{t)+2t) 

Xo(t) - 



2t {d -!)(! + xo{t) +t) + txo{t) ■ 



As xo{t) > for all t E (0, 1], and > 1 by hypothesis, it follows that x'Q{t) < for aU t E (0, 1], 
and limt_>.o^ a;o(^) = +00. We define :— 2:0(1); as we have noted above a:o(l) > for all ■& > 1; 
and therefore, c^i > 0. In summary then, XQ{t) E [c^, 00) for all t E (0, 1] and all 1? > 1. 

For any r > O) + 1, there exists <r-i G (0,1] such that a;o(tr-i) = r — 1. Hence, for any 
r>ca + l, 

/ sup G{x,t)dt^ / sup G{x,t)dt+ / sup G{x,t)dt. (B.8) 

Jo x>r — l Jo x>r—l Jt^-i x>r — l 

Now, for t E [tr-i, 1], we have that xo{t) E [a;o(l), a;o(tr-i)] — [ci?,f — 1], and so the point xo{t) 
at which the maximum value of the function x 1— > G(a;, t) is attained over the interval (0, 00) does 
not belong to the open interval {r — \, 00). Hence, 

^ sup G{x,t)dt= [ G{r-l,t)dt< [ e-'=*('-i)"'~'(r + l)''dt < (^+1)'' 



x>r—l 
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Thanks to our assumption that > 1, it fohows that the right-hand side converges to as r ^ oo, 
and therefore, 



hm / sup G{x, t) dt = 0. 

Jtr-i x>r-l 



(B.9) 



On the other hand for t € (0,tr-i], the (unique) stationary point xa{t) of the positive function 
X I— >■ G{x,t) over the interval (0,oo) belongs to the open interval (r — l,oo), and the maximum 
value of a; i-> G{x,t) is then equal to G{xo{t),t). Hence, 



sup G{x, t) dt < 

x>r—l 



G{xQ{t),t) dt. 



(B.IO) 



It remains to understand the behaviour of 



/ G{xoit),t)dt 
Jo 



in the limit of r — >■ oo. We begin by noting that, as r — oo, necessarily i^-i ^ 0+ (cf. the 
definition of i^-i above, and recall that xq is a strictly monotonic decreasing function of t on 
(0, 1], with limt^o+ xo(t) = +00). Extracting xo{t) from the final factor on the left-hand side of 
(|B.7p and passing to the limit t — >■ 0+, it follows that 



Equivalently, 



Thus, 



lim t[xo{t)] 

t — ^U-j- 



lim t^^ xo(t) 



2i5-l 



?9 



c(2t9- 1)' 



c{2d- 1) 



lim 1 20-1 (^(xnft) ^) = e 20- 

t-)-0+ 



c(2i?- 1) 

This implies the existence of a positive constant C(i9) such that 

< G(a;o(i),t) < C(i?)i-^ Vie (0,1]. 
As > 1, and therefore 2^-1 < 1, we then deduce that 



lim 



G{xo{t),t)dt= lim / G(a;o(s),s)ds = 0. 



+ Jo 



Returning with this information to (IB.lOp . we have that 

lim / sup G{x,t)dt = (). 
Using ((BJ)) and ([RTT]) in (lR8l) implies that 

lim / sup G{x, t) dt — 0. 

''^ooJo x>r-l 

Hence, also, 

lim / 7i(r,i)dt = 0. 
Thus we deduce from (|B.6P and (|B.12p and the definition (|B.4p of (5(r) that 

I? > 1 =^ lim (5(r) = 0. 



(B.ll) 



(B.12) 



(B.13) 
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Having established (IB.5|) and (jB.13p . we suppose that is a closed and bounded set in the 
normed linear space iJ^^.(R'*) n L^j^^i^^ |^2i)^,/. (R'*), the norm of the space being 



|2 



Iffi (R^) 



The boundedness of M. then means that there exists Ct > such that 



for all ^ e . 

One can now easily show that satisfies conditions (i) and (ii) of Proposition [Appendix B.l 
with G = and p = 2. Indeed, (i) is a direct consequence of the Rellich-Kondrachov theorem 
(cf. Adams & Fournier T, p. 168, Theorem 6.3, Part I, eq. (3)), applied on a bounded ball Di {jn) 
of radius m centred at the origin, thanks to the fact that on Di ^jj^^ the weight functions Mi 
and (1 + \qi\f''^ Mi are bounded above and below by positive constants, whereby Hl.j{Di (,n)) and 
^fi+lq.iy^Ahi^^-X-'n)) coincide with H^{Di^i^rn)) and L2(A,(m)), respectively. 

To verify condition (ii) of Proposition [Appendix B.l[ we take S > and recall the definition 
(|B.4p of S{r). Thanks to (jB.13p . there exists a positive integer m such that 2C^S{m) < S. It then 
follows from (IB.5|) that 



[ e-(^lg'l')'(l + |(7.|)'''[^(9.)]'dg,<5. 



z = l,. 



Thus we have verified condition (ii) of Proposition [Appendix BTT| 



it follows that M. is compact in L?, , ,„ |>2d (K'^)- This completes the proof of the compact 



Therefore, as M is closed in Lj^^^^ ^^2i> ) , M. is relatively compact in |^2flj(_f 

2 ~ 

d 



embedding 
for ?? > 1 and i ^ I, . . . ,K. □ 



In the next section, we extend this argument to the case of K coupled dumbbells, K > 1. 



Appendix C Hi 



M\ 



i^'^) n L^^^|^|^2fl Af (^^'^) -^(i+|g|)2^M(^^'^) ^^>^ 



Theorem Appendix C.l For > 1, the following compact embedding holds: 



Hj^,j{R^ ) n £(i+lgl)2tfj,/(M ) -^(i+igD^f AfV 



Proof We proceed in the same way as in the previous section. 

Once again, since qi exp(— Cis'') < Mi{s) < Ci2 exp(— Cis'') as s — )■ +00, with d > 1, where 
Cii, Ci2, i = 1, ■ ■ ■ , K, and Ci are positive constants whose actual values are of no relevance in the 
argument that follows, we shall, for the sake of clarity of exposition and ease of writing, set these 
constants to 1, and assume below that 



M,(g,):=c (^Ig'l')', 
where q = {qi, . . .^qx)^. 



« = 1, 



and 



M{q) ~ Mi{qi)---MK{qK), 



75 



Suppose that (p € Hl^{W^^) n L^^^|^|j2^;i^(K^''), define := [^(g)]^, and consider tlie 

integral 



I.»{r,t) := / exp 

/_D(r-) 

wiiere < t < r, d > 0, and 



:= {g g : > r}. 
By performing the change of variable 



we deduce that (1 — (t/kl)) \q\ = I^L ^^id therefore |g| = \q\ + t, with < t < r and > r. 
Also, \qi\ = (1 — (i/l?!)) = {\q\/{\q\ + t)) \qi\ tor i = 1, . . . ,K. Analogously as in the previous 
section, 



dq 



\Q\+t 



K{d-1) 



dq. 



Thus, 



I^{r,t) 



D('--t) 



exp 



+ t 



2« 



K{d-1) 



dq. 



We deduce by bounding the integrand of this integral and renaming the dummy integration variable 
q into q that 



h{r,t)<^{r,t) J^^^ ^^exp -Y.i^-^\qJ^ j 



{l+\q\f'>^q) dq, 



where. 



j{r,t) := sup 



Taking in particular t = 1 yields that 

hir,!) < lirA) / exp 
Next, observe that 



exp 




'\q\+t' 


2t5' 


exp 


-Ef=i(^l?iP)'' 





1 + 



2i? 



Igl+^ 



if(d-l) 



(l + |g|r$((7)dg. 



(C.l) 



Lexp[-E(^I..P) 



exp 



/DM 



(i + kl)^" 

1 01.0 



io \q\ 



dq 



2i? 



do/ 9 



dt 



kl 



d(7 ^ dt 



(l + |g|/|y,%)kd5i dt, (C.2) 
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where 



7i(r,t):= sup 



The Newton-Leibnitz formula implies that 
and we thus deduce from (jC.l|) and (jC.2[) that 



exp 




'\q\+t' 




exp 




1 





(l + lgl+<)^^ 



igi+^' 



1 d 



dt 



kl 



^ (l + |gir$(g)dq 



<7(r,l) / e^-'-V^'-.- ; il + \q\r''^q)dq 



< 5{r) 



{l + \q\f\V,^{q)W dg 



dt 



(l + |g|r<i>(g)dg 



(l + |g|ny,$((j)k dg 



where 



(5(r) := max-^ 7(r, 1), / 7i(r,i)dt 



(C.3) 



(C.4) 



We recall that $ := [^]^; hence, — 2(p ^ q(p. Use of Young's inequality 2ah < + 6^ under 
the second integral sign in (jC.3p with a = (1 + |q|)''|^| and 6 = lYgi^lfi yields that 



DC-) 

< 2(5(r) 



*' ) {l + \q\Y%${q)Ydq 



(C.5) 



The inequality (IC.SP is the analogue of inequality (8) on p. 575 in the work of Hooton 
Next, we show that if i? > 1 then lim.r^oo 6{r) — 0. Let us suppose to this end that Q < t < 1, 
I? > i and 2 < r; then, 



Also, 



K 



5:ig.r> n.ax k.r = 

^ — ' ^ l<i<K 
i=l 

Hence we deduce that 



max Iffj 

l<i<K ~ 



> 



\qr + 2m\q\ 



1 ^ 



i5l„l2i9 
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This implies that, for any > i, 

lim 7(r,l) = 0. (C.6) 

r— )-oo 

Similarly, forO<t< 1, 2<r, and assuming that 'd > we have that 

l^ir,t) < (2'<i^-^^-'3') sup e-''^' ''^^^''{1 + k\ + tf 



< 



Arguing in an identical manner as in the previous section, we then deduce that if > 1 then 

lim / 7i(r,t)dt = 0. (C.7) 



Thus we deduce from (|C.6|) and (|C.7I) and the definition (IC.4P of 8[r) that 

t9 > 1 =^ lim (5(r) 0. (C.8) 

Having established (jC.5[) and (|C.8I) . we suppose that is a closed and bounded set in the 
normed linear space i/|^(R^'*) n Lj-|^^|^|^2^^^^(IR-'^''). The boundedness of means that there 

exists C| > such that 

/ e-^-^(^lg'l')''(l + |g|)2''[^(9)]2dg+ / e" (^'g'l')" |y,^(g)|,\ dgl < Q 
for all g . 

As in the case of i^T = 1, one can now easily show that satisfies conditions (i) and (ii) of 
Proposition [Appendix B.l with G = R^'* and p = 2 (g = 2 in the notation of Hooton) . Indeed, 
(i) is a direct consequence of the Rellich-Kondrachov theorem (cf. Adams & Fournier [IJ, p. 168, 
Theorem 6.3, Part I, eq. (3)), applied on a bounded ball of radius m centred at the origin, 

thanks to the fact that on the weight functions M and (1 + jgD^'^Af are bounded above and 
below by positive constants, whereby H\^(I)(^-^) and i^-|^^|^|^2i)ji,j(£'(m)) coincide with i?^(i?(,„)) 

and L'^{D(^„i)), respectively. 

To verify condition (ii), we take 6 > 0. Noting (|C.4|) and <\C.8\i . there exists a positive integer 
m such that 2C|(5(m) < S. It then follows from (jC.SP that 

/ e-^-^(^lg'l^)\l + k|f^[^(<z)]^dg<<5. 
Thus we have verified condition (ii) of Proposition [Appendix B.l[ 

-2 fnKd 



in L^j^^|g|-|2tf ^/(R^'''). This completes the proof of the compact embedding 



Consequently, A4 is closed and relatively compact in i^j^^|^|^20j;^(R ), and thereby compact 

2 ( 

(i+lgD^tfAf V 

HIA^^'^) n i(i+|g|)2<)ji_f(R''^'^) L^^_|_|q|j2flAf(l^''^'') 
for §>!. □ 

Remark Appendix C.l When i? = 1, i/ie Maxwellian M becomes a normalized Gaussian. It 
can then be shown using Gross' logarithmic Soholev inequality 121 1 and the logarithmic Fenchel- 
Young inequality that 

Hence in particular H^.^IR^'^) is contained in -^(i+|g|)2j(j (R^''). The quest ion arises whether this 

inclusion is in fact a compact embedding. We suspect that the answer is negative. On the other 
hand, one can still show that iJ^j(R^'') ^->-> L\^(M.-^'^) for alii) > i. This we shall do in the next 
section. 
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Appendix D HIj{D) ^ LIj{D), > i 

Let £> Di X • • • X Dk, where D, ^ R'^ , i ^ 1, . . . , K , and write M{q) := Mi(gi) • • • MxigK)- 
We shall prove that, for i9 > i, 

Hl{D) ^ LIj{D). (D.l) 

We begin by noting that from Theorem 3.1 in the work of Hooton [53] and a slight modification of 
Example on p. 5 76 therein — the only difference being that there the weight-function is exp(— jg^p'') 

whereas here it is equivalent to exp (^|<ZiP)''^ — we have that H^{Di\ /i^) is compactly embed- 
ded in L^{Di; jii) for > \, where is the measure on Di defined by d/i^ := cxp dg^. 
Hence, by (|1.4ap . we have that H\j\Di) is compactly embedded in L\j,{Di)] i.e., 

HliSD,)'-^ Llj^{Di), t = l,...,K; z9 > i. (D.2) 

The proof of the compactness of the embedding of H\j{D) into L\,j{D) then proceeds identically 
as in [10], using Theorem 2.4 in the work of Opic [M]. We shall prove (jP.ip for the case oi K ~ 2, 
with D = Di X D2 and M{q) = Mi{qi)M2{q2) ■ For K > 2 the proof is completely analogous. 

Let u e H\.j{D). As M = Mi x M2, it follows from Fubini's theorem that, for almost all 
qi e -Di, 

w(gi,-)eiLp2) and d'^u{qi,-)(^LlijD2), 

where a is any d-component multi-index with < |a| < 1. Fubini's theorem also implies that, 
given ip2 G C^{D2) and a d-component multi-index ck2, with < |a2| < 1, we have 



(-1) / 7/(qi,-)5"V2d<Z2 

D2 ~ 



tfi dqi. 



J Di UD2 

for all ipi e C§°{Di). Therefore, 9"^ [u{qi, •)] = d^°-°'^'^u{qi, ■) in the sense of weak derivatives on 

D2 for almost all qi G Di. As 9*^°'"^^w((7i, •) belongs to L\j^{D2) for almost all qi G I?i we have 
that 

u{qi, ■) e HIj^{D2) for almost all qi G Di. (D.3) 

Analogously, 

u{-,q2) G H\,j^{Di) for almost all 92 G -D2- 

As each of the partial Maxwellians, Mi and M2, is bounded from above and below by positive 
constants on compact subsets of their respective domains, there exists a sequence (I?^ („) : n G N) 
of open proper Lipschitz subsets of Di, i = 1,2, such that 

00 

A,(n) C A,(n+i), ?^ G N, IJ A,(n) = A and (A,(«)) ^1/; (A,(«)); 

e.g., A,(ti) = ^(Oi"); the compact embeddings stated here follow by the Rellich-Kondrachov 
theorem (cf. Adams & Fournier [1], p. 168, Theorem 6.3, Part I, eq. (3)) applied on A,(n)i 

.2 



i = 1,2, n G N. Letting, for n G N, (^X j^j^ A,(n)j £ -D, and noting that, by Appendix A 

in [10] ■ i?(„) is a Lipschitz domain, the above properties get inherited by D(n) from A,(n)- 

00 

i:'(„) C I?(„+i), riG N, y and i/i^ (£>(„)) (£>(„)). 

n=l 

Let D^^"^ \ D,.(„) and := \ £>(„)• It follows from Opic [34], Theorem 2.4, that the 

above compact embeddings on members of a nested covering imply the following characterizations 
(the first, for i G {1,2}): 

H\u{D,) ^ Lli^iDi) ^ lim sup / A^i, / \\u\\]j, ^ . - 0, (D.4) 
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Hl,{D)^Ll,iD)^ lim sup / d/i/h||k(D) = 0, (D.5) 

where dfii :— Mi{qi) dqi, i = 1,2, and d/i :~ M{q)dq. By virtue of (|D.2p . the left-hand side of 
(jD.4p holds; hence, its right-hand side also holds. Using (jD.3|) and (|D.4p with i = 2, we deduce 
that for any 5 > there exists n — n{6) ^ N such that 



dii = 



Di 



(") 



u'^iqi, ■)dfi 



2 



u{qi,-)d^2+ / |Vq2u(gi,-)l d^2 
r>2 ~ JD2 



d/ii 



and similarly for /^(i)^^,^ d/Lt. Then, as D^"-^ = {Di x Dj"^) '-^ (-D^"^ x D2), the right-hand side 
of (jP.Sp holds; therefore, so does its left-hand side; hence (|D.1|) . 

Appendix E Compact embeddings of Banach-space valued 
Sobolev spaces 

We begin by proving the following theorem, which can be seen as a generalization of the classical 
Kolmogorov-Riesz theorem to the case of Bochner spaces L'^{M.'^; E), 1 < r < 00. 

Theorem Appendix E.l (Kolmogorov— Riesz) Let 1 < r < 00 and suppose that Eq and Ei 
are Banach spaces, with Eq Ei . Suppose that T is a su bset of L^R'^;Eq), such that: 

(i) T is hounded in L^{W'-;Eq); 

(ii) for every (5 > there exists an R > such that, for every / G J-", 

X\>R 

(Hi) for every S > there exists p > such that, for every f ^ T and every y G with \y\ < p, 

ll/fe + y)-/(3;)llk dg;<<5'-. 
Then, J- is totally bounded in U'{W'';Ei). 

Proof Suppose that F C U'iW''; Eq) satisfies the three conditions of the theorem. For any 5' > 0, 
choose ^ > such that (2"^ + 1)~5 < ^S' . Then, for such a. S > 0, choose i? > as in the second 
condition, and p as in the third condition. 

Suppose that Q is an open cube centred at the origin of such that |y| < for all ?/ G Q. 
Let further Qi, ■ ■ ■ ,Qn be mutually nonoverlapping translates of Q such that 



contains a ball with radius R centred at the origin of M."^. For any / G i'"(M''; £'0), we define 

I 0, otherwise. 
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Hence, by Jensen's inequality, 

~ " I 0, otherwise, 

which, in turn, implies that H^^q J||l'-(q,:£;o) < II/IIl'-(q,;_Eo)- Let us define, for / e L''{R'^;Eo) 

N 



k=l 



Then, 



N 

-P/llL'-(R<i;_Eo) ^ \\PQkf\ 



k=l 

N 

k=l 

Thus, the mapping P : f ^ L^{R'^;Eq) h- Pf e L'^{R'^;Eo) is a bounded linear operator, whose 
norm is equal to 1. 

Let us take any i e {1, . . . , N} and consider the set Ai := {PQifix) : f & x £ Qi}. Noting 
that for / £ L^{R'^; Eq) fixed the value of PQif{x) is independent of a; G Q,;, and using (jE.2p . we 
have that 

N 

WPqJWeo m - WPqJWlw.) <T.\\PqJ\\Uq.;e.) < WfWlw^.) £ {1,...,N}. 

k=l 

(E.3) 

As J- is, by hypothesis, a bounded set in L^{R'^;Eo), it follows from (|E.3|) that Ai is a bounded 
set in Eq. Since i?o ^-H- we then deduce that Ai is totally bounded in Ei. Thus, for each 
6' > 0, there exist a positive integer Mi = Mi{5') and balls Aij, j = 1, . . . ,Mi{6'), in Ei whose 
centres Cij, j = 1, . . . , Mi{S'), are elements of Ai, whose radii (in the norm of Ei) are smaller than 
6'/{2{N\Q\)-) for each j = 1, . . -,^^{6'), and whose union covers Ai, in the sense that 

: / e ^, x e go c U Aj. (E.4) 

In other words, for every each i £ {1, . . . ,N} and for each Pq^J , where f £ T, there exists a 
j = j{f,i) £ {1, • . • ,Mi{S')} and Cij(/,i) G A, and a ball in Ei of radius (57(2(A^|g|)^) centred 
at Cij(^f^i), such that 

Let XQi denote the characteristic function of the set Qi, i = 1, . . . , A^. Hence, by the definition 
(lElj) of P, 

\\Pf-Y.^Q^^^^U,^)\\L'■iV^^■,E^) < ■ (E.5) 
i=\ 

Noting again that Pf{x) — PqJ{x) for all a; e gj, i = 1, . . . , A^, it follows from (|E.4p that 

N N M,(S') 

{Pf{x) : f £F,x£\Jq,}^\J U a,,. 

i—1 2—1 J — 1 
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For any x contained in the complement of [j^^i Qi, we then have, by the definition of P, that 
Pf{x) = 0; therefore, 

N M,{S') 

: / e X e M"^} c {0} u U U A,,cEi. 

i=l 3=1 

Consider the (finite) set 

:={0}u|^XQ,A.(/) : A,(/) C,,(/,,) , i = 1, . . . , TV, f e C E, 

of cardinahty < 1 + J^fLi M,{S'). 

Then, for f G JT, we have from (jE.SP that 



N 



N 



i=l 



i=l 



< I1/-^'/IIl^(k-b,) + ^'^'- 



(E.6) 



It remains to bound the first term on the right-hand side or (jE.6[) . 

It foUows from property (ii) and the definition of P that, for any / € J^, we have 

JV 

\\f-Pf\\lw,) < '^'■ + E/ ll/fe)-^/(s)llkdS 



N 



dx. 



By applying Jensen's inequality again, performing a change of variable, and noting that x— z (z 2Q 
when X, 2 € Qi, we have that 



\f - Pf\\ln^^;E,) < S'-+Y,J^ J^J^\\f{x)-,f{zWE^dzdx 



< 5'' + 



1=1 
1 

W\ 



2Q JR'' 



l/(s)-/fe + y)llB, dxdy 



< + I '5''dy = (2'^ + l)<5^ 



thanks to property (iii). Thus, by our choice of 6 for a given 6' at the start of the proof, 

1I/-P/IIl-(r-£0 < (2'^ + l)^<5< 



(E.7) 



By inserting (IE.7P into (jE.6p we deduce that for any 5' > Q and any / G there exists an element 
of the finite set £5', which in the norm of U'{W^; Ei) is at a distance < 5' from /. Hence £5' is a 
finite 5' -net for J" in Ei), whereby T is totally bounded in L''(M''; £'1). □ 

In a complete metric space a set is totally bounded if, and only if, it is relatively com- 
pact. Since L^(M''; E) is a Banach space, the words totally bounded in Theorem [Appendix E.l| are 
synonymous to the words relatively compact. 
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In the sequel, if necessary, we shall consider any function / defined almost everywhere on an 
open set $7 C M'^ to be extended by zero to the whole of M''; i.e., we shall introduce the function 
/ defined for almost all a; G M'' by 



fix) := Xn{x)f{x) = 



/(s) if « e 17, 
otherwise. 



Instead of writing f{x) we shall often simply write f{x) also when x ^ f2. 

Theorem Appendix E.2 Let il be a bounded open set in M'', 1 < r < oo, and let Eq and Ei be 

Banach spaces, with Ea ^-H- Ei . Suppose that 

(i) J- is a bounded subset of L^{fl;EQ); 
(a) For every 5 > there exists p > such that 

for each f E J- and all y G M'' with \y\ < p. 
Then, T is relatively compact in £''(i7; E\). 

Proof The set xo.^ satisfies the conditions of Theorem [Appendix E.l[ condition (ii) of Theorem 
[Appendix E.l being satisfied trivially since f2 is assumed to be a bounded subset of here. □ 

Theorem Appendix E.3 Suppose that VI is a bounded open Lipschitz domain in Mf^ . Let Eq 
and El be Banach spaces, with Eq ^-h- Ei . Suppose further that 1 < p < oo and J- is a bounded 
subset of 

W^'P{VL- Eq, El) {v e LP{n; Eq) : D^v € LP{fl; Ei) Va G N^to such that \a\ = 1}. 

Here \a\ := |g|f ^ = X^iLi '^i the length of the multi-index a. Then, J- is relatively compact in 
L^{n; El) for all r such that 1 < r < oo and 7 > ^ ~ ^• 

Proof In a metric space a set is relatively compact if, and only if, it is sequentially relatively 
compact. It therefore suflices to prove that any bounded sequence C W^'P{Q; Eq, Ei), 

i.e., such sequence that there exists a cq > for which 



S^P\\fn\\w^.Pin;Eo,Er) ^ SUp | WfnW'l^^n-.Eo) + J2 W^" f^\\lp{n;E,) ] ^o, (E.8 

\a\=l 



neN riGN 



forms a relatively compact subset of L^ (fl; Ei) , with r as in the statement of the theorem. 

We shall first suppose that I < p < d. Then, by the Sobolev embedding theorem, we have that 
W^'P{n;EQ,Ei) W^'P{n;Ei) ^ L''' [0.; Ei) , where r* := clearly, r* > 1. Thus, 

sup||/„||Lr*(0;£;,) := ci < cx). (E.9) 

neN 

As p > 1, it follows from (jE.Sp and Holder's inequality that 

sup \\fn\\w^:Hn;Eo,Ei) ■= C2 < I^^I^^^Cn < OO. 

Now, let (5 > be arbitrary but fixed, and let il* C il* C be a domain such that 
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Let p > be such that x & ft* and y € M.'^, \y\ < p, implies that x + y£Q. Assume that p < 
Then, by Jensen's inequaUty and the triangle inequality, 



\fnix + y) - fn{x)\\Ei 



^ dx. 



{x + ty)yi At 







i=l 



< 



El 



1 d 

E 

1=1 



dfn 

dxi 



\yz\dt, 



El 



which then implies that 



\\fn{x + y) - fn{x)\\Ei dx < \y\ \\ fn\\w^-^n;Ei} < P C2 < - , 



and therefore 

\\fn{x + y) ~ fn{x)\\Eidx < 



si\n* 



\fnix + y)\\Ei dx + 



n\n' 



\fn{x)\\Ei dx 



+ / \\Mx + y)- fn{x)\\Eidx<S. 

Hence, the assertion of the theorem for r = 1 follows from Theorem [Appendix E.2| Thus we have 
shown that the bounded sequence {fn}?^Li in W^'P{^l;Eo,Ei) has a subsequence {fn,,}'kLi that 
is convergent in L^{n;Ei). 

Let 1 < r <r*. Then, noting that < ^-r^j < 1, using Holder's inequality and (|E.9p . we have 
that 



Hence {/rifcjfc^i is a Cauchy sequence in Ei); since the latter is a Banach space, the sequence 

{/".ib}fe°=i is convergent in L^{il;Ei). Thus, combining the outcomes of the cases r — 1 and 1 < 
r < r*, we have shown that the bounded sequence {/n}^i in W^'^{il.; Eq, Ei) has a subsequence 
{/rifclfc^i that is convergent in L^{fl; Ei) for 1 < r < r* , whenever 1 < p < d. 

We have thus shown that a bounded set J- in W^'^ifl; Eq, Ei) is relatively compact in L^{fl; Ei) 
for 1 < r < r* , whenever 1 < p < d. 

Suppose now that p> d, and is a bounded set in W^''''{il; Eq, Ei). As is a bounded set in 
M'^, the set T is then automatically a bounded set in W^''^{fl; Eq, Ei) for all s G [1, d) by Holder's 
inequality, and the stated result follows for all r g [l,oo) from what was proved above. □ 

Next, we shall extend Theorem [Appendix E.3| to the case of three Banach spaces, Eq, E and 
El, where Eq E ^ Ei and Eq is reflexive. We shall confine ourselves to considering the case 
when 1 < p < oo, as this range of p will suffice for our purposes here. 

Theorem Appendix E.4 Suppose that Eq, E and Ei are Banach spaces, with Eq ^->-> E ^ Ei 
and Eq reflexive, and let be a bounded open Lipschitz domain in R'^ . For 1 < p < oo, we define 



W^'P{n;EQ,Ei 
Here, again, |g| := 



{v e LP{Vl; Eq) : D^v € LP{VL; Ei) Va G N|o such that \a\ = 1}. 

J = X^iLi length of the multi-index a. Then, 

W^'P{n;EQ,Ei)^ L^in-E). 



Proof Suppose that J" is a bounded set in W^'P{Q; Eq, Ei). As Eq E ^ Ei, it follows 
that Eq El . We thus deduce from Theorem [Appendix E.3| that T is relatively compact in 
LP(0;i?i). In order to complete the proof, we need to show that T is in fact relatively compact 
in LP{U;E). Once again we shall rely in the proof on the fact that in a metric space relative 
compactness and sequential relative compactness are equivalent. 
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Suppose that {fn}n=i is a bounded sequence in W^'P{n; Eq, Ei). Since W^'P{n; Eq, Ei) C 
LP{n;EQ), the sequence {fn}^=i is bounded in LP(r2;£'o). As 1 < p < oo and Eq is reflexive, 
also LP{n;Eo) is reflexive. Therefore {fn}^=i has a weakly convergent subsequence {/ntlfc^i in 
LP(f7; i?o); we denote the weak limit by fif\f^,^ LP{il; Eq). 

Further, if follows from the discussion in the first paragraph of the proof that {fnk}kLi is rela- 
tively compact in L''{i^; Ei). We can therefore extract a subsubsequence {/„j.^ }iZi from {fnk}k^i: 
which strongly converges in £^(57; Ei) to a limit / e LP{n; Ei). It follows from the uniqueness of 
the limit that = /, and therefore / e LP{i^; Eo){c LP{n;E)), in fact. It remains to show that 
converges to / in LP{Q.]E). 

As Eq E ^ El we deduce from Ehrling's Lemma (cf. Temam 39 Ch. Ill, Lemma 2.1) 
that for every 5 > there exists cs > such that, for all v ^ Eq, 

\\v\\e < 5\\v\\eo+cs\\v\\e,. (E.IO) 

As / - /„,^ is known to belong to LP{n;Eo){c LP{n;E) C LP{n;Ei)) for all / > 1, we deduce 
from that 

II /(S) - (S)IU < SWfix) - + cs\\f{x) - {x)\\e„ for a.c. x e Q. 

Thus, by the triangle inequality in LP{D,), 

11/ - /nfc, \\LP{n;E) < S\\f - /nfc, IIlJ'(0;_Eo) + C^H/ - /„^^ \\LP{n:Ex)- (E.H) 

Since {/nt^ is a bounded sequence in LP{il; Eq) and / G LP{il; Eq), there exists a positive 
real number K such that 

\\f-fn,,\\LPin-,Eo)<K V/>1. 

Further, for 6 > fixed (but otherwise arbitrary) and therefore eg (as above) also fixed, there 
exists Iq = lo{S) such that 

11/ ^ /nfc, llLP(n:Bi) < — yi>lQ. 

Consequently, by substituting the last two inequalities into (jE.lip . we deduce that 

11/ - fn,, \\L.(n,E) <SK + S= {K + 1)5 yi > Iq. 
Thus we have shown that for each S > there exists = ^o(<5) such that 

\\f - fn,^\\LHn;E) < {K + l)S ^1 > Iq. 

In other words, {/n^,}^! converges to / strongly in LP{V,; E). □ 

Appendix F Hlj{n x D) ^ L|^(Q x D), 'd > ^ 

With the compact embedding of H\^{D) in L\j{D) established in Section [Appendix^ the proof 
of the compact embedding H\^{^1 x D) '-H- L\j{Q, x £)),-(? > i, proceeds, verbatim, as in [IJ]: first 
the isometric isomorphism of L\j{^) and L?{Q.;L\j{D)) is proved using the separability of the 
Hilbert space L\j{D); then the isometric isomorphism _ff°^^(il x D) and LF' [Vt\ H\j{D)) is shown, 
where 

Hl^{n X D) := {V e Lli{Vt x D) : V qV G Llj{Vt x £»)}, 
and the isometric isomorphism of H^f{Q, x D) and H^{Vl\L\j{D)) is also shown, where 

Hlf{n X £>) := {V e iM(17 X D) : VxV G ^^(ll x D)}. 

With these definitions, we then identify the space L\j{il x D) with L^(r2; L\j{D)) and the space 
Hl^{nxD) = Hl;j"{nxD)nH°^^{nxD) with H^in-, LlfiD))nL'^{n; HIi{D)). Upon doing so, the 
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compact embedding of Hlj{rt x D) into L\,j{il x D) directly follows from the compact embedding 
of H^{n;Llj{D))r\L'^{n;Hlj{D)) into L'^iO.; L\j{D)), implied by Theorem [Xjjpendix E.4| above, 
thanks to the compact embedding of Eq := H\,j{D) into E = Ei = L\j{D). 

London & Oxford, 16 August 2010. 
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